THE ANTINUCLEON-NUCLEON INTERACTION 
AT LOW ENERGY: ANNIHILATION DYNAMICS 



Eberhard Klempt 

Helmholtz-Institut fur Strahlen- und Kemphysik 
der Rheinische Friedrich-Wilhelms Universitat 
NuBallee 14-16, D-53115 Bonn, Germany 

Chris Batty 

Rutherford Appleton Laboratory 
Chilton, Didcot 
Oxon, 0X11 OQX, United Kingdom 

Jean-Marc Richard 

Laboratoire de Physique Subatomique et Cosmologie 
Universite Joseph Fourier-CNRS-IN2P3 
53, avenue des Martyrs, F-38026 Grenoble Cedex, France 

February 5, 2008 



Abstract 

The general properties of antiproton-proton annihilation at rest are presented, with special 
focus on the two-meson final states. The data exhibit remarkable dynamical selection rules: some 
allowed annihilation modes are suppressed by one order of magnitude with respect to modes of 
comparable phase-space. Various phenomenological analyses are reviewed, based on microscopic 
quark dynamics or symmetry considerations. The role of initial- and final-state interaction is also 
examined. 
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1 Introduction 

1.1 Annihilation in hadron physics 

Annihilation is a fascinating process, one of the most interesting in low-energy hadron physics, in 
which matter undergoes a transition from its baryon structure to one consisting solely of mesons. 
In the early days of antiproton physics, antinucleon-nucleon (NN) annihilation was considered by 
analogy with positronium annihilation in QED, and described as a short-range process mediated 
by baryon exchange. Nowadays the quark model offers a drastic alternative, where the so-called 
"annihilation" does not imply actual annihilation of all incoming quarks and antiquarks, but simply 
results from their rearrangement into quark-antiquark pairs. Were quark rearrangement to be the 
leading mechanism, NN annihilation would be better considered by analogy with rearrangement 
collisions in atomic or molecular physics. Intermediate scenarios are however conceivable, where 
some of the incoming quarks and antiquarks annihilate, and new quark-antiquark pairs are created. 

This review is part of a project devoted to strong interaction physics with low-energy antipro- 
tons, as measured at the LEAR facility of CERN. A first part [1] was devoted to NN scattering and 
to antiprotonic hydrogen and deuterium. The present review covers the general properties of anni- 
hilation, the results on two-meson final states and their phenomenological analysis. A third article 
will concentrate on meson spectroscopy, as studied from multimeson final states of annihilation. 

1.2 Historical considerations 

Detailed studies of antiproton-proton annihilation at rest were carried out in the 1960's, and the 
results are still significant for studies of annihilation dynamics and meson spectroscopy. These 
experiments were performed at the Brookhaven National Laboratory (BNL) and at CERN in Geneva 
by stopping antiprotons in bubble chambers. Analysis methods and early results were reviewed in 
detail by Armenteros and French [2], but many important results were not included. Later reviews 
[3-5] focused primarily on new concepts and developments and did not aim at a comprehensive 
description of all experimental data. The physics results on pN annihilation obtained from bubble 
chambers filled with H2 or D2 are included in our report. 

In the 70's, NN physics was dominated by claims for narrow baryonium states, which were 
not confirmed by more systematic searches. The motivation for quasi-nuclear NN states and for 
multiquark states preferentially coupled to NN and the experimental results have been reviewed 
extensively in several articles [6-10]. 

Research on pp annihilation was resumed in 1983 when LEAR came into operation. The Asterix 
collaboration investigated annihilation from P-states of the pp atom formed in H2 gas with a In 
electronic detector. The focus of the research was dynamical selection rules which will be discussed 
in some detail in Sec. 7. A broad resonance, called AX(1565), possibly a quasi-nuclear state, was 
discovered. The search for narrow states produced in annihilation at rest continued both at LEAR 
and KEK, eventually yielding negative answers. 

In more recent years, two 47r spectrometers. Crystal Barrel (PS197) and Obelix (PS201), took 
data on pp annihilation at LEAR. The Crystal-Barrel research activity was directed towards anni- 
hilation at rest and in flight. Obelix investigated antiproton and antineutron [11] interactions at rest 
and with very low momenta. Nuclear physics was also an important part of the Obelix program. 

The experimental progress was accompanied by active groups of theoreticians trying to under- 
stand the basic mechanisms responsible for annihilation. From a theoretical point of view annihi- 
lation is a very complicated process which is Ukely driven by both the underlying quark dynamics 
and by conventional hadronic interactions. If for instance, the NN potential is attractive in one par- 
tial wave, and repulsive in another, one expects annihilation from the former to be enhanced, and 
annihilation from the latter to be suppressed. Similar remarks hold for final state interactions with, 
in addition, the possibiUty of interferences between, for instance, primary p mesons formed by qq 
pairs and p mesons built from tttt final-state interactions. An accurate description for all annihi- 
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lation rates seems therefore to be unlikely. It is hoped, however, that the leading mechanisms of 
annihilation will be identified, in particular to explain the so-called dynamical selection rules, the 
observation that some annihilation modes are strongly suppressed in certain partial waves, while still 
being allowed by energy and quantum-number conservation. 

1.3 Outline 

This review begins, in Sec. 2, with a presentation of the beams and detector facilities used to measure 
annihilation properties. In Sec. 3, we briefly summarise the properties of the mesons seen in anni- 
hilation experiments. Kinematics and conservation laws are reviewed in Sec. 4. The main features 
of annihilation, as seen in various experiments, are presented in Sec. 5, while Sec. 6 is devoted to 
a thorough review of the rates into various two-meson final states. The dynamical selection rules 
are presented and discussed in Sec. 7. Section 8 contains a critical survey of various approaches 
to annihilation mechanisms, and an analysis of what is learned from the systematics of two-body 
branching ratios. Some conclusions are presented in Sec. 9. 

1.4 A guide to the literature 

The physics mediated by antiprotons has been presented at many Conferences, in particular the NAN 
conferences, the LEAR Workshops, and the LEAP conferences resulting from the merging of these 
two series, as well as at some Schools [12^0]. 

The early review by Armenteros and French [2] remains a reference for early annihilation data. 
Before the completion of LEAR measurements and the final analysis of the data, important review 
articles became available; some concerning general aspects of LEAR physics [3, 10], whilst oth- 
ers specialised more on the annihilation process [4, 5]. Antineutron physics, including antineutron 
annihilation, is reviewed in [1 IJ. A review devoted to annihilation in flight appeared recently [41]. 
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2 Beams and experiments 

2.1 Early experiments at CERN and BNL 

Following the discovery of the antiproton in 1955, p beams were rapidly developed and a first survey 
of antiproton annihilation on protons or neutrons was possible, as early as in the 1960's, by stopping 
antiprotons in hydrogen- and deuterium-filled bubble chambers. These experiments demonstrated 
that pp annihilation is a powerful tool to discover meson resonances, even though only limited 
statistics were achieved. Some of the early results are still important, and it seems appropriate to 
include a short discussion as to how the data were obtained. 

Two main experiments were carried out at that time: the first one at Brookhaven by a group 
from Columbia University and the other at CERN by a CERN-College de France collaboration. The 
experimental conditions were closely similar and it is sufficient to discuss just one of them. 

The CERN bubble chamber, built at Saclay, had an illuminated volume of 80 cm in length, and of 
30 cm in height and depth. Antiprotons from a separated antiproton beam of momentum 700 MeV/c 
were moderated in a Cu degrader and stopped in the target. The chamber was situated in a magnetic 
field of 2.1 T. Due to the momentum spread in the incident beam and multiple scattering in the 
degrader and target, the stopping distribution was rather wide. A cut on a minimum track length 
of 5 cm guaranteed a minimum momentum resolution; the average track length was 16 cm. The 
intensity of the antiproton beam was adjusted to allow for several (3 or 4) annihilation events for 
each bubble chamber expansion. Three stereoscopic pictures were taken of each expansion to enable 
a three-dimensional reconstruction of the tracks. 

Scanning the films and reconstructing events was a major enterprise. The spatial coordinates 
of four points for each track were measured from the films, with a precision of 80 /xm. From the 
coordinates the charged-particle momenta were determined. We estimate the momentum resolution 
for 928 MeV/c pions from the reaction pp n^n~ to be 25 MeV/c. 

A total of 1.6 X 10^ events were recorded at CERN, 7.5 x 10^ at BNL. These numbers exceeded 
the scanning capabilities available at that time, and only a fraction of the data was analysed: about 
80,000 events at CERN and 45,000 events at BNL. From the momentum of the incoming antipro- 
tons their range was estimated and compared to the true range; thus contamination due to in flight 
annihilation could be avoided, or at least reduced. 

The Brookhaven results are documented in Refs. [42-50]. At CERN, more aspects of the anni- 
hilation process were investigated, leading to a larger number of pubUcations. An incomplete list of 
publications in refereed journals includes [2,51-76] for H2-filled and [77-84] D2-fiUed chambers. 
Further publications discussed the interpretation of these results. 

In the 1970's, a first set of counter experiments were performed at BNL to study 7-rays from 
antiproton annihilation. Antiprotons from a separated beam were stopped in a liquid H2 or D2 
target. Photons were detected by their conversion in Cu(Pb) plates sandwiched between scintillation 
counters. In some experiments a Nal detector, surrounded by scintiUation counters, was used to 
measure 7-rays with better resolution. Data from the Rome-Syracuse collaboration, taken with the 
D2-filled BNL bubble chamber, were analysed in parallel. Results can be found in [85-101]. 

2.2 Experiments at KEK 

This experiment was designed to search for narrow lines in the momentum distributions of 7r° and r] 
from pp annihilation. The initial aim was to find narrow multiquark or quasi-nuclear bound states 
[102-105]. Later, frequencies for annihilation into two narrow mesons were determined with both 
H2 [106, 107] and D2 [108] targets. 

The layout of the experimental set-up is shown in Fig. 1 . A full description of the detector can 
be found in [ 1 02 , 1 03 1 . Antiprotons at 5 80 MeV/ c, produced at the KEK 1 2 Ge V proton synchrotron, 
were degraded in a graphite slab and stopped in a liquid H2 target of 14 cm diameter and 23 cm in 
length. The p beam used double-stage mass separation and a contamination ratio of e/xTr/p ~ 8 was 
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Figure 1 : Side and end view of tlie KEK detector. 



obtained with a typical stopping intensity of 270 p/synchrotron pulse. The charged particles pro- 
duced in pp annihilation were detected with scintillation counter hodoscopes and tracked with cylin- 
drical and planar multiwire proportional chambers whose total coverage was 93% x An sr. Photons 
were measured with a calorimeter consisting of 96 Nal(TI) crystals surrounded by 48 scintillating 
glass modules and assembled into a half barrel [109J, see Fig.l. The geometrical acceptance for tt" 
increased from 10.5% at a tt" energy of 500 MeV to 14.5% at 900 Me V. The overall energy reso- 
lution at FWHM for photons was approximately AEj/Ej = 6.2%/(i?-y in GeV)^/"^, for energies 
above 80 MeV. 

Events were recorded when the detector signalled that a slow antiproton was incident on the Hq- 
uid H2 target and one or two photons were measured in the Nal detector. A fast cluster counting logic 
counted the multipUcities of charged and neutral clusters separately, and if they satisfied preselected 
criteria the data was recorded. In some of the later experiments [107] an additional small (1.3% x Air 
sr.) BGO detector surrounded by Nal modules was used. No cluster counting logic was used in this 
case and the energy resolution (FWHM) was estimated to be AE^/E^ = (S.S%/{E^ in GeV)^/''. 
As the Nal photon spectrometer has less than 27r acceptance, for measurements of two-body branch- 
ing ratios it is not possible to detect both mesons. In this case the existence of the second meson and 
its mass are deduced from the inclusive energy spectrum recorded for a single 7r° or rj. 

2.3 Cooled antiproton beams 

Early experiments with electronic detection techniques, including those carried out at KEK, used 
partially separated secondary antiproton beams produced from an external target. These beams 
were characterised by a relatively low rate of stopped antiprotons over a large volume and with a 
contamination of unwanted particles efiir/p ~ 100 in the earliest experiments to e/^7r/p ~ 8 in the 
more recent ones. This situation was transformed by the availability of cooled antiproton beams at 
CERN, together with the construction and operation of the LEAR facility. 

A description of the cooled antiproton beams used at CERN has been given in a previous review 
article [1]. For a proton beam of 23 GeV/c incident on a Be target, antiprotons are produced with 
a broad maximum in momentum at 3.5 GeV/c. The use of cooUng allows these antiprotons to be 
decelerated to low momenta whilst keeping the same flux. Additionally, coohng gives the antiproton 
beams a small size and a reduced momentum. 

The LEAR facility was constructed at CERN to handle pure antiproton beams in the momentum 
range from 105 MeV/c to 2000MeV/c with small physical size and a typical momentum spread 
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of Ap/p ^ 0.1%. This small momentum spread for low momentum protons gave a very small 
stopping region in liquid H2 and D2 targets and also enabled gas targets to be used. The use of gas 
targets is particularly important since the fraction of P-state annihilation is considerably increased 
in gaseous H2 targets due to the reduced effect of Stark mixing. An ultra-slow extraction system 
enabled essentially DC beams to be produced, with spills typically lasting several hours. Typical 
beam intensities were in the range 10* to 10^ p/sec. The beam purity was 100%. 

The LEAR project was approved by CERN in 1980, and in July 1983 the first antiproton beams 
were deUvered to users. After a break in 1987, to construct a new Antiproton Collector (ACOL) 
which resulted in a flux gain of a factor of 10, the facility was operated until the end of 1996, when it 
was closed for financial reasons. The Asterix, Obelix and Crystal Barrel (CBAR) experiments were 
all carried out at LEAR; Asterix in the first, Obelix and Crystal Barrel in the second phase. 

2.4 Detectors at LEAR 

2.4.1 PS 171: The Asterix experiment 

Liquid targets were used in both the bubble chamber and counter experiments described earlier 
(Sec. 2.1 and 2.2). In liquid H2 or D2, annihilation occurs at rest and is preceded by capture of 
an antiproton by a hydrogen or deuterium atom. Collisions between the protonium atom and H2 
molecules induce transitions from high orbital angular momentum states via Stark mixing; and this 
mixing is fast enough to ensure dominant capture from S-wave orbitals. In H2 gas, the collision 
frequency is reduced and P-wave annihilation makes significantly larger contributions. In particular 
at very low target pressures the P-wave fractional contribution is very large. Alternatively, rather 
pure samples of P-wave annihilation can also be studied by coincident detection of X-rays emitted 
in the atomic cascade of the pp system (which feed mostly the 2P level). 

The Asterix experiment was designed to study pp annihilation from P-wave orbitals by stop- 
ping antiprotons in H2 gas at room temperature and pressure and observing the coincident X-ray 
spectrum. The detector, shown in Fig. 2, consisted of the following main components: 

1. A gas target of 45 cm length and 14 cm in diameter contained the full p stop distribution for 
an antiproton beams at 105 MeV/c. 

2. The target was surrounded by a X-ray drift chamber (also used to improve the tracking ca- 
pabiUty and for particle identification via dE/dx). The energy resolution of the detector for 8 
keV X-rays was about 20%. Pions and kaons could be separated up to 400 MeV/c. The target 
and X-ray drift chamber were separated by a 6 /um aluminised mylar foil to guarantee gas 
tightness and good X-ray transmission even at low energies. 

3. Charged particles were tracked in a set of seven multi-wire proportional chambers, partly with 
cathode readout to provide spatial resolution along the wires. The momentum resolution for 
pp — !■ 7r+7r~ events at 928 McV/c was 3%. 

4. A one-radiation-length lead foil in front of the outer chambers permitted reconstruction of the 
impact points of photons. 

5. Two end-cap detectors with three wire planes and cathode readout on both sides gave large 
solid-angle coverage. A lead foil was mounted behind the first chamber. The end cap detectors 
were used to identify 7's but not for reconstruction of charged tracks. 

6. The assembly was situated in a homogeneous magnetic field of 0.8 Tesla. 

With the experimental resolution of the detector, there was nearly no background for fully- 
constrained final states and up to 14% for final states with one missing tt^. 

The main data sets taken with the Asterix detector consisted of 1.38 x 10^ events with two long 
tracks (passing at least the first five chambers) without triggering on X-rays, 2.13 x 10^ events with 
two long tracks with a trigger on X-rays, and 1.89 x 10^ events with four long tracks and with the X- 
ray trigger. The "long-track" criterion guaranteed that the particles reached the outermost chambers 
and gave optimum momentum resolution. The X-ray enhancing trigger had an efficiency of 25%; 
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Figure 2: Side and front view of tine PS1 72 Asterix detector. 



one quarter of the triggered events had — after all cuts — an identified low-energy X-ray. There was 

a contamination from Bremsstrahlung X-rays of about 15% in the X-ray data sample. 
The detector is fully described in [1 10]. Physics results were published in [11 1-126]. 



2.4.2 PS 201: The Obelix experiment 

The layout of the ObeUx spectrometer is shown in Fig. 3. A full description of the detector can be 
found in [127, 128]. It consists of four sub-detectors arranged inside and around the open-axial field 
magnet which had previously been used for experiments at the ISR. The magnet provides a field of 
0.5 T in an open volume of about 3 m^. The subdetectors are: 

1 . A spiral projection chamber (SPC): an imaging vertex detector with three-dimensional readout 
for charged tracks and X-ray detection. This detector allowed data to be taken with a large 
fraction of P-wave annihilation, and to measure angular correlations between X-rays from the 
pp atomic cascade and annihilation products. 

2. A time-of-flight (TOP) system: two coaxial barrels of plastic scintillators consisting of 30 
(84) slabs positioned at a distance of 18 cm (136 cm) from the beam axis; a time resolution of 
800 ps FWHM is achieved. 

3. A jet drift chamber (JDC) for tracking and particle identification by dE/dx measurement 
with 3280 wires and flash-analog-to-digital readout. The chamber was split into two half- 
cylinders (160cm in diameter, 140cm long). The intrinsic spatial resolution was cr^ = 12 mm, 
ar<p = 200 /xm; the momentum resolution for monoenergetic pions (with 928 MeV/c) from the 
reaction pp — > tt+tt" was found to be 3.5%. 

4. A high-angular-resolution gamma detector (HARGD) [ 1 27] . The calorimeter consisted of four 
modules made of layers of 3 x 4 lead converter foils with planes of limited streamer tubes 
as the active elements. Twenty converter layers, each 3 mm thick, were used corresponding 
to a total depth of about 10 radiation lengths. Due to their excellent spatial resolution, good 
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Figure 3: Schematic view of tlie Obelix experiment set-up. Tlie numbers indicate the main compo- 
nents of the apparatus: the Open Axial Field magnet (1), the SPG (2, 4), the TOF (3), the JDC (5), 
the HARGD (6). 

The detector system allowed a variety of targets to be used: a liquid H2 target, a gaseous H2 
target at room temperature and pressure, also a target at low pressures (down to 30 mbar). The wide 
range of target densities could be used to study in detail the influence of the atomic cascade on the 
annihilation process. The H2 could also be replaced by D2. A further special feature of the detector 
was the possibility to study antineutron interactions. The n beam was produced by charge exchange 
in a liquid H2 target (positioned 2 m upstream of centre of the main detector). The intensity of the 
coUimated beam was about 40 n/lO^p of which about 30% interact in the central target. The ri beam 
intensity was monitored by a downstream n detector. 

The Obelix Collaboration had a broad program of experiments covering atomic, nuclear and 
particle physics [128]. The main results can be found in [1 1, 128-162]. 

2.4.3 PS 197: The Crystal Barrel experiment 

The main objective of the Crystal Barrel experiment was the study of meson spectroscopy and in 
particular the search for glueballs (gg) and hybrid (gqq) mesons produced in pp and pd annihilation 
at rest and in flight. Other objectives were the study of pp and pd annihilation dynamics and the 
study of radiative and rare meson decays. A particular feature of the experiment was its photon 
detection over a large soUd angle with good energy resolution. Physics results are published in 
[4,41,163-217] 

The layout of the Crystal Barrel spectrometer is shown in Fig. 4. A detailed description of the 
apparatus, as used for early data-taking (1989 onwards), is given in [218]. To study annihilation 
at rest, a beam of 200 MeV/c antiprotons, extracted from LEAR, was stopped in a 4 cm long Uquid 
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hydrogen target at the centre of the detector The whole detector was situated in a 1 .5 T solenoidal 
magnet with the incident antiproton beam direction along its axis. The target was surrounded by 
a pair of multiwire proportional chambers (PWC's) and a cylindrical jet drift chamber (JDC). The 
JDC had 30 sectors with each sector having 23 sense wires at radial distances between 63 mm and 
239 mm. The position resolution in the plane transverse to the beam axis was a — 125 fim. The 
coordinate along the wire was determined by charge division with a resolution of = 8 mm. This 
gave a momentum resolution for pions of cr/p ~ 2% at 200 MeV/c, rising to 7% at 1 GeV/c for those 
tracks that tracked all layers of the JDC. The JDC also provided tt/K separation below 500MeV/c 
by ionisation sampling. 

The JDC was surrounded by a barrel shaped calorimeter consisting of 1380 CsI(Tl) crystals in 
a pointing geometry. The Csl calorimeter covered the polar angles between 12° and 168° with full 
coverage in azimuth. The overall acceptance for shower detection was 0.95 x 4tt sr Typical photon 
energy resolutions for energy i? (in GeV) were crs/i? = 2. and cr^ e — 1.2° in both polar 
and azimuthal angles. The mass resolution was a = 10 MeV for tt" and 17 MeV for rj 2j. 




1m 



Figure 4: Overall layout of the Crystal Barrel detector showing (1) magnet yoke, (2) magnet coils, (3) 
Csl barrel, (4) jet drift chamber, (5) proportional chamber, (6) liquid hydrogen target, (7) one half of 
endplate. Left - longitudinal cross section; Right- transverse view. 

In 1995 the PWC's were replaced by a microstrip vertex detector (SVTX) consisting of 15 single- 
sided silicon detectors, each having 128 strips with a pitch of 50 /im running parallel to the beam 
axis [219,220]. (See [220, Fig. 1] for an overall view of the detector) As well as giving improved 
identification of secondary vertices, this detector provided better vertex resolution in r, (j) and im- 
proved momentum determination with a resolution Ap/p for charged tracks of 3.4% at 0.8 GcV/c 
and 4.2% at l.OGeV/c. 

To study annihilation in hydrogen gas, the liquid target was replaced by a 12 cm long Mylar 
vessel with 230/xm thick walls and a 195/xm thick entrance window, containing hydrogen gas at 
room temperature and 12 bar pressure. A 55/im thick Si detector was used to count the incident 
105 MeV/c antiproton beam. 

A particular feature of the detector system was a multi-level trigger [218] on charged and neutral 
multiplicities and on invariant mass combinations of the neutral secondary particles. This allowed 
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the suppression of well-known channels and the enhancement of rare channels of specific interest. 
The PWC/SVTX and the inner layers of the jet drift chamber determined the charged multiplicity of 
the final state. Events with long tracks could be selected to give optimum momentum resolution by 
counting the charged multiplicity in the outer layers of the JDC. A hardwired processor determined 
the cluster multiplicity in the Csl barrel, whilst a software trigger, which was an integral part of the 
calorimeter read out system, allowed a trigger on the total deposited energy in the barrel or on the 
TT^ or T] multiplicity. 

Typical incident beam intensities were lO^p/sec at 200MeV/c for stopping in liquid H2 or 
105 MeV/c when a 12 bar gas target was used. For experiments to study interactions in flight, larger 
intensities in the range from 10^ to 10® p/sec were required at beam momenta in the range 600 to 
1940MeV/c. 

A convenient summary of the data taken by the experiment, both at rest and in flight, for liquid 
H2 liquid D2 and gaseous H2 targets has been given by Amsler (see [4, Table 1]). Typical data sets 
contain 10^ to 2 x 10'' events. 

2.5 Future experiments 

With the closure of the LEAR facility in 1996, an era of intensive experimental study of low and 
medium energy antiproton annihilations came to an end. CERN has continued its involvement in the 
production of p beams with the construction of the AD (Antiproton Decelerator) [22 1 ] . This provides 
p beams with momentum from 100 to 300 MeV/c but without slow extraction. Slow extraction is 
not possible without making major modifications to the AD [222,223]. The space for experiments 
[224] is very limited and the experimental program is solely devoted to the production of trapped 
antihydrogen and studies of the formation and cascade in light antiprotonic atoms. 

For many years Fermilab has had the world's most intense antiproton source. However the 
opportunities for medium energy p physics have been very limited and experiments have focused 
on the observation and measurement of charmonium states. There has been no low-energy program. 
The possibility of building a new antiproton facility which could decelerate p to below 2 GeV/c for 
injection into a new storage ring has been discussed [225]. The storage ring would be equipped with 
RF to decelerate p down to the hundreds of MeV/c range. There are however, as yet, no firm plans 
to build such a facility. 

The design of the Japan Proton Accelerator Research Complex (J-PARC, [226]) is well under 
way. The J-PARC project includes studies of particle and nuclear physics, materials science, life 
sciences and nuclear technology. The accelerator complex consists of an injection linac, a 3 GeV 
synchrotron and a 50 GeV synchrotron [227]. At this latter machine nuclear and particle physics 
experiments using neutrinos, antiprotons, kaons, hyperons and the primary proton beam are planned. 
At one stage it was hoped that the LEAR facility could be moved to J-PARC. However the LEAR 
ring is now required for the injection of heavy ions into the CERN LHC. It also seems likely that 
neutrino physics and the study of rare kaon decays wiU be topics for the first experiments at JFK and 
the construction of a dedicated low/medium energy antiproton facility is now some way off. 

At the GSI laboratory in Darmstadt the construction of a new facility is planned and conditionally 
approved [228], the International Facility for Antiproton and Ion Research, FAIR. A conceptual 
design report [229] outlines a wide physics program and the envisaged accelerator complex. In 
particular a High Energy Storage Ring (HESR, [230]) will provide stored antiproton beams [231- 
233] in the range 3 to 15 GcV/c with very good momentum resolution (Ap/p « 10~^). The charm 
region is thus accessible with high rates and excellent resolution. A target inside the storage ring 
will be used, together with a large multi-purpose detector for neutral and charged particles with 
good particle identification [234] . Production and use of polarised antiprotons is a further future 
option for studying spin aspects of antiproton-proton scattering and annihilation. It is planned to 
broaden the program by including a low-energy component, FLAIR for the study of antimatter and 
highly-charged ions at low energies or nearly at rest. 
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3 Mesons and their quantum numbers 

3.1 qq mesons and beyond 

Since the annihilation process leads to production of mesons, it is useful to recall some basic defini- 
tions and properties of the meson spectrum. 

Mesons are strongly-interacting particles with integer spin. The well established mesons have 
flavour structure and other quantum numbers which allows us to describe them as bound states of a 
quark and an antiquark. These valence quarks which describe the flavour content are surrounded by 
many gluons and quark-antiquark pairs. Other forms of mesons are also predicted to exist: glueballs 
should have no valence quarks at all; in hybrids, the hypothetical gluon string transmitting the colour 
forces between quark and antiquark is supposed to be dynamically excited; and multiquark states 
are predicted, described either as states of (qqqq) or higher valence-quark structure, or as meson- 
meson or baryon-antibaryon bound states or resonances. These unconventional states are presently 
searched for intensively; they are however not the subject of this review. 

Quarks have spin s = 1/2 and baryon number B = 1/3, antiquarks s = 1/2 and B = —1/3. 
Quark and antiquark combine to B = and to a spin triplet (S = 1) or singlet (S = 0). In 
conventional mesons, the total spin S of the quark q and the antiquark q, and the orbital angular 
momentum L between q and q couple to the total angular momentum J of the meson: J = L + S. 
Light mesons are restricted to u, d, and s quarks. 



3.2 Quantum numbers 

Parity: The parity P of a meson involves the orbital angular momentum L between quark and 
antiquark and the product of the intrinsic parities which is PqPq = —1 for a fermion and its antipar- 
ticle: 

P= (-1)^+1. (3.1) 



Charge conjugation: Neutral mesons are eigenstates of the charge conjugation operator 

C = {-\f+^. (3.2) 

It turns out convenient to use the same sign convention within a multiplet. For instance, since 
C7r° = 7r°, we choose C-k^ = tt'f, and CK° = K°, while since Cp° = we adopt Cp^ = -p^ 
andCK*± = -K*=F. 



Isospin: Proton and neutron form an isospin doublet and so do the up and the down quark. We 
define antiquarks by u = Cu and d = C d, antinucleons by p = G n and n = G p. This means 
we use the 2 representation of SU(2) for {u, d} and the 2 representation for {p, n}. See, e.g., [235] 
for a detailed discussion on phase conventions for isospin states of antiparticles. We obtain: 



\I 


- hh - 


1) = 


-\du) , 


\I 




0) = 


- ^{\uu)-\dd)) 


I = 


= l,/3 = 


-1) = 


\ud) , 


1^ 


= 0,/3 = 


0) = 


- j,{\uu) + \dd)) 


1^ 


= 0,/3 = 


0) = 


- \ss) . 



The I nn ) and | ss ) states have the same quantum numbers and mix to form two physical states. 
With n we denote the two Ughtest quarks, u and d, while n stands for the neutron. 



Mesons and their quantum numbers 



15 



The G-parity: The G-parity is defined as charge conjugation followed by a rotation in isospin 
space about the y-axis, 

G = Ce'^^'y = {-lyC = (-1)^+^+^ , (3.4) 

and is approximately conserved in strong interactions. It is a useful concept since G = (—1)"" for a 
system of n-^ pions. This generalises the selection rule for c+c^ nj in QED, namely C = (— 1)". 

The a2(1320), for instance, decays into pn with p tttt, hence into three pions, and into rin 
with one pion in the final state. The a2(1320) never decays into tttt or ryTTTr: G-parity is conserved. 
The rj having G = +1 nevertheless decays into three pions; the ui has a small partial width for 
decays into two pions. These decay modes break isospin invariance; they vanish in the limit where 
u- and d-quark have equal masses and electromagnetic interactions are neglected. KK pairs may 
have G = — 1 or +1. 



3.3 Meson nonets 

Mesons are characterised by their quantum numbers J^'-'' and their flavour content. Quark-antiquark 
states with quantum numbers J^^ are often referred to by the spectroscopic notation n'^^^^Lj 
borrowed from atomic physics. In the light-quark domain, any n '^^^^Lj leads to a nonet of states. 
Based on SU(3) symmetry, we expect an octet and a singlet. However, the s quark is heavier than 
the u and d quark. This results into SU(3) breaking. The actual mesons can be decomposed either 
on a basis of SU(3) eigenstates or according to their uu, dd and ss content. 



3.3.1 The pseudoscalar mesons 

The pseudoscalar mesons correspond to J^*-^ = and n^^'~^^Lj= I^Sq. The nine orthogonal 
SU(3) eigenstates are shown in Fig. 5. The quark representation of the neutral members is 



TT 



,0 



1 . T,. /I. T, „ s /I 



{uu — dd), r]s = \ -{uu + dd — 2ss) , rji = \ -{uu + dd + ss) . (3.5) 



The actual mesons rj and tj' can be written as 



S Octet 




Figure 5: The nonet of pseudoscalar mesons. 
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It? ) =coseps|??8) -sineps|v7i) 
I ry' ) = sin Sps I r/8 ) + cos Ops 1 771 ) 

with the pseudoscalar mixing angle Ops- A mixing angle 0pg = = arctan(l/\/2), called 
the ideal mixing angle, would lead to a decoupling 77 oc \ss) and rj' oc \nn). The rj and rj' wave 
functions can be decomposed into the fin and ss basis. With O = 6pg — Ojd + 7r/2, 

\r] ) = cos Q\hn) — sin 8 | ss ) = | ) + | ss ) 
1 77' )= sin 6 I nn ) + cos | ss ) = | fin ) + Y^/ | ss ) 

The mixing angle can be determined experimentally from -q and r/' production and decay rates. It 
is shown in [236] that the quark flavour basis nn and ss is better suited to describe the data than the 
octet-singlet basis. The latter can describe data only by the introduction of a second mixing angle. 

The r] and 77' could also mix with other states, in particular radial excitations or glueballs. The 
77' is nearly a flavour singlet state and can hence couple directly to the gluon field; this has led to 
speculations that the 77' (and to a lesser extend also the 77) may contain a large fraction of glue. This 
requires an extension of the mixing scheme (3.7) by introduction of a non-qq or inert component, 
with a third state of unknown mass which could, e.g., be dominantly a gluebaU. 



\r])= X^\nn) + \ss) + ^Jglue) 

1 77')= X^>\nn) + Yr,'\ss) + Z^'|glue) (3.8) 
light quark strange quark inert 

At present there is no convincing evidence for a glueball content in the 77' wave function, and we 
will assume and Z^' to vanish, i.e., = Z^i ~ 0. 

3.3.2 Other meson nonets 

A meson nonet consists of five isospin multiplets. The pseudoscalar nonet, for instance, contains the 
pion triplet, two kaon doublets, the 77' and the 77. Well known are also the nonet of vector mesons 
with quantum numbers J^"^ = 1 (three p, four K*, and w), and the nonet of tensor mesons 
with quantum number 2++ (a2(1320), K^(1430), f2(1525), f2(1270)). In a spectroscopic notation, 
these are the l^Si and 1^P2 states. Both nonets have a nearly ideal mixing angle Oid = 35.3° for 
which one meson is a purely nn and the other one a purely ss state. These are the u) and (/)(1020), 
and the f2(1270) and f2(1525) mesons, respectively. Note that the mass difference between the ss 
and the nn state is about 250 MeV. This mass difference is due to the larger constituent mass of 
strange quarks. 

The mixing angles for these meson nonets (all except the pseudoscalar nonet) are defined as 

I a; ) = cos 6y I nn ) — sin 9y I ss ) 

_ _ (3.9) 

\^) = sin Qv I nn ) + cos Qy \ ss ) 

In Table 1 some meson nonets are collected. The assignment shown here is reproduced from 
the quark-model description of Amsler and Wohl in [237] and represents one possible scenario. In 
particular the scalar-meson nonet is hotly debated [238,239] but there are also open questions in the 
axial-vector nonet [240, 241] and for the radial excitations of vector [242] and pseudoscalar [243] 
mesons. 

3.3.3 The Gell-Mann-Okubo mass formula 

The Gell-Mann-Okubo mass formula relates the masses of a meson nonet and its mixing angle. It 
can be derived by ascribing to mesons a common mass Mq plus the (constituent) masses of the quark 
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Table 1 : The light mesons. The two mesons Kia and Kib mix to form the observed resonances 
Ki(1280) and Ki(1400). The scalar mesons resist an unambiguous classification; the scenario repro- 
duced here assumes that the an (980) and f()(980) are KK molecules or generated dynamically but 
are not qq states, that the fo(1370) is a qq state and not generated dynamically and that the fo(1500) 
is a glueball. The fi(1510) is discarded in the listing below. There are considerable difficulties for the 
nonets of pseudoscalar and vector radial excitations. 
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J 
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1=1 


7=1/2 


7 = 


7 = 


jPC 















TT 


K 


V 




0-+ 


I'So 





1 


1 




P 


K* 


<t> 




1— 


l^Si 


1 





1 




bi(1235) 




hi(1380) 


hi(1170) 


1+- 


l^Pi 


1 


1 







ao(1450) 


K5(1430) 


fo(1710) 


fo(1370) 


0++ 


l3Po 


1 


1 


1 




ai(1260) 




fi(1420) 


fi(1285) 


1++ 


l3Pi 


1 


1 


2 




32(1320) 


K*(1430) 


f2(1525) 


f2(1270) 


2++ 


l3P2 


2 





2 




772(1670) 


K2(1770) 


??2(1870) 


r/2(1645) 


2-+ 


IID2 


2 


1 


1 




p(1700) 


K*(1680) 


^!)(????) 


a;(1650) 


1~ 


l^Di 


2 


1 


2 




P2(????) 


K2(1820) 


.^2 (????) 


a;2(????) 


2~ 


13D2 


2 


1 


3 




/03(1690) 


K^(1780) 


<^;i(1850) 


W3(1670) 


3~ 


l^Dg 











2 


7r(1370) 


K(1460) 


77(1440) 


77(1295) 


0-+ 


2^80 





1 


1 


2 


p(1450) 


K*(1410) 


(/>(1680) 


w(1420) 


1~ 


23Si 



and antiquark it is composed of. The relation is written as 

20 3M^ + - AMk 
tan B = —-r • (3.10) 

AMk - ^M,y - ' 

Often, the linear GMO mass formula is replaced by the quadratic GMO formula which is given as 
above but with values instead of masses M. Note that in the limit of chiral symmetry quark 
masses are proportional to the mass square of the meson masses. The quadratic GMO formula reads 

„ ^ 3M2 + Ml - 4A7| 

tan^0 = -nr75 ^ttt, ■ (3-11) 

4M|. - 3M2, -Ml 

Table 2 gives the mixing angles derived from the Unear and quadratic GMO formula. 

Table 2: Mixing angles of meson nonets 



Nonet members 


€^ linear 


®quad 


7r,K,77',77 


-23° 


-10° 




36° 


39° 


a2(1320), K^(1430), f2(1525), f2(1270) 


26° 


29° 


P3(1690), K^(1780), <^3(1850), a;3(1670) 


29° 


28° 



3.3.4 The Zweig rule 

The "quark Une rule", or Zweig rule, is also called the "OZI rule", after Okubo, Zweig and lizuka, 
or even the "A-Z rule" to account for all the various contributions to its study. See, e.g., Ref. [244], 
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for a comprehensive list of references. This rule has played a crucial role in the development of the 
quark model. 

For instance, the 0(1020) is a vector meson with isospin / = 0, seemingly similar to aj(780), 
but much narrower, in spite of the more favourable phase-space. It decays preferentially into KK 
pairs, and rarely into three pions. The explanation is that the 0(1020) has an almost pure ss content, 
and that the decay proceeds mostly with the strange quarks and antiquarks flowing from the initial 
state to one of the final mesons, as per Fig. 6, left, while the process with an internal ss annihilation 
(centre) is suppressed. The decay 4> tttttt is attributed to the main ss component being slightly 
mixed with a qq component, which in turn decays into pions by a perfectly allowed process (right). 

The rule is quite strictly observed, as the non-strange width of the 0(1020) is less than 1 MeV! 
Even more remarkably, the rule works better and better for the charm and beauty analogues, for 
which the decay into naked-flavour mesons is energetically forbidden: the total width is only about 
90 keV for the J/*(cc) and 50 keV for T{bb). 



s 



s 




Figure 6: Connected (left) and disconnected (centre) contribution to 0(1020) decay. The latter con- 
tribution to 0(1020) ^ TTTTTT can be described as an allowed decay from a small impurity in the wave 
function (right). 

The magnitude of OZI violation in mesons is primarily described by the mixing angle, though an 
OZI violation from the decay amplitude cannot be excluded. The vector mesons, e.g., have a mixing 
angle of Qy — 39° which deviates from the ideal mixing angle 8id hy S ~ 3.7°. The physical 
0(1020) is then written 

I 0) = cos(5 1 ss) + sin(5 1 nn6). (3.12) 

If there are no s-quarks in the initial state we expect the production of mesons to be suppressed 
compared to uj production by sin^ S ~ 0.042. Indeed, a (jy/uj ratio of 0.0032 ± 0.0004 was found 
at Argonne in the reaction Tr^p 0(w)p at 6 GeV/c [245]. This has to be compared with the 
extremely large values (up to 0.6 !) found in WA56 data on 7r~p — > (0(w) + 7r)p at 12 and 20 
GeV/ c [246]. In pp annihilation at rest, the 0/a; production ratios were found to depend strongly on 
the pp initial state and on the recoiling particles. These results will be discussed further in Sec. 8. 

There is a wide consensus that ideal mixing and OZI rule can be, if not derived, at least justified 
from QCD, but there are different approaches: the 1/Nc expansion, where Nc is the number of 
colour degrees of freedom, cancellations of loop diagrams, lattice simulations, instantons effects, 
etc. The generally accepted conclusion is that ideal mixing is nearly achieved for mesons, except in 
the scalar and in the pseudoscalar sectors. See, e.g., Ref. [247] and references there. 

3.3.5 Meson decays 

The decays of mesons belonging to a given nonet are related by SU(3) symmetry. The coefficients 
governing these relations are called SU(3) isoscalar factors and listed by the Particle Data Group 
[237]. We show here two simple examples. 

A glueball is, by definition, a flavour singlet. It may decay into two octet mesons, schematically 
1 ^ 8 X 8. The isoscalar factors for this decay are' 

(glueball) ^ (KK , tttt , i^sVs , KK) - -^=(2 ,3,-1, -2)^/2 . (3.13) 

'The parenthesis reads (>/2, \/3, —1, — \/2), and similai'ly for Eqs. (3.14). 
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Hence glueballs have squared couplings to KK, tttt , rj^ris of 4 : 3 : 1 . The decay into two isosinglet 
mesons Tjirii has an independent couphng and is not restricted by these SU(3) relations. It could be 
large leading to the notion that the r]' couples strongly to glueballs, and that r]' are gluish. The decay 
into rjirjs is forbidden: a singlet cannot decay into a singlet and an octet meson. This selection rule 
holds for any pseudoscalar mixing angle: the two mesons tj and rj' have orthogonal SU(3) flavour 
states and a flavour singlet cannot dissociate into two states which are orthogonal. 

As a second example, we choose decays of vector mesons into two pseudoscalar mesons. We 
compare the two decays K* Ktt and p tttt. These are decays of octet particles into two octet 
particles. This 8 — > 8 x 8 coupling is either symmetric (8i — > 8 x 8) or antisynnmetric (82 — > 8 x 8) 
under the exchange of the final-state particles. The two pseudoscalars in K* or p decay having orbital 
momentum i = I, one should use the antisymmetric flavour couphng, 82 ^ 8x8, whose isoscalar 
factors are 

(K*) ^ (Ktt Krj ttK r]K) = -^{3 ,3,3, -3)^2 

(3.14) 

(p) (KK TTTT r]TT TTT] KK) = ^=(2 ,8,0,0, -2)^/2 



0.40 (3.15) 



/12 

Hence we derive K* Ktt + ttK oc 6, p ^ tttt oc 8, or 

rK.^K,.+,.K 6 /0.291 



Tp^^^ 8 V 0.358 



The latter factor is the ratio of the decay momenta q to the 3rd power The transition probability is 
proportional to q; for low momenta (or point-hke particles), the centrifugal barrier scales with q'^^ 
where £ is the orbital angular momentum. 

From data we know that the width ratio is 0.34 and so the relations are fulfilled at the level 
of about 20%, a typical magnitude for SU(3) breaking effects. We have neglected many aspects: 
the transition rates are proportional to the squared matrix element (given by SU(3)) and the wave 
function overlap. The latter can be different for the two decays. Mesons are not point-like; the 
angular barrier factor should hence include Blatt-Weisskopf corrections [248]. An apphcation of 
SU(3) to vector and tensor mesons can be found in [249]. 
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4 Kinematics and conservation laws 

In this section, we discuss the kinematics of annihilation into two, three, or more mesons, and the 
selection rules due to exactly or approximately conserved quantum numbers. 

4.1 Kinematics 

We consider annihilation at rest. The formalism below can be used for annihilation in flight in the 
centre of mass, if 2m is replaced by s^/^, where m is the proton or antiproton mass, and s the usual 
Mandelstam variable s = (pi + p2)^ built of the four-momenta of the initial proton and antiproton. 
If mi, 1712, etc. denote the mass of the final mesons, annihilation is possible if 

2m > mi + m2 + • • • , (4.1) 

Up to 13 pions could be produced. 
4.1.1 Two-body annihilation 



Pi = [Ei,Pi] 



P2 = [£'2,P2] 



P = [2m, 0] 



Figure 7: Kinematics for two-body annihilation 

The notation is defined in Fig. 7. From the energy-momentum balance rewritten as p2 ~ P ~Pi 
and squared, one gets 



4m^ — m2 + rrii 
El — — 



\Pl\ = \P2\ 



4m 

[4m^ — (mi + m2)^] ^ [4m^ — (mi — m2)^] 



(4.2) 



4m 

In case of two identical mesons the momentum can be written in the form 



\pi\ = y m'^ - ml . (4.3) 

The occurrence of two narrow mesons is easily observed due to the narrow momentum distribution of 
the two produced particles. Other annihilation modes such as pp pn or a2(1320)7r involve short- 
lived resonances. In these cases, the measurement of the annihilation frequency is more complicated. 

Due to its relatively large mass, the antiproton-proton system possesses a large variety of two- 
body annihilation modes. In Table 3 we list the momenta of typical reactions. For broad resonances 
the momenta are calculated for the nominal meson masses. Some annihilation modes like pp 
wf2(1270) are at the edge of the phase space and only the part of the £2 (1270) below 1094 MeV is 
produced. They may, nevertheless, make a significant contribution to the annihilation process. 

The mean decay length, Di — j/3ct — hp/ (mF), varies over a wide range. For pp iT^rj, the 



-^0 



TT . 



mean path of the is Z)/ = 1 nm within the 77 Ufetime. Assume the 77 then decays into tt+tt 
Because of the large Di there will be no interaction between the tt*^ recoiling against the rj and 
the pions from 77-decays. In case of irp annihilations, the recoil pion travels 7.5 fm within the p 
mean live time, and rescattering of the primarily produced pion and pions from p decays is unlikely. 
The situation is different for production of two short-lived high-mass mesons. In annihilation into 
K*K* + c.c, the K*'s have moved only a mean distance of 1.3 fm when they decay and interactions 
between the kaons and pions from K* decays are likely to occur. 
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Table 3: Momenta for pp annihilation at rest into two mesons. 



Channel 



PP 



lel 


Momentum 


Channel 


Momentum 




928.5 MeV/c 


PP ^ PP 


536.3 MeV/c 




927.8 MeV/c 


LOp 


527.5 MeV/c 




852.3 MeV/c 


UJLO 


518.5 MeV/c 


K+K- 


797.9 MeV/c 




499.7 MeV/c 




795.4 MeV/c 


7rf2(1270) 


491.2 MeV/c 




773.2 MeV/c 


7ra2(1320) 


459.9 MeV/c 


TT°LU 


768.4 MeV/c 


v'p 


364.4 MeV/c 


rji] 


761.0 MeV/c 


rj'uj 


350.5 MeV/c 


VP 


663.5 MeV/c 


K*K* 


285.2 MeV/c 




658.7 MeV/c 


P0 


280.3 MeV/c 


rjuj 


656.4 MeV/c 


uxp 


260.8 MeV/c 




652.4 MeV/c 




206.2 MeV/c 


KK* 


616.2 MeV/c 


?7a2(1320) 


91.2 MeV/c 


vv' 


546.1 MeV/c 


t^f2(1270) 


See text 



4.1.2 Three-body annihilation 

Unlike the two-body case, the energy of a given particle can vary over a certain range. Even for 
equal masses, the symmetric star of Fig. 8 (left) is only a very particular case. The minimal energy 




2m , 



Figure 8: Notation for the kinematics of three-body annihilation (right). Even for equal masses, the 
symmetric star (left) is just one possibility among many others. 

of particle 1, for instance, is obviously obtained when it is produced at rest and particles 2 and 3 
share the remaining energy 2m — mi, as per Eq. (4.2), mutatis mutandis. However, the maximal 
value of El does not correspond to either particle 2 or 3 being being produced at rest. Equation (4.2), 
if rewritten as 

2mEi = Arn^ - + ml , (4.4) 
where TO23 is the invariant mass of the {2, 3} subsystem, indicates that Ei is maximal, with value 

(m2 + m^Y + "^1 



maxi?i 



4m 



(4.5) 



when m23 is minimal, i.e., m23 = m2 + ma, when particles 2 and 3 are at rest relative to each other. 
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4.1.3 Dalitz plot 

Energy conservation, rewritten for the kinetic part Ti = Ei — rrii implies that 

T^i + T2 + Ts = T = 2m - mi - m2 - ms , (4.6) 

remains constant from one event to another This property is fulfilled if Ti is represented by the 
distance of a point to the i*'' side of an equilateral triangle of height T, as in Fig. 9. We just saw that 
while Ti = is possible, Ti = T contradicts momentum conservation, which requires 

\P2 -Psl <Pl <P2+P3 ■ (4-7) 

Saturating this inequality, i.e., fixing the three momenta pi to be parallel, gives the boundary of the 
Dalitz plot, which corresponds to all possible sets {7^} allowed by energy and momentum conser- 
vation. In the non-relativistic limit, the frontier is a circle for identical particles, and an ellipsis for 
unequal masses, as shown in Fig. 9. 



Ai 




Ai 




Figure 9: Definition of 
the Dalitz plot for kinetic 
energies (upper left), and 
boundary in the non- 
relativistic limit, for a decay 
into three identical particles 
(upper right) and for par- 
ticles with masses in ratio 
[4:1:1] (lower left) or [7:4:1] 
(lower right), correspond- 
ing to the mass ratios for 
7] — » Stt, 77' — » Tinn and 
PP — > ri'v'^' respectively. 



In the relativistic case, the shape of the Dalitz plot becomes more angular For the (rare) pp — > 87 
decays, it reduces to a triangle limiting the middles of the sides. Using Eq. (4.4), the Dalitz plot can 
be rescaled to substitute the kinetic energies Tk with the invariant masses mfj , which fulfil 

m^j = 4m^ + ml + ml + ml , (mi + 7712)'^ < mi2 < (2m - m^f . (4.8) 

Also, the equilateral frame is replaced in recent literature by rectangular triangles. Some examples 
are shown in Fig. 10. 
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Figure 10: Dalitz plots with 
relativistic l<inematics: inter- 
pretation of the distances 
{Smfj represents the excess 
of mfj with respect to the 
minimal value (m^ + mj)^), 
decays 77 Sn, pp 37, 
and pp rj'rj-R. 



4.1.4 Multiparticle final state 

The results on the energy range are easily generalised to more than three mesons in the final state. 
If 2m > mi + 777.2 + • • • + rrin, annihilation into the n mesons {771^} is allowed, and the energy of 
particle 1, for instance, is bounded by 

„ 4?7l^ — (777,2 + m^ + • • • + 777,, )^ + 777? 

mi<Ei< — — ^^^-^ — ^ . (4.9) 



4771 



4.2 Phase space 



Heavier particles are less easily produced than lighter ones. Removing this obvious kinematical 
effect leads to a more meaningful comparison among various reaction rates. We give below some 
basic results. For a more detailed treatment, see, e.g., the reviews [237,250]. 

The typical decay rate of protonium, of mass 2777, into a set of 77 mesons {777;} is given by 

/(-) \4 — 371 r ^ 

^=^k^ \^m^'Y-^'ip-p^----p'^)- (4.10) 

i=i * 

If I Alp is removed, one obtains the phase-space integral. For two-body decays (77 = 2), |A^p 
contains the non-trivial dynamical variations when going from one channel to another. For ti > 3, 
|A^P also contains information on the correlation or anticorrelation of momenta, allowing tests of 
scenarios for the formation of resonances, etc. 

For 77 = 2, the phase-space integral is evaluated in the c.o.m. frame [237] as 

^ = Jr^\-M\^- (4-11) 

327r777^ 



For 77 — 3, the phase-space integral is conveniently expressed as an integral over two kinetic 
energies or, equivalently, two invariant 2-body masses, i.e., the variables used to draw the Dalitz 
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plot. The result is [237] 



dr = 



\M\^ 



dTidT2 = 



\M\^ 




(4.12) 



(27r)3l6m 



(27r)3256m3 



where an average over initial spins is implied. Unfortunately, measurements of annihilation at rest 
with a polarised protonium state have not been performed. 

Equation (4. 12) implies that for the fictitious dynamics where \M \ is constant, the Dalitz plot in 
coordinates {Ti, T2} or {m\2i ^^23} is uniformly populated. Peaks in the population density indicate 
the formation of resonances, or at least constructive interferences of several dynamical contributions 
to the amphtude. 

4.3 Conservation laws 

Strong interactions obey very strictly conservation laws due to basic symmetries: energy E, momen- 
tum p, angular momentum J, parity P and charge conjugation C are conserved, as well as flavours. 
Isospin is also a rather good symmetry. G-parity (or isotopic parity) is a combination of charge 
conjugation and isospin, G = C exp(— ^77/2). 

4.3.1 Partial waves 

The algebra of NN quantum numbers is very similar to that of qq reviewed in the section on mesons. 
The partial wave ^^+^'^'^+^Z/j, where L is the orbital momentum and S the total spin has C = 
(-1)^+'^ (if the system is neutral), P = (-1)^+^, and G = (-1)^+'^+^. This is summarised in 
Table 4, for S and P waves. 

Table 4: Quantum numbers of the S and P partial waves (PW) of the NN system. The notation is 

2/+l,2S+l r 



PW 


l.lq 3,lQ 








3,lp^ 


1.3p 3.3p 

-TO -TO 


1.3p^ 


3.3p^ 


l,3p2 


3.3p^ 


jPC 


0-+ 0-+ 


1~ 


1~ 


1+- 


1+- 


0++ 0++ 


1++ 


1++ 


2++ 


2++ 




0+ 1- 


0- 


1+ 


0- 


1+ 


0+ 1- 


0+ 


1- 


0+ 


1- 



4.3.2 Simple rules 

The simplest and most effective selection rule comes from G-parity. Only half of the partial waves, 
for instance, are candidates for annihilation into five pions, those with G = — 1. Once G-parity is 
obeyed, one can generally match any set of quantum numbers by cleverly arranging the spins and 
internal orbital momenta of the mesons. Important exceptions are observed in channels involving a 
small number of spinless mesons or identical particles. We review these selection rules below. The 
reasoning was elaborated in the 60's (see, for instance, [25 1 , 252]) to identify the quantum numbers 
of meson resonances from their observed decays into a few mesons. The large angular acceptance 
and improved resolution of modem detectors forces us to consider higher multiplicities. 

4.3.3 Two spinless mesons 

For NN^ two scalars or two pseudoscalars, only natural parity P = (— 1)'^ is allowed. For n^n^, 
J has to be even; / + J is even for 7t~tt~^, while there is no correlation between I and J for K~K"'". 
For NN — > K°K°, one detects Kg or Ki and CP invariance provides selection rules for KgKg or 
KgKi channels. The results are summarised in Table 5. Similarly, for NN scalar + pseudoscalar, 
only states with urmatural parity P = (— 1)'^+^ contribute. 
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4.3.4 Identical vector mesons 

For two identical vector mesons, the Pauli principle does not lead to further restrictions, once G 
or C conservation is enforced. Total spin S = or 2 of mesons impUes even i and opens the 
= 0+, l"*", 2+, . . . channels, while spin 5=1 requires an antisymmetric space configuration 
and thus = 0^, 1^, 2^, . . . Note, however, that ^ "'Pi uu (or (p(f), or p'^p^) involves a total 
spin S = 2 and an angular momentum £ = 2 in the final state. 

Table 5: Allowed decays from S and P-wave protonium states into selected two-meson final states 
(FS). 



FS 



^(')^(') 
K-K+ 

7rOa;(<^) 



I'^So ^'^So I'^Si 3.3Si l.lpi 3,lp^ l,3p^ 3,3p^ l,3p^ 3,3p^ l,3p^ S.Sp^ 



V 



V V 

V V 



V 
V 

V V 



V 



V 
V 

V 



V 



V 



V 
V 

V 

V 

V V 

V V 



V 
V 

V 

V 

V V 

V V 



V 



V 



4.3.5 Symmetric multi-7r° states 

More delicate is the case of three or more 7r°. This is no longer an academic problem, since these 

channels are seen in modern detectors. For instance, ^'^So has adequate C and G for decaying into 
47r°. One can thus say that ^'^Sq 477° is not forbidden by charge conjugation. To claim that it 
is actually allowed, one should exhibit at least one example of a four-body wave function that is 
together symmetric and pseudoscalar At first sight, this seems impossible. In fact, one can build 
such a wave function, with the desired coupling + i?2 + ^3 = of the internal orbital momenta, 
but with £i, £2 and £3 not vanishing, as it will be explained shortly. As a first step, let us consider the 
7r° + tt" + 7T° case. 

The method adopted below is simple, but a little empirical. We refer to Ref. [253] for a group- 
theoretical treatment. We tentatively write down minimal polynomials in the Jacobi variables with 
the required quantum numbers. Their angular-momentum content is then the lowest term in a sys- 
tematic partial-wave expansion. 



4.3.6 Three tt" 

There is a copious literature on 3-body wave functions and their permutation properties. For instance, 
in the simple quark model of baryons, one should write down spatial wave functions with well- 
identified permutation symmetry, to be associated with spin, isospin and colour wave functions, to 
form a state with overall antisymmetry. Let p be the parity of the orbital wave-function. In the 
harmonic-oscillator scheme [254], the symmetric states are labelled as [56, J^], and the allowed 
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values of in the lowest multiplets with < 3 quanta of excitation are 0+, 1 , 2+, . . . [254]. 
There is no pseudoscalar {J^ = 0~), since to get J = 0, the two internal orbital momenta £i and £2 
should be equal, andthusp = (-1)^1+^2 = +1 x^e states Jp = 1+ and 2" are allowed [251,252], 
but with complicated wave functions, since the coupling of internal momenta, 
vec£2 = J is achieved with ii+ £2 = 8 and 5, respectively. 

To show that = 0+, 1~ and 2~ are allowed for three bosons, it is sufficient to give explicit 
examples. The Jacobi variables 

P = P2-pi, A = (2p3 - pi - P2)/V3 , (4.13) 

are built out of the individual momenta pj. 1 <-> 2 exchange and circular permutation results in 

(A,p) ^ (A,-/5), {-\±V3p,tV3X- p)/2 . (4.14) 

A constant, or + A^, is scalar and symmetric, leading to the existence of = 0+ spatial wave 
functions. The vector {Jp = 1~) 

(A2-p2)A-(2A-p)p*, (4.15) 
is also symmetric. The vector and axial vector 

t7 = (A^ - p2)p + (2A-p*)A, W = pxX, (4.16) 
are both antisymmetric and thus 

V+W+ (4.17) 

has J = 2, = 2, and orbital parity p = —1, and is symmetric. In Eq. (4.17), we use the standard 

notation V+ = + iVy)/\/2. 

For three pions, the parity is P = —p, once the intrinsic parities are taken into account. Hence 
= 0~, 1+ and 2~ are allowed, i.e., ^'^So, ^'^Pi and ^'^P2 can decay into 37r°. 



4.3.7 NN 47r° 

One first eliminates all channels but those with C = G = +1 and thus is restricted to / = 0. 
While it seems = 0+ or 2+ are obviously allowed, it seems at first rather difficult to enforce 
all requirements of permutation symmetry for 0" and 1+. For instance, with the system of Jacobi 
coordinates consisting of x (x p2 — Pi, y pi — and z oc p4 + ps — p2 — Pi, 1 ^ 2 and 3^4 
symmetries are simply translated into an even behaviour in x and ij, but the effect of transpositions 
like 1 -(-^ 3 is not easily written down. The task is simplified by using the variables 

U = {P4 + Pi - P2 - Ps) , 

V = {P4+P2-P3- Pi) , (4.18) 

W = {P4+P3-Pl- P2) , 

since the wave function simply has to survive the changes u -^-^ u, u —v,v <-> w, etc. The 
following wave functions can be written down 





= 0" 


in' 


- v^){v^ - w^){w^ 


— u^)u- {v x w) , 




= 0+ 


1, 








= 1+ 

or 


in' 


- v^){v^ - w^){w^ 


— u^){u xv + vxw + wxu), 






in' 


— v'^){w X u) X (w 


X v) + (i>^ — w'^){u X v) X (u X w) 
+ {vP' - u^){v X w) X {vx u) 




= 2+ 


4 


+ vl + wl. 





(4.19) 
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In summary ^'^So, ^'^Po^ ^'^Pi, and ^'^P2 can decay into four neutral pions. One can probably 
establish these results in a more physical way, by symmetrising amplitudes written as a product of 
terms describing the successive steps of sequential decays, provided there is no cancellation when 
the overall Bose-Einstein symmetry is implemented. For instance, ^'^Sq ^ 'k^&2 is allowed, with a 
relative angular momentum ^ = 2. Also a2 ^ f27r°, with = 1, and in turn, £2 ^ tt^tt". 

4.3.8 Five or more tt" 

For pp — > Stt", ^'^So and the triplet states ^'^Po,i,2 are allowed. We refer here to the paper by 
Henley and Jacobsohn [251]. Symmetric states are found in any Jp, sometimes with many internal 
excitations. In particular, it is difficult to produce five 7r° from = 0~, corresponding to ^''^Po. 

The pattern is seemingly generalisable for any number n > 5 of identical pions. No state is 
strictly forbidden to decay into n. tt", but the transition is sometimes suppressed by the requirement 
of having many internal excitations in the final wave function. 

4.4 Isospin considerations 

We sunamarise here some results on isospin symmetry and its violation in pp annihilation. 

4.4.1 Relations for two-body annihilation 

Simple relations can be written down for two-body decays. Consider, for instance, two-pion events 
with a trigger on a X-ray, ensuring that annihilation takes place from a P-state of protonium. Isospin 
symmetry presumably holds for such pions since their energy is large compared to the tt* — mass 
differences, and Coulomb effects are likely to be negligible in the final state. Then tttt states from a 
= 0+ or 2+ state are pure / = 0, and 

AFp(pp ^ TT+TT") = 2 AFp(pp ^ 7r°7r°) , (4.20) 

where the subscript P indicates annihilation from P-states only. 

4.4.2 Relations for three-body annihilation 

For more than two particles in the final state, there are usually more than one isospin amplitude 
contributing to the transition. Consider for instance 

np K+K°7r" , K+K~7r+ , K°K°7r+ , (4.21) 

for which one can identify two amphtudes, one with (KK) in a / = state, i.e., 

A = ^(np ^ (KK)V+) , (KK)° = ^ ^ ^ , (4.22) 



and another one where (KK) has / = 1 
^1 = ^ 



np ^ ^K+K^Tr^ - ^(K+K-7r+ + K''K\+) 
V 2 2 



(4.23) 



One can extract the contribution of ^0 and Ai to each of the final states and, after squaring, deduce 
(up to a common phase- space factor) that 

AF(K+K-7r+) + AF(KOkV+) - + W , 

2 (4 24) 

AF(K+K°7r°) = . 
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Hence 

AF(K+K-7r+) + AF(K"K°7r+) > AF(K+K°7r") , (4.25) 

the difference being the / = contribution. In other words, one gets often inequalities instead of 
equalities. 

4.4.3 Isospin equalities for tliree-body iinal states 

However, one may sometimes get equalities, either by straightforward Clebsch-Gordan recoupling 
or by more sophisticated methods. An example is antinucleon annihilation into two pions on deu- 
terium. Lipkin and Peskin [255] have shown that 

cr(pd 7r"7r+n) = i a(pd 7r"7r°p) + 2 cr(pd 7r°7r°n) . (4.26) 

This is based on examination of the dependence of the amplitudes on the isospin projection 1^ of the 
outgoing pions. Alternatively, one can write down the amplitudes Ai corresponding to a (tttt) pair 
in isospin I, and obtain (again, to an overall phase-space factor) 

l\Ao\' + l\A,\', 

l\Ao\^, (4.27) 

in the limit of isospin symmetry where, in particular, the initial state has pure / = 1/2. 

Here, we are dealing with integrated cross-sections. More subtle effects can be observed if one 
considers the rate for a given set of relative angles. 



cr(pd — » TT TT+n) = 

cT(pd ^ 7r°7r°n) = 
cr(pd 7r~7r°p) = 



4.4.4 Charge content of iinal states 

As the number of mesons increases, it becomes more difficult to set limits on the relative abundance 
of different charge states. Consider for instance in rijr = 4 pions, pp — > 7r+7r+7r~7r~, 7r+7r~7r°7r° 
and Tr^TT^Tr^Tr". Since intermediate isospin coupling / > 2 are allowed, there is no unique wave- 
function corresponding to the total isospin I = or 1. The relative abundance depends on the 
detailed dynamics. Some general results have however been obtained. See, e.g., Ref. [256]. 
If n+ is the average number of tt"*", etc., with the obvious relations 

n+ + no + ri- = , n+ = n_ , (4.28) 

an / = initial state, which is isotropic in isospin space, will lead to 

n+ = n- = no = n-^/S . (4.29) 

For the isospin 1=1 case, it is found [256] that the ratio 

R{n^) = , (4.30) 

n+ + n_ 



fulfils 



R{2) = 0, l<i?(3)<^, 

n — 2 n — 2 

< Rin-jr) < -r- ^ if tItt > 4 and is even , (4.31) 

4n„ + 2 - ^ - 2n^ + 2 ' ^ ' 

- < R(nTr ) < "^^^ ^ ^ if > 5 and is odd . 
4 - ^ ^ - 2n^ - 2 
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4.4.5 Isospin mixing in protonium 

Isospin symmetry is known to be approximate. Isospin breaking effects are, indeed, observed, such 
as the masses 771(77") and m(7r+) being different, or the transition ^ 7r° + 3/^ not vanishing. 

The long-range part of the protonium wave function is built by the electromagnetic interaction 
and thus corresponds to the simple isospin combination (using our convention) 

, ^ |7 = 1)-|7 = 0) 

|pp) = J ^-^ (4.32) 

If this content had remained the same at short distances, and the strength of annihilation been inde- 
pendent of isospin, pp annihilation would always contain 50% 1 = and 50% 1 = 1. 

However, at short distances, the protonium wave function is distorted by the onset of strong 
interaction. In particular, charged mesons can be exchanged. More generally, the NN interaction 
contains an isoscalar and an isovector part. This latter component induces pp ^ fin transitions. 
Potential model calculations, e.g. [257-260], found that this effect is very large, especially for triplet 
P-states. For instance, it is found that ^Pq is dominantly / = at short distances, i.e., in the region 
where annihilation takes place. 

How reliable are these predictions of potential models? On the one hand, they are based on 
the G-parity transformation applied to the part of nuclear forces which is best established, pion 
exchange. These long-range forces predict a hierarchy of energy shifts that is well observed in 
experiments on protonium spectroscopy, as reviewed in [1]. The hierarchy of widths, in particular 
r(^Po) being larger than other P-state widths, is also a prediction of potential models, and this is 
a crucial ingredient of cascade calculations used to extract the branching ratios from measurements 
done at various values of the density (see Sec. 6). On the other hand, when comparing, e.g., tt^tt^ 
and rjTT^ frequencies, one does not observe the hierarchy one would infer from potential model 
calculations. Hence the problem of the isospin content of annihilation remains open. It will be 
further discussed in Sec. 8. 

The isospin distortion of protonium is described by writing the reduced wave function of a typical 
protonium state as 

*=!M|pp) + !!^|nn), (4.33) 
r r 

with 

f+00 



i\u{r)\^ + \w{r)\^)dr = l. (4.34) 

In potential models, one can solve the coupled (pp, nn) equations, rearrange u(r) and w(r) into 
components uj of given isospin /, to separate the width into its / = and 1=1 components, 

^i = - j \ui{r)\^ 9my(r) dr , (4.35) 

where QmV could depend on isospin 7, though, for simplicity, this was not introduced in early 
optical models. In this approach, the short-range part of the potential, including the annihilation 
component Sml/, is tuned to reproduce the scattering data. 

4.4.6 Isospin content in antiproton-deuterium annihilation 

The problem of the isospin content in pd annihilation has been less frequently discussed. One can 
write the wave function as 

* = u| ppn ) + w\ nnn ) = a\ (NN)7=o n ) -|- 6| (NN) j=i n ) , (4.36) 

where u and w, or a and b, depend on the relative distances of the particles, through a set of Jacobi 
variables x and y. If dr denotes the integration over these variables, then the isospin distortion can 
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tentatively be expressed in terms of the probabilities a and /3 being non equal, where 



a = 




(4.37) 



Note that one cannot exclude a more complicated scenario where the radial profiles of |ap and 
|6p are rather different, in which case the isospin distortion would depend of the part of the wave 
function that is most explored. For instance, one may argue that a decay with two Ught pseudoscalar 
involves a large momentum and thus the short-range part of the wave function, while two heavier 
resonances are produced with small momentum by the external part of the annihilation region. 
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5 Global features of NN annihilation 

Annihilation can be described by a few simple variables of a statistical nature such as the mean 
number of pions and their respective multiplicity distribution, the fraction of events in which an 77 is 
produced or the fraction with strange mesons in the final state. The inclusive momentum distribution 
can be used to argue for a thermodynamic picture of annihilation. Pion interferometry can give in- 
formation about the correlation of pions and hence the size of the fire-ball formed by the annihilating 
proton and antiproton. 

The frequencies of annihilation modes such as pp — > 7r+7r~7r+7r~7r'^ are needed when annihila- 
tion frequencies for intermediate 2-body modes Hke pw are to be determined. In this section we give 
a survey on these global aspects of annihilation. 



5.1 Pionic multiplicity distribution 

In bubble chamber experiments, global annihilation frequencies can be determined by a thorough 
scan of events. Special care is needed to avoid contamination by Dalitz pairs (e.g., from 7r° — > 
7e+e^) faking charged pions. A correction needs to be applied for charge-exchange scattering at 
the end of the antiproton range simulating zero-prong annihilation. Such scans were performed 
for bubble-chamber experiments at Brookhaven and CERN. The scan at CERN gave the following 
results for the charged-particle multipUcity distribution: 

Table 6: Charged-particle multiplicity distribution (from [261]) 

prongs 4.lj:°:2% 

2 prongs 43.2^55;?% 

4 prongs 48. 6^1];^% 

6 prongs 4.lt°J% 



From Table 6 we derive the mean number of charged pions per pp annihilation to be 

n,± = 3.054^0:036- (5-1) 

In events in which the momenta of all particles are measured, a new set of kinematical variables 
can be derived that automatically satisfy energy and momentum conservation. The new momenta are 
then improved in accuracy. In events with all particles reconstructed, four constraints can be used; 
such fits are called fits with four-constraints or 4C fits. Due to the nature of bubble chamber exper- 
iments, only charged particles are detected. But kinematical constraints allow the three-momentum 
of one unseen neutral particle to be reconstructed in a IC kinematical fit. 

The bubble chamber data were split into classes with defined number of charged tracks and fitted 
to different kinematical hypotheses. When no visible Kg was present, the tracks were assumed to 
correspond to pions. A kinematical fit identified with high reliability events without missing parti- 
cles (4C events); events with one missing (IC events) contained up to 12-14% contamination of 
27r0 events. Events not passing the 4C or IC hypothesis were called missing-mass events. The distri- 
bution of pionic states to pp annihilation at rest is shown in Table 7, which is based on publications 
of the Columbia group, on a compilation of (partly unpublished) CERN results [261], and on Crystal 
Barrel data, also partly unpublished. The three frequency distributions are differently normalised: 
the CERN data exclude only events with an detected Kg. The BNL distribution is corrected for all 
annihilation modes containing kaons. Crystal Barrel data are given as branching fractions deter- 
mined from exclusive final states. Hence data containing an r] are not included for r] 77, 7r"'"7r~7 
decays. With an estimated inclusive r] frequency of 7% (see below), about 3.2% of all annihilation 
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Table 7: Annihilation frequencies of pp annihilation at rest in liquid H2 into pionic final states (in units 
of 10"^), from [2,48,215]. Events with more than one n° cannot be reconstructed without neutral 
particle detection; they are listed as missing mass (MM) events. In the Crystal Barrel experiment, 
events with rj —^ 77, 7r+7r-7, w —^ 7r°7 and a fraction of 77' decays do not lead to multi-pionic final 
states, they contribute to the Crystal Barrel missing mass events. Their fraction is estimated in the 
text. Final states with more than 5tv'^ are difficult to reconstruct. Their contribution is estimated by 
assuming that one (^) or two (^) rj mesons were produced and decayed into Stt". Their rate was 
measured from the 77 -» 27 decay mode. 



Final state 


BNL 


CERN 


Crystal Barrel 


all neutral 


32 ±5 


41 


35 ±3 


2^0 






0.65 ±0.03 


3^0 






7.0 ±0.4 


4^0 






3.1 ±0.2 


5^0 






9.2 ±0.4 


67rO (1) 






0.12 ±0.01 


7^0 (1) 






1.3 ±0.1 


SttO (2) 






0.012 ±0.001 


9n° (2) 






0.025 ±0.003 


non-multipion 






15 ±5 




3.2 ±0.3 


3.33 ±0.17 


3.14 ±0.12 




78 ±9 


69.0 ±3.5 


67± 10 


7r+7r-27r° 






122 ± 18 








133 ± 20 


77+77-477° 






36 ±5 


77+77-577° (1) 






13 ±2 


7r+77-MM 


345 ± 12 


358 ±8 


65 ± 20* 


277+277- 


58 ±3 


69 ±6 


56 ±9 


277+277-77° 


187 ±7 


196 ±6 


210 ±32 


277+277-277° 






177 ±27 


277+277-377O 






6±2 


277+277-MM 


213 ± 11 


208 ±7 


30± 15* 


377+377- 


19 ±2 


21.0 ±2.5 




377+377-770 


16 ±3 


18.5 ± 1.5 


1 40 ± 3* 


377+377-MM 


3± 1 


3± 1 




Sum 


954 ± 18 


986 ±6 


970 ± 58 


* Including final states with open strangeness and other non-muWpion events. 
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events are lost. The inclusive uj rate is not known; if one uj is produced in 10 annihilations, there is a 
fraction of 0.9% of all annihilation not leading to one of the final states listed in Table 7. Similarly, 
1.5% of all annihilations are missing if the (unknown) inclusive rj' rate is 2%. Including the 5.4% 
kaonic annihilations (see below), we expect altogether 11% of all annihilations not to contribute to 
the frequencies in Table 7. We assign (educated guess) 6.5% to 2-prong, and 3% to 4-prong and 
1.5% to zero-prong annihilation, respectively. The total Crystal Barrel 6-prong yield is taken from 
the BNL and CERN results. 

We notice a very significant even-odd staggering in the frequencies of multi-7r° final states: 
annihilation into 2n7r° is reduced since their production from S-states is forbidden or suppressed. 
As seen in Sec. 4.3, ^^'^Si and "^'^Si cannot decay into any number of tt*' by charge conjugation 
invariance, "^'^Sq is forbidden by G-parity for any 2n7r° mode; ^'^Sq cannot decay into 27r°, and 
its decays into 47r° or into higher even modes are anyhow suppressed by the internal orbital barrier 
required to match parity and Rose statistics. 

From Table 7 we derive frequencies into multipion final states. They are presented in Table 8 
and compared to the frequency distribution estimated by Ghesquiere invoking arguments based on 
isospin invariance [262]. The frequencies of [262] are normalised to unity; in the central column, 
only multi-pion final states are included and the expected sum of all contributions is 88%. 



Table 8: Pionic multiplicity distribution. 





From Table 7 


From [262] 


2 pions 


0.38 ±0.03% 


0.38 ± 0.03% 


3 pions 


7.4 ±0.3% 


7.8 ±0.4% 


4 pions 


18.1 ±1.8% 


17.5 ±3.0% 


5 pions 


35.2 ±3.7% 


45.8 ± 3.0% 


6 pions 


23.3 ±2.8% 


22.1 ±1.5% 


7 pions 


3.3 ± 0.3% 


6.1 ±1.0% 


8 pions 




0.3 ±0.1% 



The mean number of pions per annihilation into multi-pionic events is estimated to 

= 4.98 ± 0.35 , n^± = 3.14 ± 0.28 , n^o = 1.83 ± 0.21 . (5.2) 

Using Table 6 and the number of photons per annihilation rij = 3.93 ± 0.24 [263] (and allowing for 
small fractions from 77 — > 77 and similar decays), the numbers in (5.2) can be refined to 

= 4.98 ±0.13 , n^± = 3.05 ± 0.04 , n^o = 1.93 ± 0.12 . (5.3) 

The inclusive 77 rate was determined to be 

Tin = 0.0698 ± 0.0079 , (5.4) 

per annihilation [107]. 

Figure 1 1 shows the pion multiplicity distribution from pp annihilation at rest and a comparison 
with a Gaussian fit. The mean number of pions is now 5.03 ± 0.05, the width is 1.13 ± 0.07. Amado 
et al. [264] have shown that a fit with a Poisson distribution constrained by energy and momentum 
conservation reproduces both the mean value and the variance. 
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Figure 1 1 : The pion multiplicity distribution (in %) from Crystal Barrel data, Table 8. 



5.2 Inclusive spectra 

The inclusive momentum spectrum of charged particles (mostly pions) is shown in Fig. 12. The data 
are from the Crystal Barrel collaboration. The distribution reveals no significant structure, except 
for pTT production which identifies itself as a peak at 773 MeV / c in the momentum distribution. The 

signal seems much more pronounced; this is an artefact of the experimental resolution which is 
better for the recoiling 7r° than for tt^ mesons. 

The absence of narrow signals in the momentum spectrum - which would indicate production of 
narrow quasinuclear bound states - is evident. Experiments in the early phase of LEAR confirmed 
the absence of narrow states against which charged or neutral pions would recoil [1 1 1, 265, 266]. 
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Figure 12: The momentum distribution of charged and neutral pions. 



From the number of events in Fig. 12 and the number of annihilation events, the average multi- 
plicities 

= 5.19 ±0.15 , 7v± 3.12 ±0.12 , n^o = 2.07 ± 0.08 , (5.5) 

are found, consistent with those given in Eqs. (5.2) and (5.3). 

The fit in Fig. 12 corresponds to a Maxwell-Boltzmann momentum distribution as proposed 
by Orfanidis and Rittenberg [267]. In the high-momentum range the fit follows the experimental 
distribution adequately thus defining a temperature of about 120 MeV. There is a mismatch at low 
momenta which is more pronounced for charged than for neutral pions. This is not due to a reduced 
efficiency of the detector for low momentum particles. Otherwise, the frequencies (5.5) would be 
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smaller than those in (5.2) and (5.3); rather, pions at low momenta do not follow a simple thermal 
distribution law. 

The temperature of 120 MeV should not be interpreted as an annihilation temperature, nor its 
inverse as an annihilation range of 1 .7 fm. Pion interferometry gives a similar range; we will discuss 
the reasons for this wide range below. 



5.3 Pion interferometry 

In the fireball picture of pp annihilation, pions are emitted stochastically from an extended source. 
It can be argued that the pion fields emitted from different space-time points 1 and 2 superpose and 
interfere. The probabiUty to detect a pion is then given by 



(5.6) 



It is worthwhile to recall that such interferences were first used in astronomy, by Hanbury-Brown 
and Twiss (HBT) [268], to determine the diameter of Sirius from a measurement of the correlated 
intensity fluctuations in two optical telescopes. G. and S. Goldhaber, Lee and Pais [269] applied the 
effect to estimate the size in space-time of the source emitting pions in pp annihilation. 

Different correlation signals have been defined to extract the correlation between pions due to 
the HBT effect. We mention here the two-pion correlation function 

^ P2iPl,P2) 

C{pi,P2) = —, T , (5.7) 

Po(Pl,P2) 

where ^2(^1,^2) represents the distribution of two mesons correlated due to the HBT effect, and 
P2{Pi,P2) an uncorrected sample. The uncorrelated distribution can be chosen as the product of 
two single particle distributions ^2(^1,^2) = Pi{Pi)pi{P2), where 

1 da • 1 r, 

'^^'^^=admj(2E-y ^ = (33) 
1 da 
P2{P1,P2) = -d3pj{2Ei)d3P2/{2E2) ■ 

The correlated sample may be pion pairs of equal charge; for the uncorrelated pion pair, pions of 
different charge or from different events can be used. 

In pp and pd annihilation, the HBT effect was used by different groups. The precise source 
parameters depend on the data (annihilation into 2tt~^2tt~ or 27r+27r^n7r") and the model (with 
unlike-sign pions or different events to determine the uncorrelated dipion spectrum). The results 
given in 1270] and [120] are not compatible with each other within the quoted errors; the value for 
the Bose-Einstein correlation parameter, 

rBE = CTBE = 1.5 ± 0.3 fm , (5.9) 

seems to be a fair estimate of the dimension of the pion source in space and time. It corresponds 
approximately to the pion Compton wave length A,r- Again, this is not the size of the annihilation 
source. Rather, it is the size of the source from which pions are emitted. Pions may be emitted as 
initial- or final-state radiation (even though there is no experimental support for these processes); 
but certainly, pions are produced in secondary decays over a wide range of distances. Equation (5.9) 
indicates that interference between pions is strong for distances corresponding to A^r . This is a trivial 
statement as long as the annihilation range is small compared to A^r- 

The result (5.9) is in rather strong disagreement with recent findings by Locher and Markushin 
using CPLEAR [271, 272] and Crystal Barrel [273] data. These authors go beyond the conventional 
HBT analysis by plotting the double-differential cross-sections as a function of the invariant masses 
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M+_|_ and M , where the sign stands for the charge of the dipion. In a first article, devoted to pp — !■ 

277+ 27r~ , a source dimension of the order of 0.4 fm is quoted, in fair agreement with the value we will 
derive in Sec. 8 from the systematics of two-body annihilation frequencies. Locher and Markushin 
warn the reader that the interpretation of the observed correlation signal by the conventional HBT 
effect is questionable. They observe strong enhancements at low values of -M++ and M__, and 
notice that the signal for pion pairs of large momenta (~ 800MeV/c) - which should be sensitive 
to the source dimension of 0.4 fm - gives a coherence which is unreasonably large within the HBT 
framework. In the analysis of the pp 2it^2it~'k^ reaction, they show that the enhancement 

at low values of M++ and M can be simulated as the effect of the trivial need to symmetrise 

isobar amplitudes. Best suited is an interference between pan and ppn amplitudes. Such an ansatz 
reproduces also the momentum dependence of the correlation signal. The author thus refuse to 
provide simple numbers on the size and life time of the hypothetical pp fire-baU. Also in the case of 
the reaction pp Att^, an interpretation within the conventional isobar model is possible, and there 
is no evidence for an additional signal due to the HBT effect. It may be worthwhile to recall that the 
HBT effect has not been taken into account explicitly in partial-wave analyses of bubble-chamber 
or LEAR data; it is only partly accounted for by Bose-Einstein symmetrisation of the amplitudes. 
There is no parameter to describe the source dimension. 

5.4 Strangeness production 

The full data sets at BNL and CERN were scanned in searches for Kg decays into tt+tt^, leading 
to secondary vertices. These data samples comprised 40,000 and 20,000 events, respectively, with 
strange particles in the final state. In the presence of one Kg, the ionisation density of the tracks was 
used to identify the charged-kaon track. Table 9 gives annihilation frequencies with an observed Kg. 
From Table 9 the fraction of events with at least one neutral Kg is determined to be 



The rate for Ki production is obviously the same. The total yield of strange particle production 
has to include final states with K+K~ pairs. Based on a scan searching for high ionisation-density 
tracks, Armenteros et al. found a contribution of (6.82 ± 0.25)% to pp annihilation. In the CERN 
fist of pionic annihilation modes, these events are included in the missing mass class of events. The 
BNL group assigns (4.6 ± 1.8)% of all annihilations to strangeness production, which agrees with 
the estimate (4.74 ± 0.22)% of Batusov [274] assuming that annihilation frequencies for pp — !■ 
K+K~n7r are similar to those into K°K°n7r. A statistical treatment of these three numbers being 
not plausible, we quote here their linear mean and spread 



In short, one event out of 20 contains strange particles in the final state. 
5.5 Annihilation on neutrons 

Antiproton-neutron or antineutron-proton interactions at rest offer additional opportunities to study 
annihilation dynamics. Obviously, both systems have isospin 7=1. Under the hypothesis that 
annihilation at rest takes place when the NN system is in S-wave, G-parity fixes the total spin: a 
triplet 5=1 has G = (— 1)^+'5+^ = and leads to an even number of pions, while the spin 
singlet system gives an odd number of pions. 

Since there exists no free-neutron target, the cleanest way to study pure isovector annihilation 
is to build a antineutron beam line (produced by charge exchange pp fin) at very low energies. 
The Obelix collaboration has made extensive use of this possibility [11] to study specific reactions 
as a function of the n beam momentum. These results will be presented in Sec. 6. The topological 



BR(pp ^ Kg anything) = (1.55 ± 0.06)% . 



(5.10) 



BR(pp kaons + anything) = (5.4 ± 1.7)% . 



(5.11) 
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Table 9: Frequencies of pp annihilation at rest in liquid H2 into kaonic final states (in units of 10~^) 
from [2]. The data are corrected for unseen decay modes of the Kg and for the Ks reconstruction 
efficiency. K° stands for the sum of Kg and Ki. For the derivation of the KsKitf" frequency see [21 5]. 



Final state 


BNL 


CERN 


K+K- 


1.10 ±0.10 


0.96 ± 0.08 


KgKs + KiKi 


0.010+nin 


0.008 ± 0.008 


KsKi 


0.71 ±0.10 


0.80 ±0.05 


(K,Ks + KiKi)7rO 


1.46 ±0.20 


1.56 ±0.12 




0.67 ±0.07 


0.67 ±0.07 


(K,Ks + KiKi)MM 


1.28 ±0.16 


1.42 ±0.26 


KsK±7r=F 


4.25 ±0.55 


4.25 ±0.20 


(KsKs + KiKi)7r-7r+ 


4.02 ±0.52 


3.90 ±0.46 


KsKi7r-7r+ 


2.41 ± 0.36 


2.26 ±0.45 


KOK±7r=F7rO 


8.94 ± 1.06 


9.38 ± 1.10 


{K,K, + KiKOtt-tt+ttO 


2.98 ± 0.44 


2.20 ±0.28 


K0K±7r=F7r-7r+ 


0.59 ±0.08 


0.71 ± 0.07 


K0K±47r 


~ 


~ 



Sum 28.4 ± 1.5 28.1 ± 1.4 



annihilation frequencies agree within errors with those for pn [275] even though in the latter case, 
one has to worry about the role of the proton or neutron surviving the annihilation process. Global 
features of annihilation were derived by use of a deuterium target. 

Naively, it may be expected that pn annihilation on a deuteron can be viewed at as two-step 
process: 

pd (pn) + p , or pd ^ pp + n , (5.12) 

with subsequent annihilation of the pn or pp system. This is, however, a crude approximation. 
Figure 1 3 shows the neutron momentum distribution for the reaction pd n27r° from the Crystal 
Barrel experiment [177]. The neutron momentum distribution does not follow the Hulthen function. 
There is a significant excess of high-momentum neutrons. The reason can be understood once the 
p7r° invariant mass is plotted, see Fig. 13, right. Baryon resonances are produced, e.g., in the reaction 

pd ^ A°Tr°; A" ^ n7r° . (5.13) 

This atmihilation mode is called a Pontecorvo reaction; the surviving nucleon, proton or neutron, 
can be produced at a rather high momentum. In this case, a high-momentum A(1232) is produced. 
The systematics of Pontecorvo reactions will be discussed in Sec. 6.7.3. 

In the low-momentum part in Fig. 13, the neutron has acquired httle momentum; it acted as a 
spectator and was not involved in the annihilation process. A cut on this momentum, at 200 to 250 
MeV/c in the analysis of bubble chamber data or at about 100 MeV/c in data taken at LEAR, makes 
sure that the annihilation took place on a quasi-free nucleon. 

In bubble chambers, antiprotons annihilating on neutrons lead to odd numbers of visible pion 
tracks; the tracks of 'spectator' protons with very low momenta (less than 80 MeV/c) are not separa- 
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ble from the primary ionisation spot produced by the stopped antiproton. Protons with momenta up 
to 250 MeV/c are easily identified by their large ionisation losses. Thus it was easy to separate anni- 
hilations on quasi-free protons and neutrons. Bizzarri [81] deduced probabilities Sp = 0.571±0.005 
and Sn = 0.429 ± 0.005 for annihilation on protons and neutrons, respectively. Antiprotons stopped 
in liquid D2 annihilate more often on protons than on neutrons, with the ratio 



Sp/s„ = 1.331 ±0.019. (5.14) 



If the isospin content is denoted 



+ 5p , Sn — ^ (5.15) 



we have — sJj/2, thus ensuring model-independent relations due to overall isospin conservation, 
such as 

BR(pd TT-^p) _ ^ 

BR(pd^^Oc.n) " ^' 

Then > s^, more precisely, Sp — (1.662 ± 0.038) Sp, is needed to get the observed Sp/sn- This 
means that in D2, isoscalar pp annihilation is significantly more frequent than isovector annihilation. 

Table 10 lists the global annihilation frequencies for pd annihilation into pions and kaons for 
events with spectator protons. Frequencies for final states with more than one 7r° are from [209,216] 
the frequency for pn 7r~7r° is derived in [209,276]. The Crystal Barrel collaboration [209,216] 
used no spectator cut to determine frequencies. 
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Table 10: Frequencies of pn annihilation at rest in liquid D2 into pionic final states [78] and into final 
states with strangeness [79] after a cut on the spectator momentum of the proton at 250 MeV/c . 



1 lllcll MlILC 


P T'/""/ 1 1 1 n c \/ ( 1 n 1 
1 1 L U LIL 1 Iv. V \ 1 11 /(J J 




-la A-\- f] K 


TT TT 




-7r~9-7r0 
TT ^TT 


U.UO HZ U.U 1 


7r~47r^ 

/I rtH 




^/l /I 


59.7 ± 1.2 


27r-7r+ 


1.57 ±0.21 


27r-7r+7r° 


21.8 ±2.2 


27r-7r+27r° 


7.0 ± 1.1 


37r-27r+n7r° 


23.4 ±0.7 


37r-27r+ 


5.15 ±0.47 


37r-27r+7r0 


15.1 ± 1.0 


47r-37r+n7r° 


0.39 ±0.07 


Sum 


95.5 ± 1.5% 



HI 11 '"1 1 ^ 1 '1 1 P 
1 illcll MlILL 


1- r/Ti 1 ipncx' /in Id '^'^ 
1 ILULlLllL-y yi-ii -'-V ^ 


IV Iv 


14 7 -1- 9 1 
1 ZIZ z . ± 


XV XV 71 


'^fi -1- 4 9 

OU.U ZIZ 




14 7-1-90 

L^. 1 ZIZ Z.U 




91 9 _|_ q c 

Zl.Z ZlZ O.D 


K°K+7r-7r- 


24.8 ±2.6 


KOR-tt+tt- 


34.2 ±3.5 


K,Ks7r-7rO 


25.6 ±2.8 


KOK+Tr-TT-Tr" 


1.6 ±0.9 


KsK-7r+7r-7r° 


33.6 ±3.8 


K^K-w 


35.0 ±5.2 


KsKsTT+TT-TT- 


2.8 ± 1.2 


KsKlTT+TT-TT- 


1.9± 1.2 


Sum 


2.5 ± 0.1% 
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6 Annihilation into two mesons 



6.1 Introduction 

The study of NN annihilation at rest into two mesons is a rich source of information about annihila- 
tion dynamics. The role of symmetries, the topology of the dominant quark diagrams, the violation 
of the Zweig rule, etc., can be inferred from the knowledge of the two-meson frequencies. Detailed 
information on the density dependence of annihilation frequencies is required if these frequencies 
are to be assigned to specific states of the pp or pd atom. This arises from the sensitivity of two- 
meson annihilation to the initial state of the nucleon-antinucleon system. For example, from Table 
5, we see that the process pp tt+tt" takes place from initial S- and P- states whilst the reac- 
tion pp — > TT^TT^ only originates from P-states. At first glance, it would seem that the fraction of 
P-state annihilation could be obtained from the relevant tt+tt" and tt^tt" annihilation frequencies 
using isospin invariance and simple arithmetic. The derivation of the P-state fraction is actually more 
subtle and unfortunately, more complex. In particular it depends on details of the atomic cascade 
process and on the amount of Stark mixing. On the other hand, a thorough analysis of these cascade 
effects enables us to obtain the fraction of pp annihilation into a variety of channels from different 
fine-structure states of the pp system. The information which is directly given by the tt+tt" and 
tt^'tt" frequencies is the fraction of annihilations for the tt+tt" and tt'^tt" channels which take place 
from S- and P-states. This is not the same as the fraction of S- and P-state annihilations for the pp 
system. 

We first start with a few definitions. The annihilation frequency AF (ch, p), often also called the 
branching ratio in the literature, is the probability that a particular channel ch will be produced in a 
pp or pd annihilation at rest. It is usually the quantity measured experimentally and is a function of 
the target density p. We restrict the term branching ratio BK{ch, "^^^^Lj) to the probability that the 
channel ch is produced by an annihilation from the initial state '^^^^Lj of the pp or pd atom. These 
branching ratios are independent of the atomic physics effects occurring during the cascade of the pp 
or pd atom and hence do not depend on the target density. It is this quantity which normally should be 
compared with predictions from theoretical models. In Sec. 7, we shall define dynamically corrected 
branching ratios, DR, including phase space corrections, orbital-angular-momentum-barrier effects, 
and a phenomenological factor favouring the production of high-mass mesons. 

Many of the earlier experiments measured two-body annihilation frequencies by the detection 
and reconstruction of a single it^ or with a small solid-angle detector and observed peaks in the 
resulting tt" or rj inclusive momentum spectrum to identify the second meson (see Table 3). Since 
the annihilation frequencies are small, these early data are often statistically weak or subject to very 
considerable uncertainties in the background subtraction from the inclusive spectra. In more recent 
experiments, using detectors covering a large solid angle, both of the mesons in the two-body events 
are fully reconstructed and any background is much reduced. 

At a given target density p, the annihilation frequency for the channel ch is given by 

where Cdet is the detection efficiency given by 

_ Nj^gjch) 

-''tot 

Ntot is the total number of antiprotons stopped in the target, corrected for pile-up, and A^ewts (ch) 
is the number of events, corrected for background, attributed to channel ch. N^f'' and N^^{ch) 
are the corresponding numbers for the Monte-Carlo simulation and reconstruction. V is the fraction 
of antiproton events which stop in the target volume and Fi and F2 are the correction factors for 
the probability of the observed decay for each of the two mesons in channel ch; for example F = 
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0.3943 ± 0.0026 for r] — > 27 in the case where the two photons are measured so as to reconstruct 
the 77 meson. 

Here we will particularly discuss two-body annihilation frequencies for channels involving nar- 
row (T < lOMeV) mesons (tt, K, 77, t]' , oj and (6). Measurements of annihilation frequencies for 
broader states from inclusive spectra usually have considerable uncertainties due to the need for 
background subtraction. In addition, reflections can distort the spectrum in an unpredictable way. 
For measurements of two-body annihilation frequencies of broad states a full partial wave analysis 
involving all intermediate states is required. For completeness, we include two-body branching ra- 
tios involving the K* (F = 50MeV) and p (F = 151 Me V) mesons. Some two-meson annihilation 
frequencies obtained from partial wave analyses involving pseudovector or tensor mesons are also 
discussed. 

In measurements of frequencies for two-body annihilation from pp atoms, the selection of events 
is straightforward, requiring just two back-to-back particles each of the same well-defined momen- 
tum. For pd atoms the selection is less straightforward due to the presence of the recoil nucleon 
which, in some cases, can have momenta up to 1.2 GeV/c [84, 209]. Measurement of their frequen- 
cies requires that all events over the full range of spectator momenta are included. An additional 
effect is that the spectator momentum distribution depends [84] on whether the annihilation occurs 
from an atomic S- or P-state of the pd system and whether the p is in an S- or P-state relative to 
the nucleon with which it annihilates. In particular the momentum distribution for P-state pN anni- 
hilation from an atomic S-state is particularly broad. It was pointed out [277] that tight cuts on the 
momentum of the spectator nucleon therefore suppresses P-state annihilation from atomic S-state 
orbitals. 

For pd annihilations at rest, the annihilation frequency is determined by normalising either to the 
total number of annihilations on a deuteron or to the number of annihilations on a proton or neutron. 
In the following, subscripts d and N (p or n), will be used to distinguish between frequencies mea- 
sured for annihilation on deuterons or on a nucleon (proton or neutron) in the deuteron. Sometimes, 
it has been assumed that there is equal probability for annihilation on a proton or neutron in the 
deuteron. In this case the branching ratio for pd annihilation would be obtained simply by dividing 
by a factor of two. Here, we use the probabilities Sp — 0.571 ± 0.005 and Sn = 0.429 ± 0.005 de- 
duced by Bizzarri [81], as discussed already in Sec. 5.5. The annihilation frequency on a deuteron, 
for reactions involving a spectator neutron, is then given by 

AFd (c/i, n, p) = Sp AFp (c/i, n, p) , (6.3) 

and similarly for reactions involving a spectator proton. 

In Sec. 6.2, we discuss the pp and pd atomic cascade, the effects of Stark mixing, the role 
of the fine-structure levels and the need for enhancement factors. Atomic cascade calculations, 
and in particular the prediction of enhancement factors is discussed in Sec. 6.3 for both pp and 
pd atoms. Pubhshed annihilation frequencies for pp atoms are reviewed in Sec. 6.4 which also 
contains a detailed discussion of the pp— ^tt'^tt'^ annihilation frequency. A new analysis of two-body 
annihilation frequencies for pp atoms is presented in Sec. 6.5 with an emphasis on determining 
the fraction of P-state annihilation and branching ratios for annihilation from specific atomic states. 
Section 6.6 reviews antineutron annihilation on protons whilst section 6.7.1 reviews annihilation 
frequencies from pd. The determination of the P- state fraction in pd annihilation is presented in 
Sec. 6.7.2. A compilation of data on Pontecorvo reactions is given in Sec. 6.7.3. Some final remarks 
on aspects of two-meson annihilation are made in Sec. 6.8. 

6.2 Stark mixing and density dependence of annihilation frequencies 
6.2.1 pp atoms 

Formation of the pp atom and its atomic cascade has been discussed in the earUer review [1] and 
elsewhere [278, 279]. Briefly, the capture of the p typicaUy occurs at a principal quantum number 
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n « 30 [2801 . De-excitation then takes place by a number of processes including radiative transitions 
with the emission of X-rays, and the external Auger effect involving the ionisation of a neighbouring 
H2 molecule. Finally the p reaches an atomic state with angular momentum i = or 1 when 
annihilation occurs. Annihilation from states with L > 2 can be ignored due to the negUgible 
overlap of p and p in the atomic wavefunction. 

In addition, except at very low target densities, the Stark effect gives mixing of the angular 
momentum states L at high n allowing the protonium atoms to transfer to S- and P-states where 
they can annihilate before reaching the low-n states. The Stark mixing increases with the target 
density; for Uquid targets the rate is very high. When considering frequencies for pp annihilation, 
it is important to consider the effects of the fine-structure of the atomic states. For pp atoms, the 
states with L < 2 are ^Sq, ^Si, ^Pi, ^Pq, ^Pi and ^P2 with the corresponding J^^'= 0^+, 1 , 
0++, 1++ and 2++. Table 11 shows the predictions of Carbonell et al. [259] for the widths of 
these states. They were obtained using potentials for the pp interaction due to Dover and Richard 
(DRl and DR2) and Kohno and Weise (KW). A particular feature of these predictions is the very 
liirge width of the ^Po state which has been confirmed in recent experiments. Gotta et al [281], 
assuming the widths of the ^Pi, ''Pi and ■^P2 states to be equal, obtain r(^Po) = 120± 25meV and 
r(^Pi, "^Pi, ■^P2) = 30.5 ± 2.8 meV for the states with principal quantum number n = 2, in good 
agreement with the values listed in Table 1 1 . The average value of the width of the IS state given by 
four experiments is quoted in the earlier review [1] as r(lS) = 1.060 ± 0.080 keV which again is in 
good agreement with the values of Table 1 1 . 

Table 1 1 : Widths for pp atoms as predicted by three potential models 



State 








^Po 


^Pi 


^P2 


Units 


keV 


keV 


meV 


meV 


meV 


meV 


DRl 


1.02 


0.90 


26 


114 


20 


30 


DR2 


1.04 


0.92 


28 


80 


18 


32 


KW 


1.26 


0.98 


26 


96 


22 


36 



At high n, in the case where Stark mixing is important, the fine-structure levels are continually 
and rapidly repopulated according to their statistical weight. A fine structure level with a large 
annihilation width will therefore contribute more to annihilation than would be expected from its 
statistical weight only. This effect is particularly important for the ^Pq level. Similar, but smaller, 
effects will also occur for the other fine structure P-states. The effects for S-states are typically less 
than 5%. 

These deviations of the population of the fine-structure states have been described [282] in terms 
of enhancement factors E(^^~^^Lj, p) which are functions of the initial state '^^'^^Lj and target 
density p. Values of Ei^^^^Lj^p) < 1 (> 1) correspond to a fraction of annihilations smaller 
(larger) than that expected on the basis of a purely statistical population of the level. 

The annihilation frequency AF(c/i, p) can then be written [282] in terms of the branching ratios 
BR(c/i, 2^+^X7) in the form 



AF(c/i,p) = (1 -fp(p)) [ii?(iSo,p)BR(c/i,iSo) + \e{^Su p)B-R{ch,^Sx) 



fp(p) 



AiJ(ipi, p)BR(cft, iPi) + ^i^('Po, p)BR(cft, 3Po) 



(6.4) 



+ A£;(3p^,p)BR(c/l,3Pi) + A£;(3p2,p)BR(c/i,3p2) 



where fp(p) is the fraction of P-state annihilation and the factors 1/4, 5/12, etc., are the statistical 
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weights of the states. In those cases where production of the channel ch is forbidden from an initial 
state ^^~^^Lj due to selection rules, then BR{ch, "^^^^Lj) = 0. 
The enhancement factors are normaUsed [282] so that 

^ECSo,p) + ^ECSi,p) = 1, (6.5) 

and 

^ECPi,p) + ^EePo,p) + ^ECPi,p) + ^EeP2,p) = 1 . (6.6) 

From Eq. (6.4) it can be seen that the density dependence of AF (ch, p) arises from two factors. 
The first and larger one is the density dependence of the fraction of P-state annihilation, fp (p), which 
is directly due to the Stark effect. The second and more subtle one arises from the enhancement 
factors i;(2^'+iL.7,/3). In cases where these differ significantly at a given p for various °^ Lj, the 
branching ratios for these states, BR(c/i, ^■^+^Lj) can be determined. 

Attempts have been made to determine the values of the enhancement factors from experimental 
data. A study [198] of the reaction pp rjn^n^n^ at rest, both in liquid and in gas at 12 psTP> 
where STP indicates Standard Temperature and Pressure, finds 

to be compared with the value r w 1.7 obtained using predicted values for the enhancement factors 
from Table 12 which will be discussed later. 

Salvini et al. [283] have attempted to determine enhancement factors directly from a best fit to 
two-body annihilation frequencies. The fit procedure was not straightforward since they obtained 
different local minima with more or less equivalent values. This is probably due to the fact that 
the values of BR(c/i, ^'^+^Lj) and Ei^^+^Lj, p) are strongly correlated, appearing in the form 
Ei^^+'^Lj, p)BR(c/i, in Eq. (6.4). The solutions were constrained to those that were re- 

garded as physically meaningful, the remaining parameter space was explored in detail. In most 
cases values for Ei^^^^Lj, p) were obtained consistent with unity except for Ei^Po, liq.) which 
was in the range 1.7 ± 0.3 to 3.6 ± 0.1, depending on the value for AF(7r°7r°, liq) used in the data 
analysis. The discrepancies between the different experimental values for the pp — > 7r°7r° annihila- 
tion frequency in liquid H2 will be discussed in Sec. 6.4.3. 

More recently Bargiolti et al. [284] have determined enhancement factors E{^^^^Lj, p) for all 
S- and P- fine-structure states in liquid and gas (psTp) targets by fitting measurements for two-body 
and and also for resonant two-body channels obtained from partial wave analyses of the channels 
7r+7r~7r'^, K^'^K^tt^ and K+K~7r'^ in hydrogen targets at three different densities {p = O.OOSpsTP, 
PSTP and liquid). They fix the enhancement factors at low density O.OOSpsTp) = 1. 

With this assumption they were able to obtain a good fit to the measurements with enhancement 
factors different from one for the ^Pi, ^Pg and ^Pi levels in liquid hydrogen . Comparison with 
the predictions of cascade calculations (see Table 12) shows good agreement for the ^Pq level but 
disagreement for the ^Pi and ^P^ levels with the ^Pi level, £;(iPi,liq.) = 0.09 ± 0.17, being 
particularly strongly suppressed. The reason for this discrepancy is not understood. 

An alternative approach [282] is to calculate values of Ei^^^^ Lj , p) using an atomic cascade 
calculation [278] for the pp atom whose parameters are obtained by fits to pp atomic X-ray data. 
This method will be discussed in detail in Sec. 6.3. 



6.2.2 pd atoms 

The general features of the atomic cascade for pd atoms are very similar to those for pp atoms dis- 
cussed above, except that the p is now captured at a principal quantum number n « 45 and that 
the much larger predicted widths [285] for the 2P levels in pd atoms give rise to increased proba- 
bilities for annihilation from atomic P-states. D-state annihilation also becomes significant. Wycech 
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et al. [285] predict that the widths of the fine-structure components for P-states are approximately 
equal. It is therefore to be expected that for pd atoms the enhancement factors will have values 
E(^^+'^Lj, p) ~ 1 and their effect can be neglected. 

Due to the Fermi motion of the nucleons in deuterium, the p-nucleon angular momentum can 
be different from the angular momentum of the p with respect to the deuteron centre of mass. Even 
when the p is in an atomic S-state, it may be in a P-state relative to the nucleon with which it 
annihilates. This process, sometimes called "induced" P-state annihilation, was first considered by 
Bizzarri et al. [84] and has recently been discussed by Bugg [286]. It is an assumption of these 
models that the antiproton interacts with only one of the two nucleons in the deuteron, the other 
behaving as a "spectator". As a consequence, the annihilation frequencies in deuterium are directly 
related to the frequencies for capture on neutrons and protons. For the remainder of this article we use 
the following notation: capital letters denote the initial pd atomic states whilst the subscript "ann" 
(e.g., Pann) IS uscd to dcscribc the annihilation process. This differs from the notation of Bugg [286], 
which was also used by Batty [287], where small letters denote the initial pd atomic states and capital 
letters were used to describe the annihilation process. For the pp system the fraction of annihilation 
from P-states of the pp atom fp (p), and the fraction of P-state pp annihilation fp^^^ (p) are, of course, 
the same. 



6.3 Cascade calculations 
6.3.1 pp atoms 

In the calculations of enhancement factors [282], measurements of atomic X-ray yields for pp atoms, 
covering the range of target densities from 0.016 to 10.0 psTP> were analysed using a cascade calcu- 
lation based on the method of Borie and Leon [288]. In this method the effects of Stark mixing are 
calculated using an impact parameter technique in which, in its interaction with the electric field of 
the neighbouring H2 molecules, the exotic atom is treated as moving along a straight line trajectory. 
If an atomic S-state is involved, the removal of the degeneracy due to the energy shift caused by 
the strong interaction hinders the Stark mixing, since the electric field must overcome the energy 
difference between S- and P-states. This is taken into account by using a smaller impact parameter. 
Because of the uncertainties in the absolute rate for Stark mixing, an overall normalisation parameter 
fcsTK is usually used and its value determined by fitting the X-ray yields. 

An alternative method, usually referred to as the "Mainz" model, has been developed in [277, 
289]. Here the colUsions of the pp atoms with the neighbouring H2 molecules are simulated using 
the classical-trajectory Monte-Carlo method which only allows Stark mixing when the pp atom 
experiences strong electric fields during collisions with neighbouring hydrogen molecules. This 
avoids the uncertainties associated with the straight fine approximation used by Borie and Leon [288] 
at the expense of a considerable increase in computational complexity. Further information about 
cascade calculations is given in earlier reviews [1,278]. 

Once the parameter fcsTK in the Borie-Leon model has been adjusted to fit the X-ray data [282], 
the two models give similar predictions for the overall variation in the X-ray yields as a function of 
target density. However the Borie and Leon model predicts a significantly larger fraction for P-state 
annihilation in pp atoms than the "Mainz" model as shown in Fig. 14. It has been shown [282] that 
the introduction of an additional parameter ii'o, as a further normalisation factor to the rate of Stark 
transitions between atomic S- and P-states, gives significantly improved fits to the X-ray yield data 
for pp atoms and much reduced values for the fraction of P-state annihilation in pp atoms. Setting 
Kq = 1 gives the usual form of the Borie and Leon model. 

Values of E(^^^^Lj, p) for P-states calculated [282] using the Borie-Leon cascade model with 
the DRl annihilation widths of Table 11 are shown in Fig. 15. The corresponding enhancement 
factors for the S-states are generally close to 1. Table 12 gives numerical values for the enhance- 
ment factors, for all target pressures at which annihilation frequencies have been measured, for the 
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Figure 14: Fraction of P-state annihilation predicted by the Borie-Leon model [288] giving a best fit 
to the X-ray data with Ko = 1 (dashed line) or Ko = 7.6 (full line). Also shown are the predictions of 
the Mainz [277,289] model (dotted line) and the values derived in the present work from annihilation 
frequencies and annihilation widths from the DR1 potential 

three annihilation models of Table 1 1 . The errors quoted are those due to the uncertainties in the 
parameters Kq = 7.6 ± 2.6, fcsTK = 1-19 ± 0.06 used to fit the X-ray yield data for pp atoms. 

6.3.2 pd atoms 

Calculations of the atomic cascade for pd atoms follow those made for pp atoms but with a starting 
value for n = 45 and with strong interaction widths appropriate for pd atoms. Recently proposed 
values for the energy shift and width of the IS ground state [290] AE'is = -1050 ± 250 eV, Tig = 
1 100 ± 750 eV and for the width of the 2P-state [28 1 ] Tzp = 489 ± 30 meV were used (see, however, 
also [1, page 251]). This latter value is in reasonable agreement with the average 2P width = 
422 meV calculated by Wycech et al. [285]. In the case of pd atoms, D-state annihilation has to 
be considered and as measurements of the width of the 3D state are not available, the value Fsd = 
5 /LteV from the calculations of Wycech et al. [285] was used. 

In almost all cases, the only available yield data for pd atoms is for the L X-ray hnes (nD - 2P) 
at eight target densities from 0.016 to 10.0 psTP, (where psTP is the density of D2 gas at STP). As a 
result, it is only possible to determine the value of fcsxK from the fit to the X-ray yields. Data for the 
K-hnes (nP-lS) is required if the value of Kq is to be determined. Equally good fits to the pd X-ray 
data are obtained [287] both with Kq = 1 (Borie-Leon model) and with Kq = 7.6 as determined 
from fits to pp X-ray yield data. 

The calculated fractions of annihilation, fp(p) and inip), from atomic P- and D-states as a 
function of D2 target density, are plotted in Fig. 16 for the best fits with Kq = 1 and Kq = 7.6. The 
values of fp (p) for if = 7.6 are seen to be significantly reduced at higher target densities compared 
to those obtained with Kq = 1. This is in general accord with the results for pp atoms 1 282], shown 
in Fig. 14. Also plotted in Fig. 16 are the calculated values [277, 289] obtained using the "Mainz" 
model. As in the pp case shown in Fig. 14 these values of fp (p) are seen to be significantly smaller 
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Table 12: Calculated enhancement factors, as a function of the density p (in units of pstp). 



State 


p = 0.002 


p = 0.005 


p=l 


p = 12 


p = 15 


p = 27 


Liquid 


Model 




1.044 


1.046 


1.020 


1.011 


1.012 


1.013 


1.032 


DRl 




±0.001 


± 0.001 


± 0.002 


± 0.001 


± 0.001 


± 0.001 


±0.001 


ic 
oo 


0.973 
±0.001 


0.975 
± 0.001 


1.007 
± 0.001 


1.005 
± 0.001 


1.005 
±0.001 


1.007 
±0.001 


1.028 
±0.001 


DR2 




1.138 


1.139 


1.030 


1.019 


1.020 


1.022 


1.060 


KW 




± 0.001 


± 0.001 


± 0.003 


±0.001 


±0.001 


± 0.001 


± 0.001 




0.985 


0.985 


0.993 


0.996 


0.996 


0.996 


0.989 


DRl 




± 0.001 


±0.001 


±0.001 


±0.001 


± 0.001 


±0.001 


± 0.001 


3c 
Ol 


1.009 

±0.001 


1.008 
± 0.001 


0.998 
± 0.001 


0.998 
± 0.001 


0.998 
±0.001 


0.998 
±0.001 


0.991 

±0.001 


DR2 




0.954 


0.954 


0.990 


0.994 


0.993 


0.993 


0.980 


KW 




± 0.001 


±0.001 


±0.001 


± 0.001 


±0.001 


± 0.001 


± 0.001 




0.999 


0.999 


0.974 


0.968 


0.966 


0.958 


0.856 


DRl 




± 0.001 


± 0.001 


± 0.002 


± 0.002 


± 0.002 


± 0.002 


± 0.005 


Ip 


1.000 

±0.001 


1.001 
± 0.001 


0.991 

±0.001 


0.992 
±0.001 


0.991 
±0.001 


0.988 
±0.001 


0.933 
±0.003 


DR2 




0.998 


0.997 


0.960 


0.945 


0.941 


0.928 


0.809 


KW 




± 0.001 


±0.001 


± 0.003 


± 0.002 


± 0.003 


± 0.003 


± 0.004 




1.011 


1.016 


1.288 


1.372 


1.399 


1.487 


2.556 


DRl 




± 0.001 


± 0.002 


± 0.023 


± 0.016 


± 0.017 


± 0.021 


± 0.046 


Po 


1.010 

±0.001 


1.014 
± 0.001 


1.206 
±0.015 


1.280 
±0.011 


1.302 
±0.012 


1.370 

±0.016 


2.076 
±0.027 


DR2 




1.009 


1.013 


1.227 


1.296 


1.318 


1.388 


2.176 


KW 




± 0.001 


±0.001 


±0.017 


±0.012 


±0.013 


±0.016 


± 0.032 




0.995 


0.993 


0.929 


0.894 


0.886 


0.862 


0.685 


DRl 




± 0.001 


± 0.001 


± 0.005 


± 0.004 


± 0.004 


± 0.005 


± 0.005 


^Pl 


0.993 
±0.001 


0.990 
± 0.001 


0.914 
± 0.006 


0.866 
±0.005 


0.856 
±0.005 


0.826 
±0.006 


0.641 
±0.004 


DR2 




0.995 


0.993 


0.932 


0.899 


0.891 


0.868 


0.703 


KW 




± 0.001 


±0.001 


± 0.005 


± 0.004 


± 0.004 


± 0.005 


± 0.004 




1.001 


1.002 


1.000 


1.008 


1.009 


1.011 


0.964 


DRl 




± 0.001 


±0.001 


±0.001 


±0.001 


±0.001 


± 0.001 


± 0.004 


^P2 


1.002 
±0.001 


1.003 
± 0.001 


1.016 
± 0.001 


1.029 
± 0.001 


1.032 
±0.001 


1.038 

±0.001 


1.041 
±0.001 


DR2 




1.002 


1.003 


1.019 


1.035 


1.037 


1.045 


1.058 


KW 




± 0.001 


±0.001 


± 0.001 


± 0.002 


±0.001 


±0.001 


± 0.001 
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at the higher target densities than those predicted by the modified Borie-Leon model [282, 288] with 
Ko = 7.6. 

Table 13 gives enhancement factors for P-states of the pd atom calculated from this cascade 
calculation with a value Kq = 7.6, together with the calculated P-state widths of Wycech et al. [285]. 
Even in the case of the state, which has a width 20% larger than the spin-averaged value, the 
enhancement factor is only an 11% effect in liquid deuterium. In view of the relative smallness of 
this effect and so as to make a significant simplification in the following discussion, for pd atoms 
the enhancement factors are set to 1 . Using Eq. (6.4) the annihilation frequency for production of 
the two-body channel ch in a D2 target of density p is then given by: 

AFn(c/i, N,, p) = [BR(c/i, ^So) + 3BR(c/i, ^Si)] 

f n ^^-^^ 

+ ^sti^El [3BR(c/i, ^Pi) + BR{ch, ^Pq) + 3BR(c/i, ^Pi) + 5BR(c/i, ^Pa)] , 

where Ng is the spectator nucleon. Note that AF^{ch, N^, p) is the frequency for annihilation on a 
single nucleon. 

Table 13: Calculated enhancement factors for pd atoms 



State Width Enhancement factor 



25+1^^ 


meV 


PSTP 


Hq. 


^Pl/2 


398 


0.99 


0.98 


^P3/2 


386 


0.98 


0.96 


^Pl/2 


512 


1.04 


1.11 


i^3/2 


430 


1.00 


1.01 




420 


1.00 


1.00 
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Figure 15: P-state enhancement factors from a cascade calculation [282] using Ko = 7.6 and annihi- 
lation widths from the DR1 potential 




Figure 16: Fractions of annihilation from atomic P- and D-states of pd atoms predicted by the Borie- 
Leon model [288] giving a best fit to the X-ray data with Ko = 1 { dashed line) or Ko = 7.6 (full line). 
Also shown are the predictions of the Mainz [277, 289] model (dotted line). 
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6.4 Frequencies for pp annihilation at rest 

In this section, we present a compilation of annihilation frequencies. Measurements in which peaks 
in the inclusive recoil energy spectrum are used to identify the second meson, will be indicated 
by R (recoil). Measurements in which both of the two mesons in the two-body final state are 
fully reconstructed, which are expected to be more reliable, will be denoted by E (event). In some 
cases [167,215], for all neutral channels, the absolute annihilation frequencies are obtained by nor- 
malisation relative to the frequency for another channel, usually pp — > Tr^Tr*^. These measurements 
will be indicated by N (normalisation). For the majority of entries in the tables, the detection tech- 
nique is also identified. A, C and O are used to identify the Asterix (see Sec. 2.4.1), Crystal Barrel 
(see Sec. 2.4.3) and ObeUx (see Sec. 2.4.2) detectors at the LEAR facility. Other LEAR experiments 
(see Sec. 2.3), in particular CPLEAR are indicated in the tables by C , whilst JC are experiments 
carried out at KEK (see Sec. 2.2). B is used to indicate early bubble experiments (see Sec. 2.1), 
whilst T indicates experiments carried out using electronics detection techniques with early sepa- 
rated antiproton beams (see Sec. 2.3). The notation scheme is summarised in Table 14, for easier 
reading of the following tables. 

In a few cases, the frequency has been measured for the annihilation channel (tt+tt", K+K~, 
KgKg, (f>TT, (prj, (f)U) and p'^ip) in coincidence with L X-rays, so that annihilation only occurs from 2P 
states. For these measurements, the frequency AF(c/i)x, no longer depends on the target density 
and following Eq. (6.4) can be written in the form 

AF(cft)x = ^BR(c/i, iPi) + ^BR(cft, ^Pq) + ^BR(cft, ^Pi) + ^BR(c/i, ^Pa) (6.9) 

It has been confirmed [282] that enhancement factors for 2P-state annihilation, as observed in exper- 
iments using an L X-ray trigger, are within ±1% of the value 1.0 assumed in the derivation of this 
equation. 

The labels S or P, in Table 14 and the following tables, indicate the results of a partial-wave 
analysis. 

6.4.1 Frequencies for annihilation into two charged mesons 

The measurement of frequencies for annihilation into two charged mesons is much easier experimen- 
tally than for channels involving neutral mesons where the problems of measurement and identifica- 
tion are more challenging. This is reflected in the results in Table 15, where there is good agreement 
between the values obtained from different experiments and experimental techniques. The weighted 
mean value of the annihilation frequency for pp tt+tt" from 5 experiments is AF(7r+7r~ , liq.) = 
(3.27 ± 0.05) X 10~^ with x^/^ = 3.01/4, where N is the number of degrees of freedom. For 
the pp K+K- channel, the weighted mean AF(K+K-, liq.) = (1.03 ± 0.03) x lO'^ with 
X^/N = 2.28/3. For both channels the measured annihilation frequencies vary relatively smoothly 
as a function of target density [210, Figs. 7(c) and 7(d)]. 

6.4.2 Frequencies for annihilation into two neutral strange mesons 

The measurements of frequencies for annihilation into two neutral kaons are shown in Table 16. The 
five measurements for pp KgKi in hquid H2 are in moderate agreement with each other and give 
a weighted mean value AF(KsKi, liq.) = (7.86 ± 0.40) x 10"^ with x^/N = 8.45/4, where N is 
the number of degrees of freedom. The two values obtained by the Crystal Barrel experiment [210] 
for this same channel in a gas target at 12psTP disagree by 2.4cr. The lower value seems to be in 
better agreement with the overall trend of results for the KgKi annihilation frequency as a function 
of target density [210, Fig. 7(b)]. 

The annihilation frequency for pp KgKg is small (^ 10~^) and hence difficult to measure. 
Values in a gas target (pstp). both with and without an atomic L X-ray trigger, have been obtained 
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Table 14: Notation scheme for the experimental results: name of the experiment and method of 
analysis 





notation 


meaning 




E 


event 






event with minimum bias trigger 




-L'an 


event with all neutral trigger 






^ K+K- 






<j) ^ K,Ki 


>^ 


F 


a; 7r"7 


S 


F 






Ekc 


Kg 7r"^7r~ 


ta an 


Ekii 






N 


normalisation with respect to another channel 




P 


P-wave annihilation 




R 


recoil 




S 


S-wave annihilation 




X 


X-ray trigger 




A 


Asterix experiment at LEAR 




B 


Bubble Chamber experiments 


g 


C 


Crystal Barrel experiment at LEAR 


& 

X 


K 


Experiments at KEK 


u 


C 


Other LEAR experiments (CPLEAR, . . .) 




O 


Obelix experiment at LEAR 




T 


Eleclronics experimenls using separated p beams 



by the Asterix experiment. Values for hquid H2, where the annihilation frequency is even smaller, 
have only been determined on the basis of four KgKg events observed in the bubble chamber data at 
BNL and CERN. 

6.4.3 Annihilation frequencies for pp — > tt^tt". 

As can be seen from Table 17, there are widely differing measurements for the 7r°7r° annihilation 
frequency in liquid targets, particularly so for the Crystal Barrel [167, 210] and ObeUx [158] exper- 
iments. Whilst the Obelix measurement agrees with some earlier data [106,291] where only 2 or 
3 photons were observed, it should be noted that Obehx, Crystal Barrel and the work of Devons et 
al. [89] are the only experiments to measure all 4 photons and to reconstruct events fully. The 
fact that the tt+tt" annihilation frequencies obtained by the Crystal Barrel and Obelix experiments 
are in good agreement, suggests that the origin of the problem hes in the determination of the re- 
construction efficiency for detecting -k^-k^ events. This has been discussed recently by the Crystal 
Barrel collaboration in a comment [214] on the ObeUx paper [158]. It was pointed out that there 
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Table 1 5: Frequencies for annihilation into two charged mesons 



Chan. Dens. Ann. Freq. Ref. Type 


Dens. Ann. Freq. Ref. Type 


Liq. (3.20 ± 0.30) x lO'^ [42] B E 
Liq. (3.07 ± 0.13) x lO'^ [167] C E 
n+n- Liq. (3.31 ± 0.06) x lO'^ [158] O E 
1 (4.30 ± 0.14) X 10-3 [113] >IE 
0.005 (4.26 ± 0.11) X 10-3 [134] O E 


Liq. (3.33 ± 0.17) X 10-3 [2] BE 
Liq. (3.30 ± 0.20) X 10-3 [210] CE 
12 (4.05 ± 0.23) X 10-3 [210] C E 
1 (4.27 ± 0.23) X 10-3 [133] O E 
1 (4.81 ± 0.49) X 10-3 [113] >IXE 


Liq. (1.10 ±0.10) X 10-3 [42] BE 
Liq. (0.99 ± 0.05) x 10-^ [167] C E 

K+K- 

12 (9.07 ± 0.59) X 10-* [210] C E 
0.005 (4.6 ± 0.3) X 10-* [134] O E 


Liq. (1.01 ± 0.05) X 10-3 [2] BE 
Liq. (1.10 ± 0.07) X 10-3 [210] CE 
1 (6.92 ± 0.41) X 10-* [113] >1E 
1 (2.87±0.51) X 10-*[113]>IXE 



Table 1 6: Frequencies for annihilation into two neutral strange mesons 



Chan. Dens. Ann. Freq. 




Ref. 


Type 


Dens. Ann. Freq. 


Ref. Type 




Liq. (6.10 ±0.90) x 10" 


-4 


[42J 


B Ekc 


Liq. (7.60 ± 0.40) x 10" 


-4 [2J BE 




Liq. (9.00 ±0.60) x 10" 


-4 


[176] 




Liq. (7.80 ±0.76) x 10" 


-U139]OEkc 


KsKi 


Liq. (8.64 ± 1.02) x 10" 


-4 


[210] 




12 (7.04 ± 0.74) X 10" 


-4 [210]CEK„,m6 




12 (4.89 ±0.56) X 10" 


-4 


[210] 


C Exn.an 


1 (3.60 ± 0.60) X 10" 






1 (3.50 ± 0.54) X 10- 


-4 


[139] 




0.005 (1.00 ±0.32) X 10- 


-4[139]OEkc 


KsK« 


Liq. (4.0 ± 3.0) X 10- 
1 (3.00 ± 1.00) X 10- 


6 

-5 


[42,114JBEkc 
[114] AEkc 


Liq. (7.0 ± 3.5) X 10- 
1 (3.70 ± 1.40) X 10" 


5 [215]BEkc 
-5[114]^XEkc 



is good agreement of the annihilation frequencies measured by the Crystal Barrel experiment for 
several reactions with different numbers of photons in the final state. In [167] various final states 
with 4 (or 5) and with 8 (or 9) photons in the final state were measured. 

The value of AF(7r'^7r'^, liq) has been used by the Crystal Barrel collaboration to normalise other 
annihilation frequencies with only neutral particles in the final state [167,215]. Using the fre- 
quency for AF(7r°7r°, liq) [167] to normahse the anniliilation frequency for the uiiv final state gives 
AF(a;w,liq) = (3.32 ± 0.34)% [167] and more recently AF(wa;,liq) = (2.95 ± 0.15)% [215]. 
The absolute pp^ ujiuj2 annihilation frequency with uji — > tt+tt-tt", 0J2 has also been 

measured by them [197]. AF(a;a;,liq) = (3.23 ± 0.25)% was obtained, in agreement with the 
normaUsed results, giving confidence in the measured AF(7r°7r°, liq). 

Other support for the efficiency determination for the Crystal Barrel detector can be made by 
a comparison of annihilation frequencies both for all neutral and for final states involving charged 
particles. These include measurements by the Obehx collaboration [139] for the reaction pp — > 
KgKi; Kg TT+TT- at three target densities. The Obelix value AF(KsKi,liq) = (7.8 ± 0.7 ± 
0.3) X IQ-^ is in good agreement with those measured [176, 210] by the Crystal Barrel experiment, 
AF(KsKi,liq) = (9.0±0.6)xlO-*andAF(KsKi,liq) = (8.6± 1.0) x 10-*, for the same reaction 
but with Kg ^n°n°. 
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Table 1 7: Annihilation frequencies for pp -+ 7r°7r' 



Dens. 


Ann. Freq. 


Ref. 


Type 


Ann. Freq. 


Ref. 


Type 


Liq. 


(4.8 ± 1.0) X 10"^ 


[89] 


TE 


(1.4 ±0.3) X 10"* 


[292] 


TR 


Liq. 


(6.0 ±4.0) X 10"^ 


[263] 


TR 


(2.06 ±0.14) X 10"^ 


[291] 


jCR 


Liq. 


(2.5 ±0.3) X 10-^ 


[106] 


K.R 


(6.93 ± 0.43) X 10-^ 


[166, 167] 


CE 


Liq. 


(6.14 ±0.40) X 10"* 


[210] 


CE 


(2.8 ± 0.4) X 10"^ 


[158] 


OE 


12 


(1.54 ±0.09) X 10"^ 


[210] 




(1.62 ±0.15) X 10-3 


[210] 


CE 

^ ^an 


1 


(1.27 ±0.21) X 10-3 


[131] 


OE 









Further evidence as to the consistency of the Crystal Barrel results is given by recent measure- 
ments [215] of annihilation frequencies for final states containing up to five 7r° and r] mesons, using a 
liquid target. Events with between 4 and 1 1 photons were processed and frequencies obtained by nor- 
malisation to the Crystal Barrel result [ 1 67] for the tt^tt'' channel. Including annihilation frequencies 
for the production of neutral kaons from other experiments, the identified channels from annihilation 
in Uquid hydrogen add up to a total frequency for all neutral annihilations of (3.56 ± 0.28)% per 
annihilation compared to a value of (3.5 ± 0.3)% for all neutral annihilations measured inclusively 
in the same experiment without using the tt^tt^ normalisation. Both values are in good agreement 
with the values (3.2 ± 0.5)% and (T.lt^ ;g)% obtained [2] with bubble chambers at BNL and CERN 
respectively. Using the Obelix value for the tt^tt" frequency to normalise the Crystal Barrel results 
would give a total all-neutral armihilation frequency of(1.6±0.25)%in marked disagreement with 
the other measurements. 

In the above discussion we have emphasised the internal and external consistency of the Crys- 
tal Barrel results relating to their annihilation frequencies for the 7r°7r° channel. Unfortunately no 
similar checks are available for the Obelix results, although the authors do point out [158] that mea- 
surements of pp tt+tt^tt", both with and without the tt" being detected are in good agreement 
showing that their estimation of the Obelix efficiency for detection of a single 7r° is reliable. We 
conclude however, there is a large body of evidence that the annihilation frequencies measured with 
the Crystal Barrel detector for a variety of different final states are consistent. 

The Obelix collaboration has also measured [13 1] the 7r°7r° annihilation frequency with a gaseous 
H2 target at STR This measurement has again been discussed recently by both the Obehx [158] and 
Crystal Barrel [214] collaborations. The latter suggest that, on the basis of a fit [210, Fig. 7(a)] with 
Eq. (6.4) to Crystal Barrel and Asterix data, the measured AF(7r°7r°, psxp) by Obelix [131] could 
be too low by a factor of about 1.4. 

6.4.4 Other two neutral pseudoscalar channels 

Most of the measurements for two-body annihilation frequencies into other neutral pseudoscalar 
mesons, pp tt^t/, tt"?]', 7777 and 7777', have been made by the Crystal Barrel collaboration [167, 
210,215]. As for the case of the AF(7r°7r°, liq) results, measurements for AF(7r°77, liq) made by the 
Crystal Barrel and Obelix [158] experiments again differ by a factor about 2. However it has been 
pointed out that the ratio AF{n'^ri, liq) / AF(7r"7r" , liq) obtained by the two experiments, 0.32 ± 0.07 
(Obehx [158]) and 0.303 ± 0.010 (Crystal Barrel [167]) are in good agreement. 

Whilst the recent Obehx [158] value for AF(7r°?7, liq) is supported by the measurement of Adiels 
et al. [293], it is not in agreement with other early results [106, 263] which show particularly large 
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Table 1 8: Frequencies for annihilation into two neutral pseudoscalar mesons 



Chan. Dens. Ann. Freq. Ref. Type 


Dens. Ann. Freq. Ref. Type 


Liq. (82.0 ± 10.0) x 10-^[263]TR 
Liq. (1.33 ±0.27) x 10~^[293]£R 
7r°r? Liq. (2.50 ± 0.30) x 10"* [210]C E 
Liq. (2.09 ±0.10) x 10-^[215]CEN 
12 (4.57 ± 0.30) X 10-^[210]C Ean 


Liq. (4.6 ± 1.3) x 10"^ [106]A: R 
Liq. (2.12 ±0.12) x 10-*[167]CEN 
Liq. (0.90 ±0.22) x 10-4[158]C»E 
12 (5.63 ± 0.43) X 10-4[210]C E^b 
12 (4.78 ± 0.21) X 10-4 [215]C E N 


g ^ Liq. (5.0 ± 1.9) X 10"^ [106]/CR 
^ ^ Liq. (0.98 ± 0.24) x 10"^ [215]C E N 


Liq. (1.23 ±0.13) x 10-4[167]CEN 
12 (2.03 ± 0.13) X 10-4[215]C E N 


Liq. (1.60 ±0.80) x 10-^[107]/CR 
rjr] Liq. (1.64 ± 0.10) x lO'^ [167]C E N 
12 (3.17±0.14) X 10-4[215]CEN 


Liq. (0.81 ±0.31) x 10-4[293]£R 
Liq. (1.53 ±0.08) x 10-4[215]CEN 


Liq. (2.16 ±0.25) x 10-4[167]CEN 

vv' 

12 (3.81 ± 0.28) X 10-4 [215]C E N 


Liq. (2.49 ±0.33) x 10-4[215]CEN 



fluctuations in value. This may be compared to the situation for the pp — > tt^tt" frequency discussed 
in Sec. 6.4.3 above which has been used by the ObeHx coflaboration [158] as support for their low 
value of BR(7r°7r liq.). We note that all these early results, which were obtained from measurements 
of the inclusive recoil energy spectrum, are expected to be less reUable than more recent results in 
which the events were fully reconstructed. 

6.4.5 Frequencies for annihilation into a pseudoscalar and a vector meson 

For the ojn'^.ujri and ujrj' channels, the majority of the annihilation frequencies (Table 19) were 
measured by the Crystal Barrel collaboration. The results for hquid H2 targets are only in moderate 
agreement with some early measurements. These latter results were obtained from measurements 
of the inclusive recoil energy spectrum. As the wtt^, utr] and ujt]' channels are produced solely from 
^Si and ^Pi initial pp states, the recent Crystal Barrel measurements [215] at a target density of 
12psTP, allow the separate contributions from these two initial states to be determined. 

Measurements for the (/)7r" and (prj channels (Table 19) have mostly been made by the Obelix and 
Asterix collaborations, detecting the (p meson through its decay (p ^ K+K". There are significant 
discrepancies between the results from these two experiments for the (pij channel measured with a 
gas target at a density psTP, although those for the (67r" channel are in reasonable agreement. Results 
from the Crystal Barrel experiment using a liquid H2 target, detecting the decays (j) K+K" or 
(p — > KiKs — > Kitt'^tt" are generally in good agreement with the other measurements. The Asterix 
collaboration [124] has also determined the annihilation frequencies for these channels with pp 
atomic L X-rays in coincidence, so selecting annihilation from atomic P-states. 

For channels involving the p meson (p° or p^), annihilation frequencies are listed in Tables 19 
and 20. The frequencies in Table 19 were generally obtained from measurements of the recoil energy 
spectrum or from an analysis of invariant mass distributions. Those in Table 20 were obtained by 
a partial wave analysis of the data; in some cases [62, 75] for data obtained with a liquid hydrogen 
target and assuming that S-state annihilation dominates. Whilst many of the measurements in Tables 
19 and 20 are rather old and made with the bubble-chamber technique, the results are generally 
in moderate agreement with each other and the recent Crystal Barrel values. Again an exception 
is the rather high value for the p^r] chaimel obtained at KEK [107] by measuring the inclusive 
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Table 19: Frequencies for annihilation into a pseudoscalar and a vector meson 



Chan. Dens. Ann. Freq. Ref. Type 


Dens. Ann. Freq. Ref. Type 


Liq. (23.8 ± 6.5) x IQ-^ [263] T R 
WTT^ Liq. (5.73 ± 0.47) x IQ-^ [167] C E N 
Liq. (6.00 ± 0.30) x IQ-^ [215] C E N 


Liq. (5.2 ± 0.5) x 10-3 [106] JC R 
Liq. (6.16±0.44) x 10-3[4,294]C E 
12 (4.60 ± 0.24) X 10-3 [215] CEN 


Liq. (0.46 ±0.14) X 10-2 [107] /CR 
LJT] Liq. (1.51 ± 0.12) X 10-2 [167] CEN 
Liq. (1.53 ±0.06) X 10-2 [215] CEN 


Liq. (1.04 ±0.10) x 10-2 [293] CE 
Liq. (1.63 ±0.12) x 10-2[4,294]C E 
12 (9.31 ±0.42) X 10-3 [215] CEN 


Liq. (7.8 ± 0.8) X 10-3 [167] CEN 
12 (7.32 ± 0.56) X 10-3 [215] CEN 


Liq. (8.46 ±0.97) x 10-3 [215] CEN 


Liq. (3.3 ± 1.5) X 10-4 [106] ICR 
Liq. (4.88 ± 0.32) x 10-"^ 1 148] O E^c 
(l)Tr° 1 (1.90 ± 0.50) X 10-^ [124] AE^c 
1 (2.47 ± 0.21) X 10-^ [148] O E^c 
0.005 (0.92 ± 0.10) X 10-^ [148] O E^c 


Liq. (6.50 ± 0.60) X 10-'' [176| C E^„,0c 
Liq. (5.20 ±0.60) x lO-'' [204| C E^c 
1 (2.46 ± 0.24) X 10-* [137] OE^c 
1 (0.30 ± 0.30) X 10-* [124] ^IXE^c 


Liq. (7.80 ±2.10) x 10-^ [170, 176]C E^n" 

1 (13.3 ± 1.5) X 10-^ [147] CE^c 
1 (4.10 ± 1.60) X 10-5 [124] AXE^c 


Liq. (7.10 ±0.70) x 10-^ [147] OE^^ 
0.005 (16.6 ± 2.0) X 10-5 [147] O E^c 


Liq. (1.40 ±0.20) X 10-2 [48] BE 
^ Liq. (1.6 ± 0.1) X 10-2 [106] ICR 


Liq. (1.72 ±0.27) X 10-2 [2] BE 
Liq. (1.58 ±0.12) x IO-2 [205] CE 


p±7rT Liq. (2.90 ±0.40) X 10-2 [48] BE 


Liq. (3.44 ±0.54) X 10-2 [2] BE 


Q Liq. (2.20 ± 1.70) X 10-3 [48] BE 
Liq. (5.30 ±2.00) X 10-3 [293] CR 


Liq. (9.60 ± 1.60) x 10-3 [107] ICR 


p°r]' Liq. (1.46 ±0.42) X 10-3 [295] CE 





"Updated by Amsler [4] 



energy spectrum. The weighted mean values are AF(/9"7r", liq.) = (1.57 ± 0.07) x 10 2 with 
X^/N = 1.1/3, withNthe number of degrees of freedom and AF(p077,liq.) = (3.93±0.28) x 10-3 
with x^/N = 5.7/5, where in the latter case the KEK measurement [107] has been omitted. 

The p^tt" and p'^ i] channels are only allowed from the 3Si and ^Pi initial states of the pp system 
whilst the p^7r^ channel is allowed from the ^Sq, 3Si, ^Pi, 3Pi and 3p2 states as shown in Table 5. 
The Asterix collaboration have measured [121, 122, 125] frequencies for pp annihilation at rest into 
7r+7r-7r° and into tt+tt- ;/ in hydrogen gas (psTp), both with and without a trigger on atomic L X- 
rays. In this way they obtain two data samples with different fractions of S- and P-state annihilation. 
It is then possible from a Dalitz plot analysis to obtain separate frequencies for the channels pp 
p^TT^ [297], p^TT^ [121, 122] and for pp p^rj [125] from initial S- and P-states. These values 
are hsted in Table 20. The Obelix collaboration [160] has performed a coupled-channel partial- 
wave analysis of data for ■jT~^n~TT^, K='=Ki7r^ and K+R-tt'^ at three values of the target density. 
Aimihilation frequencies for S- and P-states for pw chaimels are also hsted in Table 20 where the 
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Table 20: Annihilation frequencies for annihilation into pn, pi], and K*K from partial wave analyses 







Ann T^vf^rt 


T?pf 


Type 


Ann Frf^n 




Type 




P TT^ 


< U.i4 X iU 


Lozj 


ID tL 








J^o 


p TT^ 


no -U n Oyi^ ^ i n~2 
^u.uy=tu.u4j X iu 


n on 


>l o Ji 


(0.19 ±0.03) X 10-2 


[160] 


OSE 




p TT 


(-1 CO _L n OCA w 1 n-2 


LOZJ 














f1 fiQ + f) 221 X 1 f)~^ 

I J.. Lit? —1— yjt^^i /n ±l/ 






(1.58 ±0.09) X 10-2 


[160] 


OSE 


^ 1 




ff) Af) _|_ n nq'i y i n-2 




>4 P F 


(0.43 ±0.04) X 10-2 


[160] 


OPE 


■1^12 






1 1 991 


yd P F 


(0.77 ±0.17) X 10-2 


[160] 


OPE 


'Si 




l^o.zy HZ u.yuj a iu 


ri 9S1 




(6.40 ± 1.40) X 10-3 


[63] 


BE 




p^h] 


^O.UU ZIZ l.^U ) A lU 


L/JJ 




(3.87 ±0.29) X 10-3 


[190] 


CE 


1 


p'^rj 


1 U.c'4: ZIZ U.<JO ) A lU 


ri 9S1 


^ P F 










H '1 


(1.81 ± 0.44) X 10"'^ 


[63, 125] 


ABE 






K*K 


(24.6 ±2.4) X 10-'' 


[80] 


BE 


(16.6 ±2.5) X 10-" 


[196] 


CE 




K*K 


(18.9 ±2.0) X 10-"* 


[284] 


OSE 


(19.1 ±3.3) X 10-" 


[296] 


CE 




K*K 


(0.8 ±0.4) X 10-'' 


[80] 


BE 


(4.5 ± 1.2) X 10-" 


[196] 


CE 




K*K 


(2.8 ±0.4) X 10-'' 


[284] 


OSE 


(0.4 ±3.1) X 10-" 


[296] 


CE 




K*K 


(2.0 ±0.4) X 10--* 


[80] 


BE 


(0.9 ±0.3) X 10-" 


[196] 


CE 




K*K 


(2.3 ±0.5) X 10-^ 


[284] 


OSE 


(1.1 ±1.5) X 10-" 


[296] 


CE 




K*K 


(8.6 ± 1.7) X 10-^ 


[80] 


BE 


(6.2 ± 0.9) X 10-" 


[196] 


CE 




K*K 


(9.0 ±0.9) X 10-^ 


[284] 


OSE 


(11.9 ± 1.5) X 10-" 


[296] 


CE 



S-wave frequencies are derived from the analysis in liquid H2, the P-wave ones from the analysis at 
p = 0.005psTP- Recently the Obelix collaboration published a further paper on these experimental 
results [284], which gives annihilation fractions obtained from an average of good fits to the data. 
These latest results, which in some cases differ slightly from those presented in Table 20, give a 
more realistic estimate of the systematic errors. 

There are almost no selection rules due to quantum number conservation for pN annihilation into 
K*K + c.c, only annihilation of the ^Pq state into 3 pseudoscalar mesons is forbidden. Annihilation 
frequencies for NN ^ K*K + c.c. were determined from bubble-chamber data on pp KgK^'^Tr^ 
(2000 events) and KsKs7r° (364 events) [59]. Later [80], data on pn KsK-7r°, K^KsTr-, 
KgKiTT-, and K+K-tt" were obtained. The latter final states have low statistics only (655 events in 
four Dalitz plots), but the data are related by isospin invariance and provide valuable constraints for 
the partial-wave analysis. Since the analysis [80] includes the data from [59], we use only [80]. In 
all cases, the K„ was identified through its tt+tt- decay mode and its secondary vertex. 

The Crystal Barrel collaboration has analysed 11373 events due to pp KiK='=7r^ with an 
undetected Ki [196]. Results on K*K -|- c.c. are derived from a partial wave analysis. The bubble 
chamber and Crystal Barrel analyses assume only S-state capture. 
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The Crystal Barrel Collaboration has also taken data on the reaction pp KsK^'^tt^ by trig- 
gering on secondary Kg tt+tt" decays. About 50 k events were recorded, nearly fourfold the 
statistics obtained in [196]. Using the larger data sample, more amplitudes could be included in the 
fit. In particular, amplitudes to describe K2(1430)K production were included [296]. Most results 
are not too different from those reported in [196]. We therefore include in Tables 20 and 22 the 
results on K*K and K^(1430)K production. 

As mentioned above, the Obelix collaboration studied the K^KgTr^ and K+K~7r'' final states at 
three different target densities [160]. Revised annihilation frequencies are again given in their recent 
paper [284]. The annihilation frequencies for the K*K channel from this latter paper, averaged over 
the two final states, are again given in Table 20. 

6.4.6 Frequencies for annihilation into two vector mesons 

As shown in Table 21, only a htnited range of frequency measurements are available for annihilation 
into two vector mesons. In Sec. 6.4.3 we commented on the good agreement between the three 
measurements [167, 197, 215] of AF(cja), liq.) obtained by the Crystal Barrel experiment. 

The Asterix collaboration measured [124] the annihilation frequency for pp — » (j)u) with an 
atomic L X-ray trigger as well as the frequency in a gas target (psTp). Whilst the annihilation 
frequency has also been measured [72] for a liquid target, these three results together are insufficient 
to determine partial branching ratios, since the channel occurs from four initial pp states, "^Si, ^Pq, 
^Pi and ^P2. Note that the possibility of ^Pi — > <pco was inadvertently omitted in [259, Table 4c]. 

Table 21 : Frequencies for annihilation into two vector mesons 



Chan. Dens. Ann. Freq. Ref. Type 


Dens. Ann. Freq. Ref. Type 


Liq. (3.32 ±0.34) x lO^^ (-igy-j^ e^^^^„ 

LOU) 

Liq. (2.95 ±0.15) x 10"^ [215]C N Ea,n,a;n 


Liq. (3.23 ±0.25) x 10-^ [197]C E„c,t^„ 
12 (3.58 ± 0.19) X 10-2 [215]C N E^n,ojn 


Liq. (6.30 ± 2.30) X 10-^ [72]BE^c 
1 (4.20 ± 1.40) X lO-^[124]^XE0c,a;c 


1 (3.00 ± 1.10) X 10-^[124]>lE^c,c^c 


Liq. (0.70 ± 0.30) X 10-2 [48]BE 
1 (2.95 ± 0.72) X 10-2 [126]A S E 


Liq. (2.26 ±0.23) x 10-^ [67] BE 
1 (6.35 ± 1.14) X 10-2 [126]APE 


pV 1 (3.40 ± 0.80) X 10-^ [124] >lE^c 


1 (4.40 ± 1.20) X 10-*[124]>lE^c 


pV Liq- (3.8 ± 3.0) X 10-3 [48] BE 


Liq. (1.2 ± 1.2) X 10-3 [70]BE 


K*+K*- Liq. (1.5 ± 0.6) X 10-3 [47]BE 
K^or^o Liq. (3.0 ± 0.7) X 10-3 [47]BE 




State 


State 


K*+K'- -1.7x10-''' [192]CE 
K*''K*« ~ 1.5 X 10-3 [192]CE 


3Si -1.3x10-3 [192]CE 
3Si ~ 0.2x10-3 [192]CE 



The Asterix collaboration also measured [126] annihilation frequencies for pp annihilation at 
rest into five pions in hydrogen gas (psTp), both with and without a trigger on atomic L X-rays. In 
this way they obtain two data samples with different fractions of S- and P-state annihilation and can 
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derive separate frequencies for the channel pp — > p^ui from initial S- and P-states. These values are 
also hsted in Table 21. 

Values of AF(K*+K*-) and AF(K*0K*0) for ^Sq and ^Si states were obtained by the Crys- 
tal Barrel collaboration from a Dalitz plot analysis [192, Table 1] of pp annihilation at rest into 
KiK^7r^7r°. The total annihilation frequency for this latter final state was taken from Table 9. 

6.4.7 Two-body annihilation frequencies involving pseudo- vector or tensor mesons 

A useful compilation of annihilation frequencies obtained from Dalitz plot analyses for two-body 
pp annihilation at rest has been given by Amsler [4, Table Xlll]. This includes information for final 
states including either pseudo-vector (bi (1235)) or tensor (£2 (1270), a2 (1320)) mesons. Amsler [4] 
also discusses some of the difficulties associated with determining these annihilation frequencies. 
Note that all these states are broad with widths in the range F = 100 to F = 185 MeV. 

The Asterix collaboration measured [125] annihilation frequencies for pp annihilation at rest 
into 7r+7r~?7 and 7r+7r~?7' in hydrogen gas (psTp), both with and without a trigger on atomic L X- 
rays. In this way they obtain two data samples with different fractions of S- and P-state annihilation. 
They are then able from a Dahtz plot analysis to obtain separate frequencies for the charmel pp 
a2(1320)='=7r^ from initial ^Sq, ^Si and P-states, where the latter is averaged over all fine-structure 
states (^Pi, ■'Pi and ''P2). The same analysis gives separate S- and P-state annihilation frequencies 
for the channel pp £2(1270)77. 

Similarly, the 7r+7r~a; final state was measured [126], with ui — > 7r+7r~7r°, and the S- and P- 
state annihilation frequencies for the annihilation channel pp b|j'^(1235)7r^ were determined. 
Again the P-state fraction is averaged over all fine-structure states and the analysis does not take into 
account any effects due to the enhancement factors discussed earlier in Sec. 6.2. The 7r+7r~7r° final 
state has been measured [121, 122] and analysed in a similar way to give S- and P-state frequencies 
[297] for the annihilation channel pp £2(1270)77°. 

The Obehx collaboration [160] has performed a coupled partial- wave analysis of data for the 
channels tt^t:~tt^, K^Kitt^ and K+K^tt" at three values of the target density. An additional 
paper by the collaboration on this work was published later [284], which gives annihilation fractions 
obtained from an average of good partial-wave fits to the experimental results. These measurements 
are also listed in Table 22. 

The annihilation frequencies hsted in Table 22 include all charge states and decay modes; they 
were calculated from the final states which are also listed. This involves in some cases using the 
squares of the isospin Clebsch-Gordan coefficients to determine the total annihilation frequencies 
including all charge modes [4, Table XII]. The annihilation frequencies must also be corrected for 
all decay modes of the intermediate resonances. 

The weighted mean value for AF(£2(1270)7r", ^So) = (3.70 ± 0.47) x lO"'"' is used in the 
analysis of Sec. 6.5. Measured annihilation frequencies for a2(1320)7r show a very wide spread, 
with values obtained in ref. [70] appearing to be anomalously big. Excluding these latter results 
gives weighted mean values AF(a2(1320)7r, ^So) ^ (2.04±0.30) x lO"^ ^ith xV^ = 19.9/7 and 
AF(a2(1320)7r, ^Si) = (0.48 ± 0.09) x lO'^ with x^/N = 3.3/4. Simple (unweighted) averaging 
gives AF(a2(1320)7r,iSo) = (2.8± 1.2) x lO'^ and AF(a2(1320)7r, ^Si) = (0.59 ± 0.23) x 10"^ 
where the relatively large values for a are again an indication of the wide spread in experimental 
values. 

6.4.8 Other relevant annihilation frequency information 

In Table 23 we list a number of cases where the ratio of annihilation frequencies for two annihilation 
channels has been directly measured. We have not included values derived from separate frequency 
measurements where in some cases [167], the error is reduced since the systematic error is common 
to the two frequency measurements and so cancels out. The values of AF(KsKs, p)/AF(KsKi, p) 
at densities of 15 and 27 psTP [298] and of AF(K+K-, p)/AF(7r+7r-, p) at 0.002 psTP [134] are 
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particularly useful as they allow the fraction of P-state annihilation at these three target densities 
to be determined. In the absence of other annihilation frequency measurements at these densities, 
this determination relies on the availability of branching ratios for these channels which have been 
obtained from measurements of absolute annihilation frequencies at other target densities. 

The ratio AF(f2(1270)7r°, p)/AF{p°Tr°, p) has been estimated [297] by the CPLEAR collabo- 
ration at a target density of 16psTP and used by them together with the separate frequencies for 
annihilation from S- and P-states for the two-body channels f2(1270)7r" and p^Tr" to obtain the frac- 
tion of P-state annihilation fp(16 psTp) = 0.38 ± 0.07. This estimate is not based on a partial-wave 
analysis and does not take into account the effect of the enhancement factors. For these reactions the 
initial states are ^Sq, ^Si, ^Pi, and ^Si, ^Pi, respectively. With the "^Pq state not contributing to 
either reaction, the enhancement factors for the relevant states are close to 1 .0 and should have little 
effect on the derived fp (16 pstp)- The latter value is in good agreement with that obtained from the 
more detailed analysis of a wide range of data to be discussed in Sec. 6.5. See also Table 24 and 
Fig. 14. 

In Table 25 we list production rates for the channel pp — > r?(1440)7r+7r~ with ry(1440) — > 
K^KiTT^. The values obtained by the Obelix collaboration at three target densities [138] are in good 
agreement with those obtained previously in Uquid hydrogen using the Bubble Chamber technique 
[61] and by the Asterix collaboration [1 18] in H2 gas at STP. 

In the 77(1440)7r"'"7r~ production reaction, the dipion (tt+tt^) and the r](1440) resonance are 
produced with relative angular momentum L = [118, 135]. It is expected that the dipion system 
will occur mainly with relative orbital angular momentum £ = 0, since its invariant mass is less than 
500 MeV/c^. The assumption L = £ = then implies that 77(1440) production is only possible from 
the ^So initial state. As a consequence the 77(1440) production frequency is directly related to the 
population of the ^So state. 
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Table 22: Annihilation frequencies for annihilation into a tensor or a pseudoscalar meson from partial 
wave analyses (liquid H2) 



State 



Chan. Ann. Freq. 



Final State Ref. 



Type 








f2(1270)7r" 
f2(1270)7r0 
f2 (1270)71° 
f2(1270)7r° 
f2(1270)7r° 
f2(1270)7r° 



4.3 ± 1.2) X 10-3 
3.8 ± 1.0) X 10-3 
3.1 ± 1.1) X 10-3 
3.7 ±0.7) X 10-3 
2.7 ±0.2) X 10-3 
2.6 ±0.8) X 10-3 



TT TT IT 



K+K-ttO 



[621 
11211 
[4] 
[4] 
[284] 
[284] 



BE 
AE 
CE 
CE 
OE 
OE 



'So 
'So 



f2 (1270)7/ 
£2(1270)77 



0.15 ±0.15) X 10-3 
0.13 ±0.13) X 10-3 



TT+TT 77 
77+77-7/ 



[125] 
[75] 



AE 
BE 



^So (1525)77° 



9.38 ± 1.49) X 10" 



ttORiKi [185] BE 



'So 
'So 
'So 
'So 

•0 
•0 
'So 
'So 
'So 
'So 
'So 



'S, 



'S, 



a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 



6.3 ±0.9) X 10-2 
1.3 ±0.4) X 10-2 
4.8 ± 1.5) X 10-2 
2.2 ±0.4) X 10-2 
2.7 ±0.8) X 10-2 
3.93 ±0.70) X 10-2 
3.36 ±0.94) X 10-2 
1.55 ±0.31) X 10-2 
2.4410:^4) X 10-2 
1.2 ±0.2) X 10-2 
1.2 ±0.2) X 10^2 



27r+27r- 
K±Ks7r=F 

77+77-7/ 
K0K±7r=F 

77+77-7/ 
27r°7/ 

277°^ 

7r°KiKi 
K+Ki7r=F 
K+K-7r° 

K+K07r+ 



[70] 
[59] 
[75] 
[75,80] 
[125] 

[4] 

[4] 

[4] 

[4] 
[284] 
[284] 



BE 
BE 
BE 
BE 
AE 
CE 
CE 
CE 
CE 
OE 
OE 



a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 
a2(1320)7r 



1.7 ±0.4) X 10^2 
0.45 ± 0.18) X 10^2 
0.69 ±0.34) X 10-2 
0.31 ± 0.16) X 10^2 
0.90 ±0.34) X 10^2 
0.58 ±0.20) X 10-2 
;0.33±0.05) X 10-2 



27r+27r- 
K±K,7r+ 

TT+TT-?/ 

K"K+7r+ 
Tr+TT-?/ 
K±Ki7r+ 
K+K07r+ 



[70] 
[59] 
[75] 
[75,80] 
[125] 

[4] 
[284] 



BE 
BE 
BE 
BE 
AE 
CE 
OE 



3.3Si 

l'3Si 



K^(1430)K 
K^(1430)K 
K^(1430)K 
K^(1430)K 



5.1 ±1.8) X 10-4 
0.5 ± 1.8) X 10-4 
0.2 ±0.8) X 10-4 
2.6 ± 0.8) X 10-4 



K+K"7r+ 
K+K07r+ 
K+ROtt^f 
K+ROtt^f 



[296] 
[296] 
[296] 
[296] 



CE 
CE 
CE 
CE 
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Table 23: Directly measured ratios of annihilation frequencies 



Channels 


Dens. 


Ratio 


Ref. 


Type 


K,K,/K,Ki 


15 


0.041 ± 0.009 


[298] 


^ Ekc 


K,K,/K,Ki 


27 


0.037 ±0.002 


[298] 


^ Ekc 


K+K-/7r+7r- 


15 


0.205 ±0.016 


[299] 


£E 


K+K-/TT+TT- 


1 


0.163 ±0.011 


[129] 


OE 


K+K-/7r+7r- 


0.002 


0.102 ± 0.015 


[134] 


OE 


pto 1 pTr" 


16 


2.67 ±0.45 


[297] 


£E 


f2 (1270)^7/9"^° 


16 


0.83 ±0.06 


[297] 


/:e 


f2(1270)77/f2(1270)7r° 


16 


0.37 ±0.16 


[297] 


£E 



Table 24: Fraction of P-state annihilation as a function of H2 target density p/psxp 



Dens. 0.002 0.005 1.0 12.0 15.0 27.0 Liq. 

Pfrac. 0.89 ±0.07 0.87 ±0.02 0.64 ±0.03 0.50 ±0.03 0.48 ±0.05 0.46 ±0.05 0.125 ±0.02 



Table 25: Production rate of the final state r7(1440)7r+7r", with 77(1440) -> K='=Ki7r''' 



Dens. 


Ann. Freq. 


Ref. 


Type 


Dens. 


Ann. Freq. 


Ref. 


Type 


Liq. 


(7.1 ± 0.7) X 10-^ 


[2,61] 


BE 


Liq. 


(6.0 ±0.5) X 10-^ 


[138] 


OE 


1 


(3.0 ± 0.9) X 10-^ 


[118] 


^E 


1 


(2.9 ±0.4) X 10"* 


[138] 


OE 


0.005 


(1.0 ±0.2) X 10"^ 


[138] 


OE 
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6.4.9 Compilation of two-meson annihilation frequencies 

In view of the discrepancies mentioned above, a compilation of preferred measurements of anni- 
hilation frequencies could be misleading since it requires rejection of some data points which we 
considered to be less reliable. Obviously, personal preferences could enter here. On the other hand, 
the reader may expect a 'final' list of annihilation frequencies which can be used to compare with 
model calculations, and also we shall need a table of frequencies for further analyses. With this 
warning we present in Tables 26, 27 and 28 our choice of annihilation frequencies. Aimihilation 
frequencies listed in Table 28 are weighted means of values taken from Tables 19, 20 and 22. 



Table 26: Two-meson annihilation frequencies at density p used in the analysis. 



Channel 


Units 


P — PSTP 


Ref. 


p = 12/9STP 


Ref. 


Liquid 


Ref. 




10- 


-3 


4.30 ±0.15 


[113] 


4.05 ±0.23 


[210] 


3.14±0.12t 


[167,210] 


tt"tt" 


10- 


-3 


1.27±0.21 


[131] 


1.56±0.08t [210] 


0.651 ±0.029t [166,167,210] 


K+K- 


10- 


-4 


6.92 ±0.41 


[113] 


9.07 ±0.59 


[210] 


10.3±0.4t 


[167,210] 




10- 


-5 


3.0 ± 1.0 


[114] 






0.70 ± 0.35 


[215] 


KsKi 


10- 


-4 


3.54±0.40t 


[114,139] 


4.89 ±0.56 


[210] 


8.56 ± 0.43t 


[139,176,210] 


(7T+7T~)x 


10- 


-3 


4.81 ± 0.49 


[113] 










(K+K-)x 


10- 


-4 


2.87 ±0.51 


[113] 










(KsK,)x 


10- 


-5 


3.70 ± 1.40 


[114] 










r?(1440)7r+7T- 


10- 


-4 


2.9 ±0.4 


[138] 






6.0 ±0.5 


[138] 


(jm 


10- 


-4 


2.41±0.15t 


[124, 137, 148] 




6.50 ±0.60 


[176] 




10- 


-5 


3.0 ±3.0 


[124] 












10- 


-5 


13.3 ± 1.5 


[147] 






7.10 ±0.70 


[147] 


TTT] 


10- 


-4 






4.78 ± 0.2 


[215] 


2.09 ±0.10 


[215] 


m 


10- 


-4 






3.17 ±0.14 


[215] 


1.53 ±0.08 


[215] 


ITT]' 


10- 


-4 






2.03 ±0.13 


[215] 


0.98 ± 0.24 


[215] 


m' 


10- 


-4 






3.81 ±0.28 


[215] 


2.49 ± 0.33 


[215] 


U)Tr 


10- 


-3 






4.60 ±0.24 


[215] 


6.00 ±0.30 


[215] 


IVT] 


10- 


-3 






9.31 ±0.42 


[215] 


15.3 ±0.6 


[215] 


UlT]' 


10- 


-3 






7.32 ± 0.56 


[215] 


8.46 ± 0.97 


1215] 



t Weighted mean value 
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Table 27: Two-meson annihilation frequencies used in tlie analysis (cont). 

PSTP- 



Annihilation dynamics 
The density is in units of 



Channel Units p = 0.002 p = 0.005 p = 15 p = 27 Refs. 

TT+TT- 10^3 4.26 ±0.11 [134] 

K+K^ 10-"* 4.6 ±0.3 [134] 

K+K-/7r+7r- 0.102 ±0.015 0.205 ± 0.016 [134,299] 

K3K1 10-^ 1.00 ±0.32 [139] 

KsKs/KsKi 0.041 ±0.009 0.037 ±0.002 [298] 

r/(1440)7r+7r- 10"^ 1.0 ±0.2 [138] 

(jyn 10-^ 1.66 ±0.20 [147] 



Table 28: Annihilation frequencies for prr, p°r), p°r)', K*K, 7r°f2(1270), 7?f2(1270) and 7ra2(1320), used 
in the analysis. 



State Chan. Ann. Freq. 


State Chan. Ann. Freq. 


^So p=^7r=F (0.15 ± 0.05) X 10-2 
^Si p°n° (1.59 ± 0.08) X 10-2 
^Si p°?7 (3.95 ± 0.32) x IQ-^ 
^Si pV (1.63 ± 0.30) X 10-3 
iPi p°7r° (0.43 ± 0.04) X 10-2 
^So 7r°f2(1270) (3.70 ± 0.47) X 10-3 
^So 7ra2(1320) (2.04 ± 0.30) x 10-^ 


3'3Si K*K (20.8 ± 4.0) X IQ-^ 
I'^Si K*K (1.2 ± 1.1) X 10-^ 
3'iSo K*K (1.3 ± 0.5) X 10-4 
I'^So K*K (6.7± 1.0) X 10-4 
iPi p% (1.84 ± 0.44) X 10-3 
^So r?f2(1270) (0.14 ± 0.10) X 10-3 
3Si 7ra2(1320) (0.48 ± 0.09) x 10-^ 



6.5 Analysis of two-body annihilation frequencies for pp atoms 

The annihilation frequencies summarised in Tables 26 and 27 were analysed using Eq. (6.4). The 
enhancement factors Ei^^^^Lj, p) were fixed at the values of Table 12, calculated from cascade 
calculations. Calculated annihilation frequencies were then fitted using the least squares method to 
give a best fit to the experimental measurements by varying the branching ratios BR(c/i, "^^^"^Lj) 
and fraction of P-state annihilation fp(p). The search was constrained so that BR(c/i, '^^^^Lj) > 0. 
For those frequency measurements, AF(c/i)x, made in coincidence with atomic L X-rays, Eq. (6.9) 
was used. 

Two-body final states are only allowed from certain initial states of the pN system as shown in 
Table 5. For most channels this gives a considerable reduction in the number of free parameters 
BR(c/i, in Eqs. (6.4) and (6.9) when fitting the annihilation frequency data. A further 

simplification is obtained for the tttt system by using charge symmetry when 



BR(7r°7r°,3po) = -BR(7r+7r-, 3Po) 



(6.10) 
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and 

BR(7r°7r°,3p2) = IbRCtt+tt-, ^P^) (6.11) 

Explicit equations for the frequencies for tt+tt^, tt^tt", K+K~, KgKg and KgKi two-body annihila- 
tion channels are given in the paper by Batty [282, Eqs. (3.13) to (3.21)]. 

For this review a new analysis has been made, following that made by Batty [282] but using an 
extended set of annihilation frequencies. As well as the measurements used in [282, Table 3], the re- 
cent two-body annihilation frequencies reported by the Crystal Barrel collaboration [210], were used 
together with a number of other recent measurements. The Obelix collaboration [139] has measured 
frequencies for the reaction pp KgKi; Ks tt+tt" at three densities, 0.005 psTP, Pstp and 
liquid, whilst the CPLEAR collaboration have measured [298] the ratio AF(KsKs, p) /AF(KsKi, p) 
at densities of 15 Pstp and 27 psTP- This latter measurement enables the P-state fraction fp(p) at 
27 psTP to be determined for the first time. The production of 7](1440) at three target densities in the 
reaction pp r;(1440)7r+7r~ has been determined [138] by the Obelix collaboration. As discussed 
in Sec. 6.4.8 this production frequency is directly related to the population of the ^Sq initial state. 

Together this gives a total of 36 frequencies covering the tt+tt^, tt'^tt'^, K+K^, KsKg and KgKi 
annihilation channels and ?7(1440) production, at a total of 7 densities. In those cases where several 
measurements are available for a particular annihilation channel and target density, the weighted 
mean value was used. This reduced to 29 the number of annihilation frequencies to be fitted. Seven- 
teen parameters were varied to obtain a least squares best fit; 10 branching ratios BR(c/i, 
and 7 values of the fraction of P-state annihilation fp (p) over the range of target densities from 
0.002psTP to liquid H2. 

A least squares fit to this data, which omits the Obelix value [158] of AF(7r°7r°, liq.) discussed 
in Sec. 6.4.3, gave a best fit with a per degree of freedom x^/^ = 21.7/11. Rather poorly 
fitted were the value of AF(KsKi, 12 psTp) = (7.04 ± 0.74) x 10"^ measured by the Crystal Barrel 
experiment [210] and the value of AF(7r°7r°, psTp) measured by the Obelix experiment [131] which 
also gave difficulties in the previous analysis [282]. Values of for these two measurements were 
6.8 and 4.8 respectively. Omitting these two annihilation frequencies gave a very good fit to the 
data with jN = 8.72/9. The values obtained for the branching ratios were very similar to those 
obtained in the earlier analysis [282, Table 6]. In particular BR(K°K°, ^Pg) = (0.0 ± 0.06) x lO'^ 
and BR(K"K0,3P2) = (0.20 ± 0.04) x lO'^. The relatively large value for the '^Pa branching 
ratio is in contradiction to the results for other channels where the ■^P2 branching ratio is consistent 
with zero. The analysis was therefore repeated with the recent value [215] for AF(KsKs, liq.) = 
(7.0 ± 3.5) X 10"^ replacing the one [42, 1 14] used in the earlier analysis. An equally good fit was 
obtained with x^/iV = 8.38/9 and a value of BR(K'^K'^, ■'P2) consistent with zero. Values obtained 
for the fraction of P-state annihilation are given in Table 24 and plotted in Fig. 14. Values for the 
branching ratios obtained from this best fit to the data are given in the first four fines of Table 29. 
Note that _ 

BR(K°kO, 3Po) = 2 BR(K,K„ ^Pq) = Dw, (6.12) 

where £>io is the parameter used by Batty [282, Eq. (3.19)]. Similar expressions apply for the ^P2 
state. 

For these least square fits, enhancement factors calculated with strong interaction widths [259] 
from the DRl potential were used. Very similar results, both for the fraction of P-state annihilation 
and for the branching ratios, were obtained with the DR2 and KW enhancement factors. A fit to the 
same data, but with the enhancement factors set equal to 1 (i.e. no enhancement) gave a significantly 
worse fit with x'^/N = 26.3/9. 

As a second stage in the analysis, the range of annihilation frequencies was extended to include 
the (j)TT and (f>ri channels. For the 4>Tr channel the measurement for a liquid H2 target by the Crystal 
Barrel collaboration [176] was used, for psTP gas the average of the Obelix [137, 148] and Asterix 
[124] measurements; the Asterix [124] measurement with pp atomic L X-rays in coincidence was 
also included. For the (prj channel, the Obelix measurements [147] in hquid, psTP and 0.005psTP 
gas were used. A best fit, in which the fraction of P-state annihilation and frequencies for the 
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Table 29: Two-body Branching Ratios from S and P wave atomic states, in units of 10 ^. The last 
column ('^Po*) corresponds to a fit where it is assumed that the ■^P2 state does not contribute. 



Channel 






'Pi 


'Po 


3P1 


'P2 


3p^* 




X 


X 


X 


46.0 ±8.9 


X 


1.7± 1.8 


54.3 ± 1.5 




X 


2.79 ±0.18 


X 


92± 18 


X 


3.4 ±3.6 


108.6 ±3.0 


nri 
1 


X 


X 


X 


14.3 ±2.9 


X 


0.7 ±0.8 


17.5 ±0.9 


TTT]' 


X 


X 


X 


7.2 ±0.4 


X 


0.0 ±0.05 


7.5 ±0.5 


VV 


X 


X 


X 


11.4 ±0.9 


X 


0.0 ±0.01 


12.1 ±0.6 


rji]' 


X 


X 


X 


14.2 ±0.9 


X 


0.0 ±0.01 


14.7± 1.1 


K+K~ 


X 


1.47 ± 0.06 


X 


4.25 ± 0.33 


X 


0.0 ±0.05 


4.20 ±0.30 




X 


1.31 ± 0.07 


X 


0.41 ±0.50 


X 


0.11 ±0.14 


0.79 ± 0.15 


1 


X 


21.2± 1.1 


17.2 ± 1.6 


X 


X 


X 


X 




6.0 ±2.0 


42.4 ± 2.2 




X 






X 


LOTl 


X 


8.75 ±0.53 


12.1 ±2.8 


X 


X 


X 


X 




X 


0.87 ±0.09 


0.06 ±0.10 


X 


X 


X 


X 




X 


5.3 ±0.4 


7.4 ± 1.8 


X 


X 


X 


X 


p°r]' 


X 


2.2 ±0.4 




X 


X 


X 


X 


LOT} 


X 


23.1 ± 1.1 


7.0 ±5.0 


X 


X 


X 


X 


ujr}' 


X 


12.1 ± 1.6 


23.3 ±6.9 


X 


X 


X 


X 


dm 


X 


0.08 ±0.02 


0.73 ±0.07 


X 


X 


X 


X 


K*K(7 = 0) 


2.7 ±0.4 


0.16 ±0.15 




X 






X 


K*K(7 = 1) 


0.5 ±0.2 


2.8 ±0.5 




X 






X 


7r°a2(1320)0 


27.3 ±4.0 


X 




X 






X 


7r±a2(1320)=F 


54.7 ±8.0 


6.4 ± 1.2 




X 






X 


r/f2(1270) 


0.56 ±0.40 


X 


X 


X 






X 


7r0f2(1270) 


14.8 ± 1.9 


X 


X 


X 






X 


7r°f^(1525) 


0.38 ± 0.06 


X 


X 


X 






X 


KIK{I = 0) 


1.0 ±0.3 


0.07 ±0.24 




X 






X 


K^K(J=1) 


0.08 ±0.32 


0.7 ±0.2 




X 






X 


7?(1440)7r+7r- 


2.78 ±0.19 


S 


S 


X 


s 


X 


X 



x: Channel forbidden from this initial state. 
S: Channel suppressed by dynamical effects. 
-: Value not available. 



TT+TT , 7r°7r°, K+K , KgKg and KgKi channels were unchanged, was obtained with a very good fit 
to the data x^/iV = 10.8/11. The branching ratios for the (fm and (j)ri channels are also hsted on 
Table 29. 
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As a final stage of the analysis the measurements by the Crystal Barrel experiment [215] for the 
Trry, t]t], ttt]', rjrj', cow, lot] and cor]' channels in liquid and 12 psTP gaseous H2 targets were included. 
As these channels involve just two initial pp states, the two annihilation frequencies at different target 
densities allow the branching ratios to be determined and the value of x^/N remains unchanged. 
Again values of the resulting branching ratios are given in Table 29. As we have commented earlier, 
and as can be seen from Table 29, the values of BR(c/i, ^P2) for the annihilation channels included 
in this analysis are all consistent with zero. Fixing the values for BR(c/i, ^P2) = 0, and repeating 
the fit to the data gave an equally good fit to the data with x^/^ = 12.5/18. The corresponding 
values of BR(c/i, ^Pq) are given in the last column of Table 29. 

The p°7r° and p°T] two-body final states are produced solely from the ^Si and ^Pi states of 
the pp system. The Asterix measurements [121, 125,297] of S- and P-state annihilation frequen- 
cies for these channels (Sec. 6.4.5 and Table 20) can be used to derive the corresponding branch- 
ing ratios. The S-state annihilation frequency is derived from data measured in a gaseous hydro- 
gen target of density psTP and using Eq.(6.4) is related to the branching ratio by AF(/9"7r°)s = 
jE{'^Si, /9sTp)BR(p'^7r'^, ■'Si) and correspondingly for the p'^t] channel. The P-state annihilation 
frequency is derived from data obtained with an atomic L X-ray trigger (Sec. 6.4). Using Eq.(6.9) 
then gives AF(yo"7r'^)p = -^BR{p'^ tt'^ , ^Pi). The branching ratios obtained using these equations 
for the S- and P-state annihilation frequencies listed in Table 20 are given in Table 29. It should 
be noted that the S- and P-state annihilation frequencies are calculated neglecting interference ef- 
fects [121] from a partial wave fit to the Dalitz plot. In a similar way, branching ratios for ^Sg and 
•^Si states with 1 = and 1=1 can be obtained from Table 28 for the K*K channel and also for the 
7r°f2(1270), 77f2(1270) and 7ra2(1320) channels. These values are again Usted in Table 29. 

After the results of Table 29 were obtained, a somewhat similar analysis was made by the Obelix 
collaboration [284]. An essential difference from the present work was that they varied the en- 
hancement factors at STP and liquid H2 densities to give a best fit to the data whilst those for target 
density 0.005 psTP were fixed _E(^"''+^iyj, O.OOSpsTp) = 1- In the present work the enhance- 
ment factors were fixed at values obtained from a cascade calculation (See Sec. 6.3 and Table 12.) 
Two body branching ratios obtained by them for the (pn, K*K, pw, 7rf2(1270) and 7ra2(1320) chan- 
nels are in reasonably good agreement with those obtained in the present work (Table 29). For 
the p7r,K*K,7rf2(1270) and 7ra2(1320) channels they were, in addition, able to obtain two-body 
branching ratios for P-wave atomic states. 

The Obelix collaboration reported some interesting ratios of annihilation frequencies [160] which 
are reproduced in Table 30 and compared to the results from Table 29. 

Table 30: Ratios of annihilation frequencies measured by the Obelix collaboration [160]. 



Frequency ratio from [160] from Table 29 pp 

0.19 ±0.08 ^So 
^Pi 

^ 

AF(pp^K*K,J = 0) _ [ 0.17 ±0.03 0.06 ±0.06 % 
AF(pp^K*K,/=l) \ 0.20 ±0.01 

0.026 ±0.005 ^So 

^P2 



AF(pp^ K*K,J= 1) 
AF(pp^K*K,/ = 0) 



0.40 ± 0.07 
0.25 ±0.02 
0.10 ±0.01 



AF(pp^ f^(1525)7r°) ^ 
AF(pp ^ £2(1270)77") 



0.028 ±0.006 
0.026 ±0.003 
0.051 ±0.020 
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6.6 Antineutron annihilation on protons 

The study of np annihilation at very low li momenta into two-body final states gives access to the 
isospin 1=1 part of the annihilation operator. Such experiments were carried out by the ObeUx 
collaboration; methods and results have been reviewed in detail in [11]. Of interest here are cross 
section measurements for two-body final states. Data are reported for rip 7r°7r+ , KsK+, (f>Tr~^, 
WTT"*", r?7r+ and K*°K+, in the momentum range from 50 to 405MeV/c. Several partial waves can 
contribute to the final states; their momentum dependence can be fixed from the Dover-Richard 
potential model [300]. A system of linear equations can be formulated (and solved) which should 
determine hadronic branching ratios for the reactions considered, for S- and P-state annihilation. 
Only two results are given quantitatively: 

AFC'Si ^ 7r+7r°) = (3.1 ± 0.5) x 10"^ , 
AF(3Si ^ K+K°) = (1.3 ± 0.2) X 10"^ 

To include these results in further analyses, equivalent branching ratios are calculated. There are, 
of course, no enhancement factors due to an atomic cascade, only the normalisation 4/3 has to be 
applied (since only 3/4 of all np systems have spin 1) and we allow for the undetected Ki to arrive 
at branching ratios which can be compared to those of Table 29. 



BR(3Si ^ 7r+7r°) = (4.1 ± 0.7) x 10"^ 
BR(3Si ^ K+K°) = (3.5 ± 0.5) x 10 



(6.14) 



The collaboration also reported some ratios of annihilation frequencies [153]. In particular they 
find, integrating over the full momentum range, 

=0.63 ±0.16. (6.15) 

AF(7r+ry) 

The ratio (f)Tr^ /lutt^ was determined for the full momentum range and for three momentum bytes 
[154]. The results are reproduced in Table 3 1 . 

Table 31 : Cross section Ratio for np (/)7i-+/a;7r+. 



Momentum ( MeV/c) cr(np TT^(p)/a{np tt+w) 
50 to 405 0.075 ± 0.008 

50 to 200 0.100 ±0.017 

200 to 300 0.074 ±0.009 

300 to 405 0.062 ± 0.009 



6.7 Antiproton annihilation in D2 

The determination of annihilation frequencies for p stopping in D2 requires some remarks so as to 
appreciate the meaning of the results. We discuss these measurements considering as an example 
the two annihilation modes 

(A) pd^TT-wp, 

(6.16) 

(B) pd^7r°a;n. 
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In these reactions, the surviving nucleon can have large momenta; mostly one is interested in deter- 
mining the frequency for annihilation on a quasi-free nucleon. A nucleon not participating in the 
reaction will survive the process with a Hulthen momentum distribution. At sufficiently low mo- 
menta, the nucleon can thus be seen as a spectator particle. The experimental difficulty is then to 
determine the fraction of all annihilations in which one of the nucleons acts as spectator. 

In bubble chamber data, this task is easily met, at least for annihilations on neutrons. The incom- 
ing antiproton leaves a visible track. If a short secondary track is observed, it is likely to be a proton; 
its momentum is determined by its curvature or its range. A cut at 250MeV/c was mostly applied to 
select events with a spectator proton. Protons with momenta below 80MeV/c leave no visible track 
and produce events with an odd number of tracks. These events are retained. In events without short 
track, the proton had a long track which requires a large momentum and was not a spectator proton. 

In electronic counter experiments, the acceptance for low-momentum protons has to be well 
controlled, which is experimentally challenging. The Obelix collaboration determined annihilation 
frequencies in D2 only when using a gas target. Stopping of protons is then not a problem. The 
Crystal Barrel collaboration [209] measured first reaction (B), using the all-neutral final state, with 
all five photons (uj 7r°7 and the two tt" then decaying to 77) and the neutron detected or recon- 
structed as missing particle. Thus reaction (B) was measured over the full kinematic range. Then it 
was assumed that the proton momentum distribution in (A) is the same as the neutron momentum 
distribution in (B). Thus the Monte-Carlo proton momentum distribution is very close to the true 
one, and the cut, at lOOMeV/c, leads to the same fractional loss in Monte Carlo and real data. The 
KEK collaboration observes only one of the particles, a 7r° or rj. A finite momentum transfer to the 
spectator nucleon leads to a broadening of the meson recoil momentum distribution and, eventually, 
to a loss of the signal. A precise definition of a spectator nucleon is not possible, and the data could 
be difficult to interpret. 

The Asterix collaboration determined the annihilation frequency only for tt+tt^ and K+K~. A 
cut in coUinearity is equivalent to a cut in the neutron momentum distribution. It was checked that 
the final result did not depend on this cut. 

Finally we comment on measurements over the full kinematical range. Bizzarri et al. [83] have 
measured the frequency of annihilation into KgKi and KgKg. The purpose of this study was to 
demonstrate the importance of induced P-state annihilation. Indeed, their KgKg/KsKi ratio is about 
20 times larger than in H2. It is hence not suited for a comparison with the other results. 

6.7.1 Two-body annihilation frequencies 

Measurements of two-body annihilation frequencies for the annihilation at rest of antiprotons on 
deuterons are hsted in Table 32. The measurements by Chiba et al. [108] for the channels 7r°p°, 
'K^p~, rjp^ and r]p~ have not been included. As well as indicating the detection technique and 
method for event identification, as in Sec. 6.4, it is also stated whether the measured annihilation 
frequencies are normalised to the number of pd (d), pp (p) or pn (n) annihilations (see Sec. 6.1). 
Ratios of annihilation frequencies for two-body channels from pd annihilation at rest are presented 
in Table 33. 

A number of checks can be made on the consistency of the annihilation frequency data as has 
been discussed by Batty [287]. By charge independence [91, 255] (see Sec. 4.4.3), for pd atoms, 

AFp(7r+7r-,n,p) = ^AF^ (77-77°, p, p) + 2AFp(7r°7r°, n, p) . (6.17) 

For Uquid D2 using Eq. (6.3) with the measurement of AFd (tt^tt" , n, liq.) by Amsler et al. [177] 
and the other measurements [91,276] also from Table 32, gives a value of (4.2 ± 1.2) x lO"^ 
for the left-hand side and (6.47 ± 0.51) x 10"^ for the sum of the two terms on the right-hand 
side of the equation. This latter value becomes (6.27 ± 0.58) x 10~^ when the measurement of 
AFd(7r~7r°, p, liq.) of Abele et al. [209], together with Eq. (6.3), is used. Whilst these two latter 
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values for the right-hand side are in good agreement, there is a 2 sigma discrepancy between the left- 
and right-hand sides of Eq. (6. 17). 

A further check can be made by calculating values of 

r = AFd(7r"7r",p,liq.)/AFd(7r+7r",naiq.) , (6.18) 

from the annihilation frequencies of Table 32 for Uquid D2, using Eq. (6.3) where necessary, and 
comparing them with the direct measurements of Table 33. Using the measurements of Refs. [91, 
276], together with Eq. (6.3), gives r = (1.58 ± 0.49). Replacing the tt^tt" frequency measurement 
of ref. [276] by that of ref. [209] then gives r = (1.50 ± 0.48). Neither of these deduced values is in 
agreement with the direct measurements of Table 23, with the possible exception of the value of r 

Table 32: Annihilation frequencies for pd interactions at rest normalised to the number of pd (d), pp 
(p) or pn (n) annihilations. 



Chan. Dens. Ann. Freq. Ref. Type 


Dens. Ann. Freq. Ref. Type 


Liq. (5.90 ± 0.32) X 10-'' 1177] CE(d) 
Liq. (1.23 ± 0.55) X lO"''! 108] /CR(p) 


Liq. (6.70 ± 0.80) x lO-'* [209] C E (d) 


TT+TT-n Liq. (4.20 ± 1.20) X 10-3 [91] BE(p) 


1 (2.01 ± 0.27) X 10-3 [120] y4E(d) 


Liq. (3.61 ± 0.49) x lO-'"* [209] C E (d) 

TT 7r°p 

Liq. (2.02 ± 0.31) x 10"^ [108] /C R (n) 


Liq. (8.80 ± 1.00) x 10-^ [276] B E (n) 


7r°r?n Liq. (2.46 ± 0.12) x 10"^ [177] C E (d) 


Liq. (5.7 ± 1.8) X 10-^ [108] /C R (p) 


TT-rjp Liq. (4.06 ± 1.00) x 10"* [209] C E (d) 




7r-?7'p Liq. (2.98 ± 1.52) x 10"^ [209] C E (d) 




Liq. (6.04 ± 0.69) x lO'^ [209] C E (d) 
TT-up Liq. (6.00 ± 1.00) X 10-3 [68] BE(n) 
1 (4.97 ± 0.89) X 10-3 [136] O E (d) 


Liq. (4.10 ± 0.80) X 10-3 [78] BE(n) 
Liq. (13.2 ± 4.3) x IQ-^ [276] B E (n) 


7r"^'u Liq. (4.LS ± 0.3!)) x K)-"' |209| C E (d) 


Liq. (1.9 ± 0.9) X 10-"' 1 108| Kl R (p) 


rjujii Liq. (6.71 ± 0.81) x 10-^ [209] C E (d) 


Liq. (4.9 ± 1.2) x 10-^ [108] K R (p) 


K-Ksp Liq. (1.42 ± 0.36) x 10-^ [209] C E (d) 




KsKsH Liq. (3.60 ± 1.00) X 10-5 [83] BE(d) 




KsKin Liq. (3.60 ± 0.40) x 10-^ [83] B E (d) 




K+K-n 1 (7.30 ± 1.60) x 10-'* [120] A E (d) 




Liq. (8.80 ± 2.20) X 10-* [80] B E (n) 

TT ^p 

1 (6.62 ± 0.49) X 10--* [136] O E (d) 


Liq. (9.20 ± 1.10) X lO-'* [82] BEin) 



measured by Angelopoulos et al. [301] which is much larger than the other directly measured values 
[91,1 00] . These discrepancies between the annihilation frequencies and also between their ratios are 
discussed further in the following section. We also note, from Table 32, that the KEK measurements 
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Table 33: Ratios of annihilation frequencies for pd annihilation at rest 



Channels 


Density 


Ratio 


Ref. 


Type 


7r~7r° p/7r~7r+ n 


Liq. 


0.68 ± 0.07 


[91] 


BE 


7r~(7r°)p/7r~(7r+) n 


Liq. 


0.70 ± 0.05 


[100] 


TR 


7r~(7r°)p/7r+(7r~) n 


Liq. 


0.55 ± 0.05 


[100] 


TR 


7r~7r'^ p/7r~7r+ n 


Liq. 


2.07 ±0.05 


[301] 


jCR 


K+K~ n/TT+TT" n 


1 


0.36 ± 0.08 


[120] 


AE 


K+K~ ti/tt+tt" n 


1 


0.27 ±0.02 


[129] 


OE 


u)iT~ p/^7r~ p 


Liq. 


7.3 ± 1.5 


[82] 


BE 



[108] are generally in disagreement with other work. These measurements were therefore excluded 
from the forthcoming analyses. 

6.7.2 Analysis of two-body annihilation frequencies for pd annihilation at rest 

Before discussing the analysis of two-body annihilation frequencies in D2, a few sentences of warn- 
ing seem appropriate. A final separation of truly spectator-like events would require a full partial- 
wave analysis of the three-body final state. Such fits have never been performed. The largest con- 
tribution from 3-body Pontecorvo reactions (Sec. 6.7.3) is expected from A(1232) recoiling against 
a pion in the reaction pd ttttN. Hence we should expect deviations from simple models. A 
further problem hes in discrepancies between measurements of the same annihilation frequency ob- 
tained from different experiments. As we have seen, such discrepancies also occurred in data on pp 
annihilation frequencies which are conceptually much easier. 

An analysis of the annihilation frequency information for pd annihilation at rest has been made 
by Batty [287]. The annihilation frequencies were analysed using Eq. (6.8). Calculated annihilation 
frequencies were fitted using the least squares method to give a best fit to the experimental measure- 
ments by varying the fraction of P-state annihilation fp^„„ (p). The branching ratios BR(c/i, ^'^"'"^Lj) 
were fixed at values obtained from the analysis of pp annihilation frequencies by Batty [282, Table 
6] using either the DRl or DR2 models. For the nn system, in addition to Eqs. (6.10) and (6.11), 
charge symmetry can be used to give 

BR(7r-7r°, ^Si) = 2BR(7r+7r-, ^Si) . (6.19) 

Exphcit equations for the annihilation frequencies for the tt+tt" , 7r°7r°, tt" 7r° , K+K~ , KgKg and 
KgKi channels in terms of branching ratios and the fraction of P-state annihilation fp^„„ (p) for pd 
annihilation have been given by Batty [287, Eq. (9)]. As a first step, single annihilation frequencies 
for these channels were fitted to determine whether the best fit values of fp^„„ (p) lay, within errors, 
in the physically acceptable range < fp^„„(p) < 1- This preHminary analysis indicated problems 
with the two measured values [209,276] of AF(7r^7r°, p, liq.) and the measured value of r (see 
Eq. (6. 18)) determined by Angelopoulos et al. [301 J which were discussed previously in Sec. 6.7.1. 
The measured KgKg annihilation frequency [83] also gave values of fp^^_^(liq.) > 1. This result 
included the full proton momentum range and cannot be compared to data selecting events with a 
spectator proton. 

These 4 values were therefore not included in the fit to the measurements of Tables 32 and 33 
for these channels. A least squares fit to the remaining 6 measurements for a liquid target with 
Sp = 0.571, as discussed in Sec. 6.1, and branching ratios obtained with the DRl model, gave 
fPann(lict-) = O-^O ± 0.02 wlth a chi-squared per degree of freedom x^/^ = 10.9/5. Much of 
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the contribution to comes from the measurement [100] of r = 0.70 ± 0.05. Omitting this 
measurement gives little change in the value of fp^„„ (liq-) = 0.42 ± 0.02, but chi-squared decreases 
significantly and x^/^ = 4.8/4. Using a value of Sp = Sn = 0.5, i.e. assuming equal annihilation 
on neutrons and protons, gives fp.^j^_^(liq.) = 0.48 ± 0.02 and = 12.0/5. Repeating the analysis 
with branching ratios obtained using the DR2 model, gave the same values for fp^_^^(liq.), within 
errors, whilst decreased at most by 0.5 in absolute value. This result is not surprising, since 
the values of BR(K+K", and BR(K°K", ^■S'+iLj) obtained with these two models are 

almost identical. The same is also true for the values of BR(7r+7r^ , ''Si) and (^BR(7r+7r~ , ^Pq) + 
^BR(7r+7r~, ^P2)) which appear on the right hand side in Eq. (6.8). 

Repeating the analysis using the four measurements obtained with a gas target gave fp^„„ (pstp ) = 
0.34 ± 0.04 with xV^V = 5.7/3. With Sp = s„ = 0.5 values of fp^„„(psTp) = 0.34±"o.04 and 
x'^/N = 9.5/3 were obtained. For branching ratios obtained with the DR2 potential and Sp = 0.571 
abestfitwas obtained with fp,^„,,(/7sTp) = 0.34±0.04and x^ = 5.3/3. 

As we have already mentioned, P-state annihilation on a nucleon can occur from S-states of a pd 
atom, due to the size and momentum distribution of the nucleons in the deuteron. This is sometimes 
referred to as "induced" P-state annihilation. Bizzarri [84] first considered this process in terms of 
the spectator model, using a semi-classical method based on the Fermi motion of the nucleons in 
the deuteron and an impact parameter relative to the deuteron centre of mass. A calculation then 
gave the amount of P-state annihilation from atomic S-states to be 0.4 times the rate of S-state 
annihilation and predicted that there should be a relatively large amount of S-state annihilation from 
atomic P-states. 

If ass„„„ and app,^^^^ are used to denote the fraction of S- and P-state annihilation from atomic 
S- and P-states respectively, then the fraction of P-state annihilation is given by fpj,„„(/o) = (1 — 
fp(/9))(l — asSann) + fp(p)oipPa„„- Usuig thc values of fp(psTp) = 0.75, fp(liq.) = 0.40 obtained 
in Sec. 6.3.2, and fpa„„(/5sTp) = 0.34 ± 0.04 and fp„„(liq.) = 0.40 ± 0.02, given above, then 
values of asSa„„ = 0-53 ± 0.06 and apPa„„ = 0.30 ± 0.07 are obtained, where the errors shown are 
purely statistical. The fitted values of asSa„„ ^uid apPa„„ are sensitive to the values of fp(/9sTp) and 
fp(liq.) which were calculated using a cascade calculation as described earlier in Sec. 6.3.2. These 
were obtained assuming the Stark mixing parameter Kq = 7.6. If instead Kq = 1.0 is used, then the 
cascade calculation gives fp(/9sTp) = 0.92 and fp(liq.) = 0.64 (see Fig. 16) and a best fit then gives 
asSann = 0.46 ±0.11 and app.^^^^ = 0.32 ± 0.05. Whilst the numerical values for ass„„„ and app^^^^ 
given by these estimates should be treated with caution, they nevertheless do indicate the importance 
of "induced" S- and P-state annihilation. 

6.7.3 Pontecorvo reactions 

So far in the discussion of pd annihilation we have concentrated on the production of two-meson final 
states in which the antiproton annihilates with a single "quasi-free nucleon" in the deuteron whilst the 
other non-participating nucleon is a "spectator" particle. It was pointed out by Pontecorvo [302] in 
1956 that there is a further class of annihilation reactions in which, through a three-body interaction 
involving both nucleons, the final state consists of a single meson and a baryon. These Pontecorvo 
reactions can conveniently be considered in three classes. 

The first of these classes has a final state consisting of a single meson and a nucleon, e.g., pd 
7r~p. These reactions, in the simplest model, can be considered in terms of a two-step (rescattering) 
process, for example pd — > 7r~ (7r+n) 7r~p in the above case. Other possible reactions involve 
the production of 77, ry', w and (p mesons. The second class consists of cases where one of the mesons 
interacts with the second nucleon to form a nucleon resonance, e.g., pd — > 7r~(7r+n) 7r~A+. An 
alternative to the two-step involves mentioned above is the fireball model which pictures a compound 
system ("fireball") formed by the participating three antiquarks and six quarks, decaying statistically. 
This latter model predicts significant production of the third class of channels with open strangeness, 
e.g., pd AK°, E°K° or E~K+. These latter reactions have recently been observed by the Crystal 
Barrel Collaboration [206]. 
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A compilation of measured frequencies for Pontecorvo reactions in pd annihilation is given 
in Table 34. In general the results are in moderate agreement with each other in those cases 
where several measurements have been made for a single channel. In particular AF(7r~p) = 
(1.38 ± 0.07) X 10"^ with x^/N = 4.1/3, where N is the number of degrees of freedom. Similarly 
AF(7r°n) = (7.13 ± 0.71) x lO"'^ with = 0.6. 

A recent comparison of these results with the predictions of dynamical (two-step) and statistical 
(fireball) models has been made by the Crystal Barrel Collaboration [206]. They show that the rates 
AF(AK°), AF(S°K°) and in particular the ratio AF(I]OkO)/AF(AKO), exceed by a large factor 
the values predicted by two-step models, but are in general agreement with the predictions of the 
fireball model. 

Table 34: Frequencies for Pontecorvo reactions in pd annihilation. Some final states are produced 
via the same isobar. In these cases we give the isobar and the inverse squared Clebsch-Gordan 
coefficients to facilitate the comparison of different results. K" stands for the sum of and Ki 
modes. 



Chan. 


Dens. 


Isobar 


C.G.-2 


Ann. Freq. 


Ref. 


Type 


7r~p 


1 


ttN 


2 


(1.40 ±0.70) X 10-5 


[120] 


AE 


7r~p 


1 


ttN 


2 


(1.20 ±0.14) X 10-5 


[130] 


OE 


7r~p 


1 


ttN 


2 


(1.46 ±0.08) X 10-5 


[152] 


OE 


7r~p 


Liq. 


ttN 


2 


(0.90 ±0.40) X 10-5 


[303] 


BE 


TT^n 


Liq. 


ttN 


2 


(7.03 ±0.72) X 10-6 


[180] 


CE 


TrOn 


Liq. 


ttN 


2 


(1.03 ±0.41) X 10-5 


[304] 


ICR 


T] n 


Liq. 






(3.19 ±0.48) X 10-6 


[180] 


CE 


T] n 


Liq. 






< 8.94 X 10-6 


[304] 


ICR 


u n 


Liq. 






(2.28 ± 0.41) X 10^5 


1180] 


CE 


T]' n 


Liq. 






< 14 X 10-6 


[180] 


CE 


(j) n 


1 






(3.56 ±0.35) X 10-6 


[157] 


OE 




Liq. 


ttA 


9/2 


(2.21 ± 0.24) X 10-5 


[177] 


C E AO ^ ttO n 




Liq. 


ttA 


3 


(4.67 ± 1.66) X 10-5 


[304] 


/CE 




1 


ttA 


9 


(1.22 ±0.20) X 10-5 


[152] 


O E AO ^ TT- p 


TT- A+ 


1 


ttA 


9 


(1.01 ±0.08) X 10-5 


[152] 


O E A+ ^ ttO p 


77 AO 


Liq. 






< 6.49 X 10-5 


[304] 


/CE 


ARO 


Liq. 






(2.35 ± 0.45) X 10-6 


[206] 


CE 


SO RO 


Liq. 






(2.15 ±0.45) X 10-6 


[206] 


CE 



6.8 Discussion 

A compilation of two-meson annihilation frequencies for pp annihilation at rest has been presented, 
together with an analysis of some of the data in terms of branching ratios BR(c/i, ^^~^^Lj) for the 
channel ch from specific ^'^+^Z/j states of the pp system. The fraction of P-state annihilation as a 
function of H2 target density has also been determined. 
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A very wide range of measurements are available with some 187 annihilation frequencies cover- 
ing about 30 different annihilation channels at 7 target densities. Unfortunately in some cases there 
are very significant discrepancies between measurements for the same channel and target density; in 
one extreme case by a factor of about 40. In most of these cases, the discrepant measurements were 
obtained using peaks in the inclusive recoil energy spectrum to identify the second meson. These 
experiments should generally be regarded as less reliable than those in which the complete two-body 
event is recorded and identified. 

A particularly unfortunate discrepancy is that between the recent measurements for the pp — > 
7r°7r° channel in liquid H2, made by the Crystal Barrel and Obelix collaborations. Both experiments 
fully reconstruct the tt^tt^ events but the measured annihilation frequencies differ by a factor of 
about 2. This discrepancy is discussed in some detail in Sec. 6.4.3 where the internal and external 
consistency of the Crystal Barrel results relating to their annihilation frequency for the 7r°7r° channel 
has been emphasised. In the absence of similar checks for the Obelix measurements, the Crystal 
Barrel results are regarded as being more reliable. 

For some channels (p^n^, p^r], KgKi, tt+tt" and K+K~) produced from a liquid H2 target, there 
is good consistency between several measurements of the annihilation frequencies and the weighted 
mean value is given in the text. In Tables 26, 27 and 28 we presented a selection of annihilation 
frequencies which is used to derive branching ratios from specific atom states. 

The target density dependence of annihilation frequencies has been discussed in Sec. 6.2, where 
it was shown that the effects of Stark mixing are important. In particular, where there is a large 
contribution from Stark mixing, an atomic fine-structure level with a large annihilation width will 
contribute more to annihilation than would be expected from its statistical weight only. These devia- 
tions, which are particularly important for the ■'Pq state, can be described |282] in terms of enhance- 
ment factors, which are a function of the initial state and target density p. In Eq. (6.4) the annihilation 
frequency AF(c/z, p) is written in terms of the branching ratios BR(c/i, enhancement fac- 

tors Ei^^^^Lj, p) and the fraction of P-state annihilation fp [p). Values of the enhancement factors 
(Table 12), have been calculated from a cascade calculation using a modified version of the Borie 
and Leon model [288], whose parameters Kq and ksTK are obtained from a least-squares fit to pp 
atomic X-ray yields as a function of H2 target density. 

Using these enhancement factors, the annihilation frequencies Usted in Tables 26, 27 and 28 were 
fitted using the least-squares method by varying the fraction of P-state annihilation and branching 
ratios. The best fit values are given in Table 24 (See also Fig. 14) and Table 29. Branching ratios for 
the p°7r° and p^r] channels, also listed in Table 29, were obtained directly from measurements of S- 
and P-state annihilation frequencies. Branching ratios for ^Sq and ^Si states for K*K and K2K, and 
for 7r°f2 (1270), ?7f2 (1270), 7r°f2 (1525) and 7ra2 (1320) channels were obtained in a similar way from 
the results of partial wave analyses and are also Usted in Table 29. Branching ratios are independent 
of atomic physics effects and are the quantities to be compared with the predictions from models of 
the annihilation process and two-meson production. 

A problem with determining branching ratios by using Eq. (6.4) to fit experimental annihilation 
frequencies, should be mentioned. For a particular annihilation channel which occurs from several 
(say, n) annihilation states, measured annihilation frequencies for at least n different target densities 
(or in coincidence with atomic L X-rays) are required. However the branching ratios can only be 
determined for two sets of cases. Firstly where there is only one S- and one P-state (e.g., "^Si and 
^Pi for wtt) and the fraction of P-state annihilation is known or can be determined from the data. 
Secondly for the case where there is more than one P-state (e.g., '^Po and "^P2 for tt^tt^) and the 
enhancement factors for these states differ significantly at a given density for the relevant 
The latter condition generally implies that data for a liquid H2 target is available and that one of the 
initial states is ^Pq which has a large enhancement factor at liquid H2 density. 

To illustrate this point we consider the reaction pp (j)co which occurs from the initial states 
^Si, ^Po, ^Pi and^P2. In this case, as can be seen from Fig. 15, the enhancement factors for the ^Pi 
and ^P2 states are similar and the two branching ratios cannot be separated by a fit to the annihilation 
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frequencies. In practice, as can be seen from Eq. (6.4), the quantity determined will be 

^E{3p,,p)BR{^uj, 3Pi) + |-ij(3p2,p)BR(0u;, ^P^) , (6.20) 



where Ei^^+'^Lj, p) is the enhancement factor averaged in some way over the range of target den- 
sities. 

A compilation of two-meson annihilation frequencies from pd annihilation at rest shows that 
only a limited range of data is available and that much of it is inconsistent. A fit to the rather limited 
range of pd atom X-ray data, assuming that the value of the cascade parameter Kq is the same as 
in the pp case, gives values for the fraction of annihilation (fp (p)) from atomic P-states in gaseous 
(psTp) and liquid H2. For the case of pd atoms the values of the enhancement factors are expected 
to be close to one. This is confirmed by a cascade calculation. 

Using the method described earher for pp atoms, a fit is made to a selected set of annihilation fre- 
quencies for pd annihilation at rest with branching ratios given by the analysis of pp data discussed 
earlier. Values for the fraction of P-state annihilation (fp^^^^ [p)) are obtained for gaseous (pstp) and 
hquid H2. The difference in value between fp(p) and fp^^^ (p) is explained using a simple model in 
terms of a significant amount of "induced" S- and P-state annihilation. This latter conclusion is in 
agreement with the predictions of Bizzarri [84]. 
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7 Dynamical selection rules 

Nucleon-antinucleon annihilation reveals a great surprise: some processes which are perfectly com- 
patible with the rules of quantum-number conservation summarised in Sec. 4 (parity P, charge 
conjugation C, G-parity, among others) are observed only at a much reduced rate. These are dynam- 
ical selection rules. This property has often been emphasised; in particular Dover stressed the link 
between dynamical selection rules and the dynamics of quarks and gluons in NN annihilation [5]. 



7.1 The pn puzzle 

The best known example of a dynamical selection rule is the so-called pn puzzle. Figure 17 shows 
the Dalitz plot of the reaction j5p Tr+vr^vr" in annihilation at rest [190]. The three bands in 
Fig. 17 show production of p+Tr^, p^7r+, and p°7r" intermediate states. The intensity distribution 
along, e.g., the p+ band gives directly the p+ — > vr+vr" angular distribution in the p+ rest frame. 
On the right-hand side of Fig. 17, in (a), the angular distribution is plotted as a function of cosO 
where O is defined as angle between the p+ flight direction and the direction of the 7r+, in the p'^ 
rest frame. The decay angular distributions give access to the quantum numbers of the initial state 
feeding the final state. Table 35 lists the initial states from which annihilation into pir is allowed, 
the p — TT orbital angular momentum, and the p decay angular distribution when all interferences 
are neglected. The sign of the amplitudes is also given {p'^ii~ — p~tt^ stands for a negative sign 
between the peak amplitudes for p+7r~ and p~Tr~^ production). At lines of identical tttt masses, i.e., 
for identical Breit-Wigner phases, the interference is constructive for annihilation from isoscalar and 
destructive for annihilation from isovector initial states. 




Figure 17: The tt+tt tt" Dalitz plot in pp annihilation at rest, and p+ (a), p (b) and p° (c) decay 
angular distributions. 
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A first look at the decay angular distributions reveals two strong peaks above a sin^ Q distribu- 
tion. The peaks can be understood by inspecting the full Dalitz plot: the p"*" and p~ bands cross 
in the lower left comer of the Dalitz plot. The amplitudes for the two processes pp — > p+7r~ and 
pp p^7r+ interfere constructively and lead to marked deviations of the observed angular distribu- 
tion from the expected one. The crossing of the p+ band with the p~ or p° band leads to an increase 
of the intensity by a factor 4 because of quantum mechanical interference (the amplitudes are added). 
The dominant sin^ O distribution beneath the two interference peaks is due to strong contributions 
from the 2^+1, 2S+i^^ _ i,3g^ initial state. The ^'^Sq pp initial state generates a cos^ & which 
would lead, with a fraction of the sin^ 6 part, to a constant plateau beneath the sin^ O distribution. 

Table 35: Angular distributions for pN fm annihilation. The atomic states are represented as 
2j+i,2s+i^^ = ^'^Si with I,S,L,J being isospin, spin, orbital and total angular momenta; £ is the 
orbital angular momenta between p and n, 6 the angle between the direction of the more positively 
charged pion from p decay with respect to the p direction of flight. Forbidden transitions are marked 
by an x. The signs indicate constructive and destructive interference. 
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A second look reveals an even more important aspect: the three p peaks have nearly the same 
strength. Now, p°7r° production is forbidden from initial states with charge conjugation C = 
From initial states with C = — 1, the rates for the three charge modes are fixed by Clebsch-Gordan 
coefficients, and are identical. In the limit of exactly equal strengths, only initial states with C = — 1 
contribute. Because G = (—1)^ C, C = — 1 and an odd number of pions entails / = 0. The three 
angular distributions give evidence for a small reduction of p^Tr" relative to p+7r~ or p^7r+. The 
smallness of the reduction shows that states with C = +1 do contribute, but only by a small fraction. 

The small flat contribution to the decay angular distributions can be due to a cos^ O part or due 
to annihilation from atomic P-states; both contributions are evidently small. 

A partial- wave analysis quantifies these observations, determines masses and widths of contribut- 
ing resonances and gives fractional contributions from the different pp initial states to the 7r+7r~7r'^ 
final state. Here, only the annihilation frequencies for p7r production are given (adapted from Ta- 
ble 29). 

The pTT channel is produced from the ^J^+i.as+i^,^ = i>3g^ gp initial state with a branching 
ratio (63.6 ± 3.3) x IQ-^, and with branching ratio (6.0 ± 2.0) x lO'^ from the ^'^Sq state. The 
ratio 

^|£j|^ = 0.094 ± 0.032. (7.1, 

is rather small. There is no obvious explanation of the suppression of pn from the ''■^So state. The 
preference for annihilation via the / = initial state is a dynamical selection rule. In the literature, 
it is known as the pn puzzle. 

This selection rule is partly confirmed in annihilation from atomic P-states. Recently, the Obelix 
collaboration reported a pressure-dependent study of pp annihilation into 7r^7r~7r°, K"'"K~7r°, and 
K^'^KgTr^. Figure 18 shows the tt+tt" and tt^tt*^ invariant mass distributions for pp annihilation at 
rest in liquid H2, with ~ (87.5 ± 2.0)% S-state, ~ (12.5 ± 2.0)% P-state capture and in gaseous H2 
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Figure 18: Invariant mass square distributions for pp annihilation at rest using H2 targets of three 
different densities [160]. Left: pp ^ tt+tt-tt" , centre: pp ^ K+K-K° , right: pp ^ K^K^tt^f . LH: 
liquid H2; NP: atmospheric pressure and room temperature; LP: low pressure. 



at atmospheric pressure (with ^ (36 ± 3)% S-state, ^ (64 ± 3)% P-state capture, and at very low 
pressure, with ^ (11 ± 7)% S-state, (89 ±7)% P-state capture. 
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The comparison of the upper two plots showing the tt+tt" and t:^t:^ (squared) invariant mass 
distributions in liquid H2 confirms the observations of the Crystal Barrel data: the peak height for 
the is nearly as large as the peak height of the charged p^: isoscalar initial states must dominate 
the reaction. Since there is S-state capture dominance, the most important contribution has to come 
from the Si protonium state. The two lower plots, from annihilation at low pressure, show a more 
significant signal. The is reduced by approximately 30%. This would indicate that 2/3 of all 
annihilations at low pressure are from the isoscalar, and 1/3 from the isovector component of the 
protonium atom. Inspection of Table 35 shows that the isovector P-states have a statistical weight 8, 
the isoscalar ones a weight 3. Taking this into account, the pir puzzle is present also in annihilation 
with a large P-state capture probability, even though less pronounced. More remarkable is, however, 
the result of the partial wave analysis that the ^'^Pi level decays only weakly into p-n with £ = 0: 

AF(pp^p±7rT,iPi,£ = 0) ^ ' 

This is a surprise: as seen in Table 5, ^''^Pi is forbidden to decay into most simple two-meson modes 
(tttt, KK, ttw, . . . ), and one would naively expect this channel to make full use of the p-K mode, 
in particular in S-wave {t = 0). The decay ^'^Pi — > pw with i = seems to be dynamically 
suppressed. 



7.2 Annihilation into a2 (1320)7r 

A further dynamical selection rule similar to the pn puzzle is observed in pp annihilation into 
a2(1320)7r. Figure 19 shows the most direct evidence for this rule, from Crystal Barrel data on 
the 7r"'"7r~27r° final state. In the 3-pion invariant mass distribution, clear evidence for the intermedi- 
ate a2(1320) decaying into pw is seen. The strength of the a2(1320) in the three plots is about the 
same. Now, annihilation into a^(1320)7r^ is allowed from both the ^'^Sq and the ^'^Si initial state 
while a2(1320)7r° is forbidden from ^'^Si. Both processes require i = 2 between a2(1320) and 
the pion. The fact that the rates for production of charged a2(1320) and neutral a2(1320) mesons 
are not too different entails that a2(1320)7r production is preferred from the isoscalar ^'^Sq state and 
suppressed from the isovector ^'^Si initial state. Again, there is no known reason for the suppression 
of one initial state compared to another initial state. The suppression is due to a dynamical selection 
rule. 




Figure 19: The n^2n , tt+tt^tt", and TT~2n invariant mass distributions from pp anniliilation into the 
7r+7r~27r° final state. There are two entries per event for the 7r+7r~7r° mass plot. All three distributions 
show a broad peak due to the a2(1320)7r production. The number of neutral a2(1320) is not visibly 
smaller than that of charged a2(1320). 
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The data shown above were not analysed for the different partial-wave contributions. They are 
shown here only since they visuaUse best the new dynamical selection rule. Partial wave analyses 
were performed on bubble-chamber data, with the a2( 1320) decaying into pw, rjn, andKK. The par- 
tial wave analyses confirm the qualitative observation discussed above: annihilation into a2(1320)7r 
is strong from the / = initial state (^Sq) and weak from the / = 1 component of the ^Si state, 
with a ratio, see Table 29, 

BR(3.3Si ^ a2^) ^ 6.4 ±1.2 

BR(i^iSo ^ aaTr) 82.0 ±12.0 ' ^ ' ^ 

The a2(1320)7r annihilation mode from the isovector component of the protonium atom is sup- 
pressed by one order of magnitude compared to annihilation from the isoscalar component. This is 
the same dynamical selection rule as the pir puzzle. 



7.3 Annihilation into tttt 

There is a related observation in pp annihilation into tttt and into KK. The branching ratios into KK 
and into tttt depend strongly on the initial state. In Table 36 the relevant branching ratios are listed. 
Branching ratios determined from rip scattering are included in the table. 

The ratio of KK versus tttt production is about 1 in annihilation from the Si state and less 
than 0.1 for annihilation from ^Pq. This result was interpreted as KK suppression in P-wave an- 
nihilation [305]. The (corrected) branching ratios in Table 36 show that, instead, tttt production is 
suppressed in annihilation from the '^Si state and not KK in P-wave annihilation. This interpretation 
is also favoured by a comparison of cross sections for annihilation in flight. The cross sections for 
annihilation into tttt and into pTT are very similar in size, while at rest the annihilation frequencies 
differ by one order of magnitude. Hence we interpret the change in the KK/tttt ratio as suppression 
of the TTTT frequency in the ^Si state. 



Table 36: Branching ratios in units of 10 ^ for pp and np anniliilation at rest into tttt and into KK, 
from Table 29 and Eqs. (6.14), (6.14). An x indicates that the transition is forbidden. 
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1.31 ±0.07 
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54.3 ± 1.5 


4.2 ±0.3 


X 


0.79 ±0.15* 





*The branching ratio for pp — > KgKs + KiKj is reduced by a factor of 2 due to Bose symmetry. 



Annihilation into tttt from the ^Si initial state proceeds via the isovector component of the pp 
system: again annihilation from the isovector part of the protonium wave function into two isovector 
particles is suppressed, by about one order of magnitude. 

Obviously, there is a common property in these three cases. Annihilation into two isovector 
mesons such as tttt, ptt and a2(1320)7r is strong from isoscalar initial states. In the language of 
SU(3)f flavour symmetry, to be described later in this review, the transition from an isoscalar pro- 
tonium state to two isovector mesons corresponds to a symmetric D coupling that is large. The 
transition from an isovector protonium state to two isovector mesons, described by an (antisymmet- 
ric) F coupling, is weak. 

7.4 Annihilation into KK 

The annihilation frequencies listed in Table 36 reveal another surprising fact: the annihilation fre- 
quencies into K+K~ and KgKi differ only by 10%. Assuming an isospin-blind annihilation po- 
tential, annihilation into K+K~ or into KgKi (depending on the relative phase of the 7 = and 
1=1 amplitudes) could be suppressed in the same way as the charge exchange reaction pp fin 
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is in a scattering situation. On the other hand, if only one isospin component of the protonium wave 
function contributes to the process, the two annihilation frequencies must be identical. 

Assuming that one isospin component dominates annihilation into a KK pair, it can be identified 
by comparison with the reaction lip K+K''. The branching ratio for the latter reaction should 
be either close to zero (for an isoscalar transition operator), or close to the sum of annihilation 
frequencies into K+K" and KgKi (2.78 ± 0.09) x 10"'^ if the transition operator is iso-vectorial. 
The experimental value (3.47 ± 0.53) x lO^'^ requires I — 1. Why is the annihilation reaction 

pp (^'^Si) ^ K+K" - K°K° , strong , (7.4) 

and 

pp (^-^Si) ^ K+K" + K"k" , weak? (7.5) 

There is no known reason for the suppression of one isospin channel with respect to the other chan- 
nel. It is a dynamical selection rule. 



7.5 Annihilation into K^K^tt^ and K±Ki7rT 

A similar selection rule is seen in pp — > K^KsTt^. Figure 20 shows (on the right) the Crystal Barrel 
data on this reaction. The most striking feature are the two K* bands (K^tt^ and KsTt^). The two 
bands cross but show little intensity at the crossing on the diagonal, where the two K* masses are 
identical. Obviously, the two amplitudes (for K**^ and K*^ production) interfere destructively. 

X 10-^ 




500 1000 1500 2000 500 1000 1500 2000 ^QQ 1000 1500 2000 

.III"' ,'-1"' 3 

[MeVV] [MeVV] X 10 

, [MeV'/c*] 

Figure 20: The experimental Dalitz plot for the reaction pp -> K^K^n^ (liquid H2). The experimental 
data (from Crystal Barrel) are shown on the right; on the left are Dalitz plots expected for annihilation 
from the various possible S-wave components of protonium. 

The following discussion is restricted to S-state capture. The final state K^KsTt^ does not have 
defined G-parity, and both isospin components may contribute to the annihilation process. There are 
four initial states which may contribute to the Dalitz plot in Fig. 20 (right). On the left. Fig. 20 shows 
the Dalitz plots expected for the four allowed initial states. The ^So initial state is characterised by a 
cos^ Q angular distribution, with high intensity at the ends with cos 9 = ±1. The ^Si state produces 
a sin^ Q distribution and the intensity vanishes at the ends. The Dalitz plot expected for the isovector 
and the isoscalar component of the protonium wave function differ by their interference pattern. 

It is obvious from the comparison of data and theoretical Dalitz plots that the two theoretical 
DaUtz plots with destructive interference between the two K*K amplitudes are the ones contributing 
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most to the final state. This observation is confirmed in partial wave analyses, see Table 22, where 
the contributions to K*K production are listed. The sum of the two partial waves with destructive 
interference exceeds the sum of constructive interference by one order of magnitude. 

This dynamical selection rule can be interpreted at the level of the produced K*. For this purpose 
the decomposition of the K* contribution into the final-state particles needs to be considered: 



(7.6) 



In the Dalitz plots. Fig. 20, interference between the K*^ and K*'' is observed. Interference 
occurs in a given final state; the K* has to be expanded into its decay products to appreciate the 
meaning of the interference pattern. The next two lines give the final states produced in the above 
reactions. 
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A pair of final states in which the K* can interfere is underlined in (7.7). These two states 
have K*+K~ and K*"K" intermediate states, the first and the last entries in Eq. (7.6). The two 
annihilation modes marked strong in (7.6) produce the K*K in the form K*+K~ — K*°K° and the 
two weak modes in the form K*+K~ + K*'^K'^. It is not the isospin which drives this dynamical 
selection rule: the two initial states contributing strongly are 

^•^Si ^ K*+K- - K*°K° + c.c. strong , 

(7.8) 



while the two states 



^•^So ^ K*+K- - K*°K° + c.c. strong , 

^'^Si ^ K*+K- + K*°K° + c.c. weak , 

_ (7 9) 

^'^So ^ K*+K-+K*''K°+c.c. weak, 



contribute at most weakly. The charged and neutral K* K combinations are produced strongly with 
a relative minus sign. 

Inspecting Table 29 shows that K2 (1430)K production exhibits a similar selection rule: 

^■^Si ^ K;+K- - KfK° + c.c. strong , 

^•^So ^ K*+K- - KfK° + c.c. strong , 

_ (7 10) 

^•^Si ^ K;+K- + KfK° + c.c. weak , 

^'^So ^ k;+K- + Kf K° + c.c. weak . 

The dominance of one isospin channel for these decay modes may be surprising since the nominal 
mass of the K2 (1430) plus the mass of the K exceed the available energy. Hence rescattering of the 
final-state particles is expected to play a large role. 

Equations (7.6) suggest a pattern of strongly and weakly produced K and K* which depends on 
charge conjugation and isospin of the initial state. This generalisation can be tested in annihilation 
from atomic P-states. The Obelix collaboration has measured the density dependence of KKtt pro- 
duction, see Fig. 18. Also here, a preference for one isospin for KK* production is observed; the 
prevailing isospin channel is shown in Table 30. For the ^Pi initial state, the expectation is fulfilled, 
for the ^Pi and ^P2 not. Aimihilation from P-states seems not to follow the same simple pattern as 
observed in S-state annihilation. 
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In the K+K~7r° final state, a clear signal due to (b production is seen. The final state (f>Tr can be 
produced only from states with negative charge conjugation and with isospin 1 = 1. Thus only the 
^'^Si and ^'^Pi states can contribute. At high density, the (p is strong, but weak at low densities, 
indicating a preference for the ^''^Si initial state. If (f>'s are produced as two-step process, through 
K*K and rescattering of the two kaons into a 0, then K*K should be strong in ^'"^Si (which is 
the case) and weak in ^'^Pi. This finding is supported by the partial wave analysis of the data of 
Fig. 18 where the ^'^Pi couples weaker to K*K than ^'^Pi. This question wiU be addressed again 
in Sec. 8.7. 



7.6 Annihilation into K*K* 

The frequency for pp annihilation into K*K* was determined in the BNL bubble chamber experi- 
ment [47]: 

AF(pp ^ K*+K*-) = (3.3 ± 1.1) X 10^4 , 

_ (7 11) 

AF(pp ^ K* " K* " ) (7.3 ± 1.5) X 10"'' . 

The more recent annihilation frequencies from Crystal Barrel are at variance with these findings. 

The Crystal Barrel collaboration has studied the dynamics of annihilation into two K*, in the 
final states K,Ki7r"7r" [191] and K±Ks TT^TT^ [192]. The former is forced to proceed via the ''Si 
initial state, and is dominated by K*K* production. The two K* are produced with £ = 1 between 
them; the two K* spins add up to 5 = or 2 with about equal amplitudes and opposite phases; the 
contribution from 5 = 1 is small. The absolute annihilation frequencies were not determined. 

In the reaction K='=Kg7r^7r°, K*K* production plays a less significant role. The fractional con- 
tributions in the partial wave analysis can be normalised to the annihilation frequency of the reaction 
channel (given in Table 9) to arrive at the following branching ratios: 

BR(pp(3Si) ^ K*+K*-) ~ 1.3 X 10"^ , 
BR(pp(^Si) ^ K*°K*°) ~ 0.2 X 10"^ , 

(7 12) 

BR(pp(iSo) ^ K*+K*-) ~ 1.7 X 10-3 , 
BR(pp(iSo) ^ K*°K*°) ~ 1.5 X 10-3 . 

Aimihilation from the ^Si state requires £ = 1; spin states S = and S = 2 both contribute. As 
above, the amplitudes are similar in magnitude but opposite in phase. However, the dominance of 
K*+K*- over K* ° K* ° shows that both isospin channels contribute. 

Aimihilation from the ^Sq initial state may proceed with £ = and with one total spin only, 
S = 1. The frequency for K*K* production is larger, and the two charge modes are about equal in 
size. The measured phase between the two amplitudes is —50°. Hence we do not find the dynamical 
selection rules as in the other annihilation modes into strange final-state mesons. 

We note that in pp K*K*, the K* momenta are 285MeV/c and the two K* decay before 
having left the strong interaction volume. K*K* production constitutes only a small fraction of the 
final state and rescattering of the final-state particles presumably has a decisive influence. Hence it 
is not surprising that this reaction does not exhibit any striking dynamical selection rule. 



7.7 Discussion of the dynamical selection rules 

Table 37 summarises the most important selection rules. Weak and strong are meant as relative 
weights; annihilation into 3.2 ( 1320)7r from 7 = 1 is listed as weak even though it is a larger branching 
ratio than the one for KK annihilation. But for a given final state, a strong decay mode has a 
branching ratio exceeding its associate weak decay mode by about one order of magnitude. 

In any reasonable model, the transition rate for annihilation of the pp system to a two-meson 
final state is proportional to 
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Table 37: Summary of the most important dynamical selection rules. A cross x indicates G-parity 
forbidden reactions. 





/ = 


1=1 




7 = 


7=1 




strong 


X 


^Si^ KK 


weak 


strong 




X 


weak 








iSo^ a2(1320)7r 


strong 


X 


3Si^K*K + c.c. 


weak 


strong 


^Si^ a2(1320)7r 


X 


weak 


iSo^ K*K + c.c. 


strong 


weak 


^Po^ TTTT 


strong 


X 


^Si-^ K^(1430)K + c.c. 


weak 


strong 


3Sl^ TTTT 


X 


weak 


iSo^ K|(1430)K + c.c. 


strong 


weak 



1. the probability to find the needed isospin component in the protonium wave function (initial 
state interaction), 

2. the strength of the hadronic transition operator, 

3. the probability that the final state is formed. 

Are the dynamical selection rules an effect of initial state interaction? Annihilation of the ^'^Si 
level to pTT is strong, hence the isoscalar component of the protonium •'Si wave function must be 
enhanced. Annihilation from '^'^Si — » nir is weak, hence the isotriplet component of the wave 
function must be small at this momentum. From ^'^Si and ^'^Si annihilation into open strangeness 
we conclude that the '^Si isoscalar component must be small. A component can hardly be enhanced 
and suppressed at the same time. If the initial state interaction is responsible, there must be a strong 
momentum dependence of this suppression and enhancement. Note, however, that the protonium 
wave function, as calculated in potential models, sometimes exhibit oscillations, see, e.g., [306]. 

Table 38 Usts the dominant isospin component of protonium wave functions, the momentum at 
which it should prevail, and the final state from which the conclusion is drawn. In particular for 
the ■'Si state, an oscillating wave function would result. Translated into a spatial wave function, 
the isoscalar part should dominate at 0.21; 0.26; and 0.43 fm, the isovector one at 0.25 and 0.32 fm. 
Extremely sharp oscillations between isovector and isoscalar components would be needed. 

However, tensor forces induce ^ ''Di and '^'''Di components in the protonium wave function, 
which contribute to annihilation. Maruyama et al. [307], for instance found that in the framework of 
a specific model with planar diagrams, "inclusion of D states solve the pn puzzle." 

The role of final-state interaction in the pw puzzle was underlined by Mull et al. [308]. For the 
■^'^So : ^'^Si ratio they found about 1:5 without final-state interaction and 1:25 when it is added. 
Surprisingly, the effect is not due to the pn interaction of their model being different in these two 
states, but simply to the existence of this interaction. The effect is almost negligible for P-wave 
annihilating into pw. 

Anyhow, the interaction between final state mesons is not well known, and thus one must rely 
on models. In the case of a well pronounced resonance, one could expect all channels, e.g., tttt, KK, 
etc, to feel attractive forces. In other circumstances, channels that are weakly coupled might behave 
differently. We conclude that final state interaction cannot be neglected, but is not demonstrated to 
explain all observed dynamical selection rules. 

Alternatively, we may try to understand the dynamical selection rules as the effect of an hadronic 
transition operator. Then we find that the transition to two isovector mesons is small from the isovec- 
tor component of protonium wave functions (at least for ^Si and ^Sq) and large for the isosinglet 
component. The remaining dynamical selection rules are condensed into the observation that in pp 
aimihilation, kaons (and K*, ..) are produced in the form K+K~ and in K^K*^. Whenever we probe 
the annihilation potential with two isovector mesons in the final state, we find a large isospin 7 = 
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Table 38: Interpretation of the dynamical selection rules as an effect of the initial state interaction. 
The Table gives the prevailing isospin component of protonium at a given momentum if the origin of 
the dynamical selection rules is assigned to the initial state interaction. The reaction from which the 
dominant isospin component is determined is given in the last column. 



Lj (pp) 


dominant isospin 


2«+iL^(pp) 


dominant isospin 


at momentum 


From 


^Si 


7 = 


^Po 


7 = 


928 MeV/c 


7r7r 


^Si 


7=1 


^So 


7=1 


797 MeV/c 


KK + c.c. 


^^Si 


7 = 


^So 


7 = 


773 MeV/c 


pn 


^Si 


7 = 1 


^So 


7 = 


616 MeV/c 


K*K + c.c. 


^Si 


7 = 


^So 


7 = 


460 MeV/c 


a2(1320)7r 




7=1 




7 = 


~ 120 MeV/c 


K^K + c.c. 



contribution. When we probe the annihilation potential with final states with open strangeness, we 
find the combination pp K+K~ — K°K° to be large compared to pp K+K~ + K^K". 

We conclude that the dynamical selection rules can economically be interpreted as suppression 
of specific hadronic transition operators. The hadronic operator acts in a similar way for two pseu- 
doscalar mesons, for a pseudoscalar and a vector and for a pseudoscalar and a tensor meson. The dy- 
namical selection rules drive annihilation preferentially into specific favoured flavour combinations; 
other quantum numbers Uke spin and angular momentum of the final-state mesons are seemingly 
less relevant. 

It is not easy to find a microscopic derivation of such a selection rule. An attempt was made by 
Niskanen and Myhrer [309], who analysed quark-antiquark annihilation (necessary for NN pir) 
into one or two gluons and found, interestingly, that it is suppressed if the NN initial state is ^'^Sq. 
A non-perturbative generalisation would be desirable. 



84 



Annihilation dynamics 



8 Phenomenological analysis 

There is an abundant literature on NN annihilation. This complex process has been studied in a 

variety of models. Most detailed numerical analyses are now obsolete, to the extent that the param- 
eters of the models have been tuned to reproduce early and incomplete sets of data. However, the 
underlying mechanisms still deserve to be presented and compared to other possible mechanisms. 

In this section, we will expand on the following questions: the pion multiplicity and the clustering 
of mesons into resonances; the range of annihilation; the role of initial and final state interaction; the 
importance of symmetries; the probability of producing strange mesons and the mechanisms of OZI 
violation; the interpretation of dynamical selection rules. Ideally, answers to these questions should 
come from a full understanding of hadronic interactions which we do not yet have. 

8.1 Initial state interaction 

8.1.1 Overall suppression 

Before considering annihilation mechanisms, it is important to underline the role of initial-state 
interactions. If a specific process NN ^ mi m2 . . . is calculated in a given model without accounting 
for initial-state interactions, its rate is overestimated by orders of magnitude. In more technical 
words, a Born-approximation treatment is unacceptable, while a distorted-wave Born approximation 
can be rather realistic. In any realistic model, the NN wave function is dramatically suppressed at 
short distances by the cumulative effect of annihilation channels. 

8.1.2 Induced channels 

A serious warning by Green et al. (see, e.g., [310]), is that if initial-state interactions are taken 
seriously with our current ideas on nuclear forces, there are ample transition amphtudes NN <-> 
N*N + c.c, or NN ^ N*N*, with at least one baryon or antibaryon being a spin (N* ~ A), 
orbital or radial excitation of the nucleon. Similarly strange-meson production can proceed via a AA 
doorway or other hyperon-antihyperon intermediate states. Hence the quark content of the initial 
state does not reduce to (uuduud). This might influence the conclusions drawn about the hierarchy 
of various quark diagrams describing annihilation. 

8.1.3 Selective suppression 

Even for ratios of branching ratios, a pure Bom-Oppenheimer treatment of annihilation might be 
misleading, as different initial states do not necessarily experience the same suppression, and in a 
given partial wave, the damping of the wave function depends on the momentum range which is 
explored (see, e.g., Ref. [306], and discussions later in this section). 

The NN interaction is investigated in scattering and protonium experiments, and described in 
a number of models that combine long-range meson exchanges and short-range absorption. This 
subject is discussed in a previous review article [1]. 

As already stressed in Sec. 4.4.5, the long-range potential has a strong spin and isospin de- 
pendence, starting with the one-pion-exchange term which includes a fi.f2 <?i-<?2 operator for the 
central interaction, and something analogous for the tensor one. In potential models, a dramatic 
spin-isospin dependence is induced in the various partial waves contributing to annihilation at low 
energy. 

8.1.4 Orbital mixing 

Meson-exchange models also predict some"orbital" mixing in natural-parity states of protonium, 
due to the tensor component of the NN interaction. For instance, the authors of Ref. [311] combine 
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the so-called "A2" model (to be defined later) with initial-state protonium wave functions estimated 
from potential models, and computed the relative rates for annihilation into two mesons. Some ratios 
are found to be extremely sensitive to details of the NN interaction used to produce the protonium 
wave-function. In particular, neglecting the ^Di admixture into the ■'Si wave function sometimes 
changes the results by an order of magnitude. It was already noted by Green et al. [310, 312] that 
the ^Po model requires an orbital momentum between the annihilating quark and antiquark. For an 
overall S-wave, this is provided by departures from the harmonic-oscillator behaviour of the wave 
function^. For a spin S = 1 system, an overall J = 1 state comes either from an S-wave, or, in 
presence of tensor forces, from a D-wave. 

8.1.5 Isospin mixing of protonium 

In potential models, the effect of initial-state interaction is predicted to be particularly important 
in protonium wave functions, as compared to scattering wave functions relevant for annihilation in 
flight. This was underlined by Kaufmann and Pilkuhn [257J, and several other authors. 

In short, a pure pp state corresponds to equal weights for isospin 7 = and 1=1. However, the 
charge-exchange potential induces transitions from pp to nn, and when the pp and fin ampUtudes 
are recombined into amplitudes of given isospin 7, one component is often much larger than the 
other one. For instance, ^Pq is clearly dominated by its 7 = part. However, this firm prediction 
of meson-exchange models is not confirmed by studying the systematics of branching ratios. See, 
for instance, the discussion in Refs. [167, 306, 313|. It is somewhat of a paradox that the pion- 
exchange force nicely reproduces the pattern of fine splitting of P- levels of protonium [1] as well as 
the hierarchy of hadronic widths (see Table 11), but resists experimental checks for branching ratios. 

The expUcit calculation of Ref . [311] illustrates the concerns about isospin mixing in protonium. 
For instance, the r]p^ (7 = 1) to ir^p^ (7 = 0) is calculated to be much smaller than the experimental 
value. 

In Ref. [314], Gutsche et al. also discuss this question of isospin mixing. They estimate the rate 
for radiative annihilation pp ^ 7 + X, with X being 7r°, r], r]', p or w. The transition is sensitive 
to interference between the 7 = and 7=1 components of the protonium wave function, and thus 
probes the isospin mixing predicted by potential models. A rather good agreement is found with the 
data. However, the rate for 7^ cannot be reproduced. 

8.1.6 Checking isospin mixing in protonium 

The first attempt to deduce the isospin ratios from experiment was made in [313]. The model as- 
sumed that the transition matrix element for annihilation into two mesons is entirely determined by 
the isospins involved. There are transitions from the 7 = initial state to two isoscalars and to two 
isovectors, and from the 7=1 initial state to one isoscalar and an isovector. The transition matrix 
elements are then supposed to be independent of the quantum numbers of the initial state, apart from 
a normaUsation which could be different for different initial states. The results are certainly model- 
dependent, and the errors are large. The isospin-mixing coefficients are found compatible with the 
predictions of potential models, but, due to the large errors, also with the absence of mixing effects. 
The subject was further discussed by Dover et al. [306] and more recently by the Crystal Barrel 
coUaboration [209]. 

8.1.7 Isospin content of pp in pd 

The analysis can be extended to aimihilation on deuterium. To our knowledge, there are no published 
calculation of the detailed isospin content of pp within pd. Note that the isospin content of the 

^For instance, a J = state of three bosons has pure £ = orbital momentum between any two constituents in the 
harmonic-oscillator model, but this is not true for other type for binding interaction. 
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Table 39: Ratio of isovector to isoscalar fraction of the protonium wave function for various initial 
states ^°+^Lj. Tlie tliree theoretical values correspond to different NN potentials, as compiled in 
Ret. [315]. 



Initial Potentials Data analysis 

state KW DRl DR2 Ref. [313] Ref. [209] 





0.68 


0.68 


0.8 


72 

' ^ -0.18 


0.50 


+0.48 
-0.29 


^Si 


1.22 


0.95 


1.26 


1 -17+0.39 
-■^•^ ' -0.28 


1.17 


+0.30 
-0.23 


^Po 


0.03 


0.03 


0.05 


1.16 ±0.34 


0.41 


+0.11 
-0.09 


■^Pi 


9.4 


9.7 


6.5 


9±5 






^Pi 


0.96 


0.82 


0.61 


0.81 ±0.51 







neutral j3p state in pd is not necessarily the same as in protonium, due to the presence of the third 
hadron. 

The transition pd into a isovector meson mi and a isoscalar meson m2 proceeds via the (squared) 
isovector component of the pN subsystem. This component is 1 for pn and smaller for pp. Thus 

AF(pd ^ TT°uj + n) 8.4 ±0.4 „ ^„ , „ 
if- I L = _ n CO n no 

AFipd-^TT-Lu + p) 12.1 ±1.4 "-^^^^^-"^ 

AF(pd^^";; + n) _ ^^''^ 

AF(pd^7r-7? + p) " 8.1±2.0 -"■''^^"•^^ 

The frequencies were obtained by calculating mean values from the numbers given in Table 32. 
Those normalised to the number of pd annihilations have been multiplied by 2. The mean values do 
not include the results from [108]. The n^u! data from [68] and [78] are superseded by a reanalysis 
[276]. 

The reaction pp^ ttuj is dominated by the ^Si initial state; the ratios (8.1) suggest a mild 
deviation from a pure pp system at annihilation for which ratios 0.5 would be expected. According 
to Table 41, annihilation into ttt] goes through the ^Fq isovector component of the pp system which 
is calculated to be very small, see Table 39. 



8.2 Final state interaction 

When discussing annihilation, one can hardly forget the strong interaction of the produced mesons, 
which are likely to rescatter, form resonances, decay into pions, etc. Even lighter mesons, which 
escape faster from the interaction region, might be affected. 

Final-state interactions were introduced, e.g., by the Bonn group [308] and found to play an 
important role for the spin effects associated with NN mr or KK, and for the pn puzzle. We shall 
see later in this section that it is suggested that rescattering could explain the observed deviations 
from the OZI rule. 

Clearly, if a channel is not or only weakly populated by the main mechanism of annihilation, it 
might receive a non-negligible contribution of rescattering from another final state. In a fictitious 
world where pp does not couple to KK, the tttt — s- KK reaction would be crucial. 

On the other hand, rescattering among channels which are copiously produced, presumably has 
Uttle effect on the annihilation frequencies. The loss towards other channels is compensated for by 
the feed-back. More generally, in a regime of strong production, saturation occurs. If one adds 
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another mechanism, the production rate grows less than would naively be expected by adding the 
squared amplitudes estimated independently. 

Consider for instance the following annihilation frequencies for ui and p production: 

^fc^= 0.41 ±0.06, ^f^HP^= 0.23 ±0.06. (8.2) 
AF(pp p°w^) AF(pp 'qui) 

There is no obvious dominance of p production over w production, though the p meson is more 
likely than lo to be formed in the rescattering of primary mesons. 

In short, the final-state interaction is not expected to dominate the systematics of branching 
ratios, but should be kept in mind for channels which are suppressed. 



8.3 Pion multiplicity and two-meson doorway scenario 

An average of five pions are produced in nucleon-antinucleon at rest. In a baryon-exchange picture, 
or in microscopic quark models, it is natural to assume that a few meson resonances are primar- 
ily produced, the observed final states resulting from the decay of these resonances. A large, but 
quantitatively still unknown, fraction of all annihilation modes proceeds even via two-meson inter- 
mediate states, sometimes called quasi two-body annihilation. Insisting on two-body aimihilation as 
the dominant contribution however does not account for the large fraction of events leading to high 
pion multiplicities. 



8.3.1 The Vandermeulen model 

Vandermeulen [316] suggested that these high pion multiplicities could come from an enhanced pro- 
duction of mesons with high mass. He observed that meson resonances are preferentially produced 
with nearly the maximal mass which is allowed by phase space. The effect can be parametrised by as- 
suming that annihilation proceeds via two primary mesons and writing the frequency pp — > mi +m2 
for producing two mesons with masses m\ and m2 as proportional to 



F{q) = q exp -RsoV s - (toi + m2)^ 



(8.3) 



where q is the momentum of the mesons, as given by Eq. (4.2); ,s is the Mandelstam variable which 
is s = Arrip for annihilation at rest. With a reasonable adjustment at i?so = 1-2 GcV~^, Vander- 
meulen was able to account for many key features of annihilation in flight, covering a wide range 
of antiproton momenta. The Crystal Barrel collaboration applied (8.3) to annihilation at rest and 
derived [167] i?so = 0.83GeV~^. The model of Vandermeulen was further developed, e.g., by 
Mundigletal. [317]. 

To get more insight into Eq. (8.3), we define the annihilation amplitude by F{q) = q \ f{q)\'^. For 
the case of two mesons having the same mass, f{q), thanks to Eq. (4.3), reduces to an exponential 
function. 



f{q) = exp 



l-SO / 



(2mi 



= exp[-gi?so] 



This momentum distribution corresponds to the spatial distribution 

Sir) = z z^''''' 



i?2 ' 

J. i-cr^ 



(8.4) 



(8.5) 



of a source having a size of Rgo, interpreted as the size of the annihilation source. 

Equation (8.3) does not take the centrifugal barrier into account. The formula can be extended 
to include a suppression of high angular momenta at small momentum q 



Fi{q) = qfj{q) = qBj{q) exp [-2qR,. 



(8.6) 
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where fg{q) is the transition amphtude and the functions are given by 

B,{q) = l, B,{q) = ^l^, = z=iqR^,f. (8.7) 

The decay momenta are measured in units of 1 /Rsi, where Rgi corresponds to a strong interaction 
radius. In the zero-range approximation, i?^; 0, and the Bf {q) become proportional to q"^^'. 

Equation (8.6) can be read as the product of the two-body phase-space factor q by the probabihty 
of getting the required orbital angular momentum and linear momentum of the outgoing mesons. 

Note that the enhancement of high-mass mesons is qualitatively equivalent to the pion suppres- 
sion that was found necessary by Green and collaborators [310,312]. Schematically, annihilation 
produces qq pairs in some quantum numbers. The overlap of such qq pair with an actual meson is 
reduced in the case of the pion, due to smaller pion size, and its intricate internal structure. 



8.3.2 An illustration: the pseudoscalar mixing angle from pp annihilation. 

Already in 1983, Genz [318] suggested that the "quark hne rule" could be apphed to pp annihilation 
into two pseudoscalar mesons and the branching ratios could be used to determine the pseudoscalar 
mixing angle. Genz was not able to obtain the right result: the pre-LEAR data were of low statis- 
tical significance and had no redundancy, so that dynamical questions hke the influence of phase 
space and orbital angular momentum barrier could not be investigated. The situation improved once 
the Crystal Barrel Collaboration determined a large number of branching ratios involving rj and rj' 
mesons [167]. The Obelix collaboration extended the study to low-energy lip scattering [153]. 

Table 40 summarises the results. It gives the initial state of protonium which contributes most 
significantiy to the final state, the particle against which r]' or r] recoils, the ratio of annihilation 
frequencies, 

^ 26{x,n) AF(NN ^ 7?X) ' ^ ' ^ 



and the ratio 



_ DR(NN-.VX) 
^ 2S(^'V) DR(NN ^7iX)' ^ ■ ^ 



of dynamically corrected branching ratios, defined as 

DR = AF/f>(g), (8.10) 

where -Fi(g) has the simplified form (8.6), even for unequal masses. The factor 2''('^'''^ = 2 for 
annihilation into r\r\ and 2*^'^'''^ = 1 elsewhere accounts for the Bose symmetry of two identical 
bosons in final state. 

As seen in Table 41, annihilation into two pseudoscalar mesons proceeds dominantly via the "^Pq 
state, and we shall restrict ourselves to these data. 

We further assume that the proton and antiproton have no intrinsic s,s component and couple, 
in the notation of Sec. 3.3.1, only to the nn part of the ry and 77' wave functions. This is the so- 
called Zweig or OZI rule; its vahdity in pp annihilation (where it is also called Quark Line Rule and 
abbreviated as QLR) will be discussed later in this section. If the and i]' wave functions are written 
as in Sec. 3 as 

\r] ) = cos(eps - 6id) I nn) - sin(eps - 6id)| ss ) , ^ 
\rl) = sin(eps - 0id)| nn ) cos(©ps - 6id)| ss) , 

the ratio dx is given by dx = l/tan^(0ps — ©id)- 

Although they do not incorporate any phase-space factor, the uncorrected ratios tx are already 
reasonably consistent. This indicates that the dynamical corrections are small, i.e., that the phase- 
space effects are cancelled out by the Vandermeulen factor. This requires the size of the aimihilation 
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Table 40: Ratios rx, and ratios dx corrected with (8.6), of annihilation frequencies into rf and ?? 
mesons recoiling against the same particle X. The dominant initial atomic state is given in the first 
column, the orbital angular momentum £ between the outgoing mesons in the last column. 



Atomic state 


X 


rx 


dx 




^Po 


ttO 


0.50 ±0.10 


0.53 ±0.11 





^Po 


7r+ 


0.63 ±0.16 


0.65 ±0.17 





^Po 


V 


0.62 ±0.07 


0.69 ±0.07 









0.52 ±0.07 


0.78 ±0.11 


1 


^Si 


P 


0.42 ± 0.08 


0.60 ±0.11 


1 



source to be small, compatible with i?so = l/(2"^p)> and the interaction radius to be large, Rsi = 
1.5 fm. The systematic errors are are estimated to be about 10% . 

The mean value of the corrected ratios dx is 0.65 ± 0.07 with a = 3.3 for 4 degrees of 
freedom, leading to 

Gps = -(15.9 ± 1.5)°, (8.12) 

which is not inconsistent with other measurements of this quantity. The choice of a larger Rso or 
smaller R^i leads to larger mixing angles. The importance of the result lies less in the final number 
but rather in giving insight into the dynamics of the annihilation process. 



8.4 Dynamically corrected branching ratios 

The dynamically corrected branching ratios for different two-body final states, as defined in Eq. (8.10) 
are listed in Table 41. The corrections are empirical and take already into account parameters like 
interaction strengths, wave function overlap and finite size effects. Hence these ratios should not be 
compared to results from a full model calculating annihilation dynamics. Rather, the rates represent 
the squares of elementary transition matrix elements. The parameters of the Vandermeulen model 
were fixed using branching ratios with no open strangeness. A suppression of these final states due 
to a penalty for ss pair creation as advocated, e.g., in [5] is not yet accounted for. This aspect will 
be discussed in Sec. 8.6.3. 

The decays from ^Pq states into two pseudoscalar mesons are normalised to 7r°7r°, these from the 
^Si states into vector plus pseudoscalar to p°7r°, and from the ^Sq states into tensor plus pseudoscalar 
to 100% a^7r°. The branching ratios BR are taken from Table 29 and Eqs. (6.13). 

Table 41 provides insight into the hierarchy of annihilation modes: 

• The dynamical selection rules, discussed in Sec. 7, are confirmed after the dynamical cor- 
rection is applied. The pn puzzle manifests itself in the smallness of the branching ra- 
tio from isospin 7=1 initial states into two isovector mesons (3'^ column from bottom) 
compared to those from 1 = 1 initial states (P* column). The branching ratios called 
i;(/pp — > Iijh) = i(l 1,1) here are smaller by about one order of magnitude compared to 
t{0 ^ 1,1)- Ii,2 are the isospins of the two mesons. 

• The sum of branching ratios for kaonic decay modes t{l/2, 1/2) is of the same order of 
magnitude as t(l 1,1). In annihilation from ppS-states, one isospin component dominates; 
the "^Po initial states prefers to decay into K+K^. 

• Aimihilation modes t{0 0, 0) show no systematic behaviour. For comparison, final states 
containing an r] (//') meson, should be corrected for their ss component by multiplying the DR 
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Table 41 : Dynamically corrected branching ratios of selected two-body modes. For easier compari- 
son, the absolute values DR(^Po ^ tt'V") = (133 ± 26) 10"^; DR(^Si n^p") = (32.2 ± 1.7) 10"^; 
and DR(^So -» a2(1320)7r°) = (71 ± 10) 10"^ are normalised to 50 (to account for Bose symmetry), 
100 and 100, respectively. 



Isospin 


PP 


channel 


DR 


PP 


channel 


DR 


PP 


channel 


DR 


^ 1 -n 




TV n 


50 






100 


'So 




100 


0^0-1-0 




VV 


12.7 ± 1.0 




UJT} 


115 ±6 






3±2 


0^0-1-0 




m' 


17.5 ±1.1 




OJT)' 


89 ±12 








0-»0-|-0 












0.44 ±0.11 








1 -» 1-1-0 










■KU) 


41 ±3 




7rf2 


55 ±7 


1^1-^0 






15.6 ± 3.2 




PV 


26 ±3 








1^1 + 




7r?7' 


8.3 ±0.5 




pv' 


15.8 ± 2.8 






1.38 ± 0.23 


1^1 + 










d>n 


4.3 ±0.5 








1/2 + 1/2 




K+K" 


4.7 ±0.4 




(K*K)j=o 


0.9 ±0.9 




(K^K)/=o 


9.2 ±2.8 


1/2 + 1/2 




Ks,iK3,i 


< 1 




(K*K)/=i 


13.7 ±2.2 




(K|K)/=i 


0.8 ± 0.3 


1^1-1-1 




7r"'"7r~ 


3.0 ±0.2 






28.0 ± 4.4 






23.4 ± 4.6 


1/2 + 1/2 




K+K" 


1.6 ±0.2 




(K*K)/=o 


13.6 ± 2.2 




(K^K)/=o 


< 3 


1/2 + 1/2 




K.Ki 


1.4 ±0.2 




(K*K)/=i 


2.5 ±1.0 




(K^K)/=i 


6.5 ± 1.8 



by 1/(0.65 ± 0.07) or 1/(0.35 ± 0.07), respectively. With this correction the pp coupling to 
r/r/ and jyry' is about 1/2 of the tttt coupling, while the coupling to wry and wtj' is twice larger 
than the couphng to pn. The DR for ■qi2 is very small; of course this transition has a small 
phase space and a larger angular momentum barrier (with £ = 2). But for the corrections as 
suggested here, the DR for T/f2 remains small. 
• The annihilation modes t{l 0, 1) are smaller than i(0 ^ 1, 1) by a factor 2 to 4. 

8.5 The size of the annihilation source 

From NN scattering data, several "radii" were determined [1]. The charge-exchange reaction occurs 
at typically i?cc ~ 2.5 fm, the mean strong interaction radius is R^i = 1.5 fm and at Bun = Ifm, 
annihilation takes place. Two further length scales are given by the Compton wave length of the pion 
and proton, 1.4fmandmp^ 0.2 fm, respectively. 

8.5.1 Baryon exchange mechanism 

Baryon exchange was the first mechanism proposed for annihilation, in analogy with electron ex- 
change in e+e~ annihilation (see Fig. 21). This naively suggests a range of the order of l/(2mp) ~ 
0.1 fm, where rUp is the proton mass. With this sole value 0.1 fm for any possible range, size or 
form-factor parameter of a model, one would never reproduce the observed ratio of annihilation to 
elastic pp cross sections, nor the smallness of the charge-exchange cross-section, nor the occurrence 
of P-wave annihilation at rest, such as pp — > tt^tt^. With large form factors associated with baryon 
exchange, annihilation acquires, however, a more realistic spatial extension. 

In principle, baryon exchange is very appealing, since it uses for NN annihilation the same 
baryon-baryon-meson couplings that enter NN forces. It is not clear, however, whether the form 
factors can be extrapolated from the NN scattering region to the NN annihilation region. There is 
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Figure 21 : Positron-electron annihilation mediated by electron exchange: the range of the induced 
absorptive potential is the inverse of twice of the electron mass. 

also a warning by ChristilUn [319] that the exchange of a A or other nucleon resonance leads to a 
range and a strength comparable to those of nucleon exchange. Hence there is an ambiguity on how 
the series of exchanges should be truncated, before adjusting the parameters. 

However, baryon exchange models have been developed by several authors, in particular the 
Bonn-Jiilich group, see [308, 320], and references therein. These authors treated annihilation "adi- 
abatically", first by estimating - in their approach - the transition to a few two-meson channels and 
mimicking the remaining channels by an empirical optical potential. They then gradually increased 
the number of channels explicitly accounted for. In this framework, one can test the role of various 
ingredients such as: the coupling constants in the meson-nucleon-nucleon vertices, the associated 
form factors, the role of final-state interaction, etc. The role of A exchange has also been investigated 
in this framework [321]. 

8.5.2 Annihilation range from statistical considerations 

Another estimate of the range for annihilation was proposed by Fermi [322] . It is reproduced here in 
the notation of Amado et al. [323]. The phase-space for annihilation into n particles with momenta 
Pi and energy Ei, see Eq. (4.10), is given by integrating 



The integrals with different n having different dimensions, one should introduce a scale factor L 
such that the rate for producing n pions reads 



Adjusting the average multiplicity (n) = J2 S Rin) to be 5 gives L ~ 1.2 fm. This simple 

reasoning also reproduced the observed variance An ~ 1 of the multiplicity distribution. A more 
careful calculation accounting for the 2tt factors of Eq. (4. 1 0) would even increase the range _L [323] . 
Obviously, the scale factor L is related to the size of the fireball from which pions are emitted and 

thus to -Rsi- 

8.5.3 Quark rearrangement 

The quark model explains why a baryon and its antiparticle can energetically annihilate by simple 
quark rearrangement, i.e., why the reaction {qqq) + {qqq) {qq) + {qq) + {qq) can occur at rest. A 
reasonable phenomenology of both mesons and baryons is, indeed, achieved [324] if it is assumed 
that the interquark potential obeys the so-called "1/2 rule" 



n 



2Ei 



Qnis^iVi}) = S^{P-pi Pn)Yl 



(8.13) 




(8.14) 




(8.15) 
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This is obtained by exchanging colour-octets. A linear confinement V = Xr for mesons is better 
generalised as a F-shape interaction for baryons [325] 

Vy = Aniin(di + ^2 + ds) , (8.16) 

where di is the distance from the quark i to a junction J, whose location at minimum corresponds 
to the well-known Fermat-Torricelli point of elementary triangle geometry. This genuine 3-body 
interaction is however close to the result for the 1/2 rule, but slightly larger, since [325] 

+ r-23 + rsi) < min(di +d2+ ds) < "^(^12 + ^23 + rsi) . (8.17) 

The variational principle implies that if V^qq > J2 ^qqi'^'ij) 1'^^ ^^^^ 2M(qqq) > 3M{qq) [326]. 
Note that quarks and antiquarks are assumed here to have equal masses. If m{Q) m{q), then the 
inequality can be inverted into M(qqq) + M(QQQ) < 3Af(Qq) [327]. Very heavy antibaryons 
with three units of heavy flavour would not "annihilate" on ordinary matter. 

In the quark model, the finite size of annihilation is understood from the composite structure of 
baryons and mesons. Mesons are produced according to their ability to make a "bridge", i.e., pick 
up a quark in N and an antiquark in N [328-330]. The baryon size governs the spatial spread of the 
final mesons. 

This can seen as follows. A typical transition potential is 

(*f|0|*f) , (8.18) 

where the operator O correspond to various terms of the interaction Hamiltonian. One gets a good 
idea of rearrangement by estimating the simple overlap integral corresponding to O = 1 . Other 
matrix elements have similar shape. For describing the initial state, the individual coordinates fj, 
corresponding to the labelling of Fig. 22, can be rearranged into 




Figure 22: Notation for the rearrangement diagram. 



^ r4 + r5 + r6-ri-r2-rs 
r = ;= 

5 ri -I- r 2 -I- f 3 -I- r4 -I- fs -I- re 
H = -pz 



r2 — ri 

V2 ' 



An 



2r3 — f 1 — n 

V6 



(8.19) 



so that in a harmonic-oscillator model for baryons, the initial wave function reads, in the centre-of- 
mass 

(8.20) 



cxp 



-2(^ + 4 + '^N + A^)J F{r), 



assuming a simple factorisation of the relative wave function F and the internal quark or antiquark 
motion. 
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Similarly, the final state is described by the internal meson coordinates Xi and global meson 
coordinates i/i, i = 1, 2, 3, 

fi+3-fi ^ ri + fi+3 

X, = —j^ , m = , (8.21) 

and out of the latter, on can built R and 

'-^JT' — ' ^ 

so that, again assuming a harmonic oscillator and factorisation, the final-state wave function reads 



9/4 



= ( — 1 exp 



^/-<2 , -•2 I -'2\ 
--(Xl +X2+X3) 



G{p,X). (8.23) 



If one keeps an intermediate set of variables made of R, f, p, A and 

_ _ X2 - xi _ PN - Pj^ _ 2^3 - xi- X2 _ An - 

one can integrate over these latter variables and get 

2„9/2!^15/4 / bf'^\ /a - \ 

= nm^a + b) "^^'^ ^^'^ [-—) (- 2 + ^ • ('-"^^ 

As mentioned earUer, the r dependence is governed by the oscillator parameter b of the meson wave 
functions, while the spatial distribution of the final mesons, described by the coordinates p and A, is 
linked to the baryon size a. 

This expression (8.25) corresponds to a transition NN 3 mesons. If iterated with its conjugate, 
it gives the contribution to the absorptive part of the NN ampUtude, possibly identified as the driving 
term of the imaginary part of the optical potential. One gets an exactly separable potential 

Qm[V] cx exp(-6f V2) exp(-6r"V2) , (8.26) 

acting between an initial NN wave function F{r) and a final mesonic wave function F{r"), in 
contrast with the local character of optical potentials used in current phenomenological pictures. 
Green et al. [310, 331] and Ihle et al. [329] went a Uttle further and studied to what extent quark 
rearrangement can describe a large fraction of the observed annihilation cross-section. The predicted 
order of magnitude is reasonable, at least in the framework of simple constituent models. This means 
that quark rearrangement is hardly negUgible. It even opens up the possibiUty of attempting a first 
study of the systematics of branching ratios [332]. 

However, to account for two-body modes, kaon production and the detailed features of branch- 
ing ratios, other quark diagrams have to be included, with some, if not all, incoming quarks and 
antiquarks annihilating and some quark-antiquark pairs being created in the final state out of the 
released energy. A phenomenology has been developed, to try to extract from the data the relative 
importance of the various types of diagrams. Different authors have reached conflicting conclusions, 
thus illustrating the difficulties of the art of annihilation diagrammatics, which will be presented later 
in this section. 



8.5.4 Annihilation ranges from the Vandermeulen model 

The values i?so = l/(2mp) ~ 0.1 fm for the radius of the annihilation source, and Rgi ^ 1.5fm 
for the average interaction radius are suggested by data on annihilation into 77 and 77'. They deserve 
some comments: 
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• A small source cannot provide C > orbital angular momenta. Without the Blatt-Weisskopf 
correction, the branching ratios would scale as q'^. The momentum for cor]' production is only 
half of the momentum for wr] production. In addition, the r]' has a smaller fraction of fin 
quarks than the rj. Hence rx should be expected to be about 1/10 while experiments give 
about 1/2. The angular momentum barrier cannot be large. This is why the -Rgi in Eq. (8.7) 
must be large. One could simplify into Rgi = Ran = 1 fm but adopting a larger value for Rgi 
gives a slightly better description of the data. 

• A large source cannot provide large linear momenta, as expressed by the term exp(— 2gi?so)- 
Only for small Rso, the momentum dependence is weak enough to describe simultaneously 
data on pp irr] and pp tor]'. 

• Rso is small compared to the radii used to fit cross sections [1]. The smallness of the anni- 
hilation source i?so ^ 0.1 fm underlines the importance of the real part of the interaction: at 
Rsi ^ 1 fm strong interactions lead to a strong attraction which focuses the wave function into 
a region (of size Rgo ~ 0.1 fm) where annihilation takes place, a scenario often underlined by 
Shapiro [6]. 

These results illustrate the hot debate about the annihilation range. Shapiro [333] and others insisted 
that the annihilation range must be in the order of magnitude of the Compton length of the annihilat- 
ing baryons, independently of how the annihilating objects are constructed from their constituents. 
On the other hand, the annihilation range needs to be properly defined. Nucleons are composite par- 
ticles, as well as mesons. As soon as there is sufficient overlap of the wave functions, rearrangement 
of quarks can occur and mesons are produced. Thus annihilation sets in at large distances [334]. 

8.6 Quark diagrams 

Figure 23 shows a sample of annihilation diagrams where the flavour flow is represented by lines. 
A2, A3, R2, R3 are abbreviations to denote quark diagrams with two or three mesons or mesons 
resonances produced, with or without crossing the lines. These diagrams are not Feynman diagrams. 




R3 R2 



Figure 23: Annihilation (A2, A3) and rearrangement (R2,R3) diagrams for NN annihilation. 

there is no mathematical prescription as to how to calculate annihilation branching ratios from these 
pictures. They have to be supplemented by a model providing initial and final state wave-functions, 
and an operator describing pair creation. So the question arises whether these diagrams can be a 



Phenomenology 



95 



useful guide to annihilation processes. The main argument in favour of quark diagrams lies in the 
dynamical selection rules and the observation, at the end of Sec. 7.7, that the dynamical selection 
rules are related to the flavour flow. 

Diagrams with the same topology might be thought to be equivalent, for instance to those in 
Fig. 24 where either three mesons are created and rescatter into a resonance or two mesons are pro- 
duced and one meson undergoes a subsequent decay. In actual model calculations, with constituent 
wave functions and empirical creation/annihilation operators, the two diagrams are not necessarily 
equivalent. 



Figure 24: Possible meclianisms for production. 



8.6.1 Quark-antiquark creation or annihilation 

Hadron physics requires more than the "naive quark model", an oversimplified approximation to 
QCD, where the number of dressed constituents is frozen, as in ordinary quantum mechanics. A 
minimal extension has been proposed [335-338] to account for the decay of meson and baryon 
resonances. Usually, a "qq pair-creation operator" is introduced in an ad-hoc way, its strength is 
fitted to reproduce a first resonance width, and then the model is (rather successfully) applied to 
predict all other widths. The calculation of the hadronic widths for, say, A ^ B + C, contains an 
overlap integral involving this creation operator, and the wave functions of A, B and C computed in 
a specific quark model. For a review, see the book by Le Yaouanc et al. [336]. 

In the physics of decaying resonances, the momenta are rather low, so it is not too much of a 
surprise that a single ansatz can account for the regularities of the observed widths. The current 
prescription is that qq is created with vacuum quantum numbers, i.e., in a '^Po partial wave. This 
hypothesis is not tested in great detail, since in a decay A ~* B + C , there are not too many 
possibilities for the angular momentum between B and C in the final state. It is thus an audacious 
enterprise to use the same type of model for describing annihilation: the momentum of the emitted 
mesons is much larger; the qq creation or annihilation operator is required to work more than once; 
much more freedom is imaginable for the angular momenta, etc. Quark models of annihilation 
mostly use the ■'Pq model or a ^Si model, where the qq pair has the quantum number of an isoscalar 
vector, or a combination of both. See, e.g., the discussion by Mandrup et al. [339, 340], Niskanen 
and Myhrer [309], Dover and Fishbane [341], and Maruyama et al. [307]. 



8.6.2 Planar and non-planar diagrams 

A large fraction of pp annihilation events, called two-body or quasi-two-body annihilation, produce 
two primary mesons. These mesons may be unstable and decay with a short lifetime. The reaction 
pp pTT with p decaying into tttt is a typical example. There has been much discussion on whether 
pTT production is better described by the diagram A2 in Fig. 23, without crossing of quark Unes, or 
by R2, which avoid annihilation and recreation of one quark-antiquark pair but requires crossing of 
quark lines. 
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Besides this specific example, attempts have been made to single out a few dominant diagrams, 
leading to definite predictions for the hierarchy of branching ratios [318]. Guidance was sought from 
two main sources. 

• The OZI rule, as discussed in Sec. 3.3.4 suggests that disconnected diagrams are suppressed. 
A further step is to assume that only diagrams with planar topology (diagrams that can be 
drawn on a plane without intersecting lines) are dominant. 

• The suppression of disconnected diagrams is unanimously acknowledged to become rigorous 
in the large Nc limit of QCD, where A'c is the number of colours. However, the question 
of the suppression of connected but non planar diagrams (such as rearrangement) has been 
reanalysed by Pirner [342]. His conclusions do not support the claim that planar diagrams 
dominate. 

Models assuming dominance of annihilation diagrams over rearrangement diagrams predict the 
transition amplitudes T for annihilation into two isovector and into two isoscalar mesons to be 
the same. For instance, in these models, the rate for pp cut] is similar to pp pn, once 
some corrections have been applied. One correction is required for the ss component in the rj wave 
function. Another correction is the Vandermeulen factor Fg{q) (see Eq. 8.6). The dynamically 
corrected ratios DR are given in Table 42, with Q = 0ps — Oid + 7r/2 and 0ps = —15.9° (see 
Eq. 8.12). 

The ratios in Table 42 are largely incompatible with 1, whilst some have the right order of 
magnitude. With this warning. Table 42 could support the hypothesis that planar diagrams might in 
some cases drive the leading contribution to two-meson annihilation. 

Table 42: Ratios of dynamically corrected annihilation frequencies for pp annihilation into two 
isoscalar and two isovector mesons. 



1 DR(pp^r?a;) ^ i.77±o.lO DR(p_p ^ p p ) ^ o.04 ± 0.04 
cose DR(pp pV") DR(pp ^ ww) 

cos^Q ^T^i?'-^^"^? = 1.67 ±0.14 \ P^^^''''^ ^ f < 0.05 (la) 

DR(pp^w) cos^ e DR(pp ^ 7r"a2(1320)) ^ ' 



8.6.3 Strangeness production 

In models where A2 is small compared to R2, a small fraction of events contain a pair of kaons. On 
the other hand, if A2 is the leading mechanism of annihilation, and if SU(3) symmetry is approxi- 
mately valid, a very large number of kaons is expected. This is a serious problem for models based 
on A2, since the overall yield of strange particles is only about 5%. 

This is why, in models where A2 and other planar diagrams are assumed to dominate annihila- 
tion, it is crucial to introduce an explicit SU(3) breaking, in the form of a suppression factor A for 
ss-creation, compared to fin. A value A = 1 corresponds to the SU(3) limit. To push the fraction 
of hidden strangeness production close to the experimental value, models based on planar diagrams 
need a value as low as A ^ 0.1. 

A theoretical foundation for ss suppression was given by Dosch and Gromes [343], who showed 
that ss could be very much suppressed just above threshold, by a kind of tunnelling effect. In 
the hadronisation following high-energy reactions, the production of strange quarks is also reduced 
compared to the production of up and down quarks. This is why fragmentation models incorporate 
a strangeness suppression factor A. It increases with energy from values as low as 0.1 at a few 
GeV to ^ 0.3 at 30 GeV [344]. j3p annihilation being a soft process, a value A 0.1 seems not 
unreasonable. Such a low value has decisive consequences for the interpretation of annihilation data. 
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If A ^ 0.1 were true, SU(3) symmetry would be dramatically broken. The suppression would result 
in heavier kaon resonances being more suppressed than hghter kaons, and also KK dramatically 
increasing in annihilation in flight as compared to annihilation at rest. Hence there is a possibility of 
understanding the data by the dominance of annihilation diagrams supplemented by ss suppression. 

Fragmentation of hadrons is, however, a rather indirect way to address the question of ss sup- 
pression at low energies. Meson decays may be a better guide. The tensor meson a2(1320) has 
dominant decay modes into pn, rjir and into KK. SU(3) relates the decays into rjir and into KK. 
Their ratio, and the decay branching ratios of all tensor mesons, are fully compatible with only small 
SU(3) symmetry breaking. The authors of Ref. [249] fitted 16 decay modes of tensor mesons with 
SU(3) amplitudes allowing for ss suppression. They found A = 0.8 ± 0.2: data on tensor meson 
decays are compatible with SU(3) and rule out a substantial SU(3) sytmnetry breaking. 

Table 43: Selected J/?/> decays probing SU(3) symmetry. Data are taken from Ref. [237]. The 
experimental value is shown, and then corrected for a phase-space factor where p is the 

momentum, and I the angular momentum in the final state. The units for BR/p^^+^ are irrelevant. For 
isospin multiplets, we assume normal weights for each charge state or average over the data, and 
display this average, contrary to [237] where, e.g., pir means the sum of the three channels. 



Channel 


BR (10^3) 


BR/p2^+i 


Channel 


BR (10"-'') 


BR/p2^+i 


PP 


2.12 ±0.10 


1.7 


pw 


4.2 ±0.5 


1.7 


nn 


2.2 ±0.4 


1.5 


ojrj 


1.58 ±0.16 


0.6 


AA 


1.30 ±0.12 


1.3 


K*K + c.c. 


2.3 ±0.2 


0.9 


EE 


1.27±0.17 


1.3 


TTTT 


0.15 ±0.02 


0.04 




0.9 ±0.2 


1.1 


KK 


0.19 ±0.03 


0.06 


AA 


1.10 ±0.29 


1.2 


a2(1320)/3 


3.6 ±0.7 


3.2 


E*E* 


0.52 ±0.07 


0.8 


f2(1270V 


4.3 ±0.6 


3.6 








K^K 


3.4 ± 1.3 


2.4 



The Mark III Collaboration studied J/ip decays into a vector and a pseudoscalar meson in order 
to find "inert" or gluonic components in the r] and r]' wave functions. The result was negative. One 
of the (rather numerous) parameters was the suppression of ss pair creation compared to the creation 
of uu or dd pairs. From their fit a value A = 0.8 was deduced [345]. 

We looked for strangeness suppression in J/ip decays into baryons, into pp, nn, AA, EE, SS, 
AA, and S*S*. After correcting for phase space (i.e., after division by the respective decay mo- 
menta) the squared invariant couplings are similar in size, but scale with 0.75"= where Ug is the 
number of ss pairs created. The creation of a second ss pair in SS is certainly a non-perturbative 
process but still governed by A ^ 0.75! 

In short, the physics of tensor mesons and charmonium decay does not support ss suppression, 
and hence calls, in the case of NN annihilation, for an important role of non-planar, rearrangement 
diagrams. 

Perhaps strangeness production in NN cannot be explained in a simple uniform manner and 
should instead be examined for each type of final state. Table 44 compares (dynamically corrected) 
branching ratios for annihilation into two strange mesons with those for annihilation into two isovec- 
tor mesons. The data are grouped into processes in which the two vector mesons come from isoscalar 
(D-coupling) and isovector (i^-coupling) protonium states. D or F coupling belong to standard 
SU(3) notation [237]. 
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The ratios with D- and _F-coupling differ remarkably while being internally compatible. There 
is an astonishing consistency between these ratios when processes of the same SU(3) structure are 
compared. Again, the flavour flow is responsible for gross features of two-body annihilation pro- 
cesses. 

The ratios with F-coupling suggest that A2 might be the relevant quark line diagram, possi- 
bly with a 'penalty factor' (~ 0.5) for producing a ss pair. Then the rearrangement diagram R2 
must contribute strongly to the production of two isovector mesons from isoscalar initial states. To 
conclude, it seems that the data on strangeness production support a scenario where rearrangement 
diagrams are important. 

The problem of strangeness production is not restricted to the question of the relative impor- 
tance of A2 vs. R2 diagrams. Another mechanism has been proposed by Ellis et al. [346], where 
strange quarks and antiquarks are extracted from the nucleon or antinucleon sea, leading to specific 
signatures, for instance for the spin effects in pp AA, and for the violation of the OZI rule (on 
which more in the next subsection). Holinde et al. [347] suggested a possible role of initial-state 
interaction, via NN — > YY, to describe the pp (fxj) annihilation in flight. The mechanism can 
also be applied to a violation of the OZI rule, since the hyperons Y and Y can produce a (j) together 
with light mesons by simple rearrangement. 

Table 44: Dynamically corrected branching ratios for annihilation into strange mesons and into two 
isovector mesons in the final state for different SU{3) flavour couplings. A "-|-c.c." is implied for every 
final state with a K* or Kj. 



D-coupling 



F-coupling 



DR ("Po 



KK) 

DR (''^Po tttt) 

DR (^Si K*K) 
DR(''-'Si ~ 



■ pi") 



DRp^Si ^ a2(1320)7r) 



= 0.047 ± 0.004 



0.049 ± 0.008 



DR('Su ^ KK2(1430)) 



= 0.033 ± 0.010 



DRf Si ^ KK) 
DR (^'^Si ^ ti-tt) 

DR ('So ^ K'K) 
DR (-'''So ^ ptt) 



DR 



DR^^-'Si ^ a2(1320)7r) 



= 1.00 ± 0.20 



0.58 ± 0.13 



(«Si ^KK^(1430)) _ 



34 ± 0.14 



8.7 Violation of the OZI rule 

The quark-line rule, or OZI rule, has already been introduced for mesons, see sec. 3.3.4. To first 
order, the rule forbids production of (f) mesons from initial systems like pp with up and down quarks 
only. The vector meson mixing angle, Gy = 39° for the quadratic GMO mass formula, allows for 
the smaU cj) 7r"'"7r~7r° decay rate, and also cj) production from pp is permitted with an expected 
ratio 

(vector, expected) = tan2(ey - Oia) = 0.004 , (8.27) 
where (vector) is defined as the ratio of the dynamicaUy corrected annihilation frequencies 

while the rx (vector) are the corresponding ratios without dynamical corrections. The ratios rx and 
dx are defined analogously for the tensor mesons £2 (1270) and £2 (1525). 

Different analyses of bubble chamber data gave results in the range from 0.07 to 0.23 for the 
ratio of annihilation frequencies AF(pA'^ — > tt^) to AF(pA'^ — > iru) [348]. Based on a selection of 
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results, a mean value 0.085 ± 0.030 was given in [124]. Obviously production is much stronger 
than the anticipated ratio 0.004 suggests. 

At LEAR, (p production was studied by the Asterix [124], Obelix [132, 147, 148, 154], and Crystal 
Barrel collaborations [170, 176]. In Table 45 we give ratios rx(noiict) and (ix(nonct) for X = tt^ 
and X — rj from the two initial states '^Si and ^Pi. 

Table 45: Ratios of (fi/uj and f2(1525)/f2(1270) production. Given are the pp initial state, the recoil 
meson, the ratio without and with dynamical corrections, and the ratio expected from the nonet mixing 
angle. In the last column, the isospin component is given which gives the dominant contribution to 
the corresponding final states with kaons. 



PP 

3,lp^ 
l,lp^ 



X rx (vector) 

TT° 0.099±0.012 
V 

V 



0.005±0.009 
0.10±0.07 



dx (vector) 
0.105±0.013 
0.0038±0.0010 
0.005±0.009 
0.12±0.09 



FromGMO K*K 

0.004 J = 1 

0.004 / = 1 

0.004 / = 

0.004 J = 



PP 


X 


(tensor) 


dx (tensor) 


From GMO 


K^K 


3,lQ 

So 




0.020±0.004 


0.025±0.005 


0.012 


J = 


PP 


X 


rx (vector) 


dx (vector) 


From GMO 


K*K* 


3.1Q 

SO 


P 


0.018±0.007 


0.030±0.010 


0.004 


both 




to 


0.0095±0.0035 


0.017±0.006 


0.004 


both 




1 


0.25±0.09 


0.026±0.09 


0.004 


both 



The ratios vary over a wide range: the pp annihilation frequencies from the "^Si state to (pr], from 
the ^Pi state to 07r, and the frequency to f2(1525)7r are all of the order of magnitude expected from 
the meson nonet mixing angles. A few processes give a moderately large production rate like (pp 
and (puj. For some reactions however, the OZI violation is really large, the (p/u) ratio being about 
10% or larger. 

Three interpretations of this large excess of (p production compared to lo production have been 
pursued. Figure 25 sketches the three scenarios. 




Figure 25: Diagrams which could contribute to (p production: as shake-off of hidden ss pairs in the 
nucleon wave function, via rescattering of Kaons from secondary decays, and from formation of four- 
quark exotic states. 

Dover and Fishbane [348] link the excess in (p production to the production of four-quark exotics 
with hidden strangeness (qsqs), e.g., to the tail of the C(1480) meson [349]^. 

^The C(1480) was observed as </)7r resonance and interpreted as (nsns) resonance. However, it was never confirmed. 
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The nucleon wave function is known to contain an ss component. This is evident from deep 
inelastic scattering [350] or from the so-called 7r-nucleon a term (Tttn [351]. There is no guarantee 
however that virtual qq pairs may be shaken off in an OZI rule violating diagram. If this is assumed, 
the sharp selectivity of the processes leading to large OZI rule violating effects can potentially be 
understood [346, 352] as originating from the ss component being negatively polarized. 

The PSI group studied OZI rule violating effects in NN annihilation in a series of papers [353- 
355]. See, also, [356-358]. The large ncf) annihilation frequencies were interpreted by rescattering 
into (j) mesons of KK pairs from KK* + K*K annihilation and or pn rescattering from p'^p~. The 
large rate for 7^ channel was understood as the effect of pcf) and oj^ production, and vector meson 
dominance. The study was extended in [359] to include (t)'jT~^iT~ production in pp annihilation at rest 
and in flight. 

We notice that the ratios rx and dx are large when the (p and u) are produced from the isospin 
component which gives a large contribution to kaonic final states. Rescattering of, e.g., KK in the 
K*K final state is certainly proportional to the frequency with which K*K with the correct isospin 
is produced. Indeed, the OZI violation is found to be correlated with the K*K production strength 
of the required isospin. This observation supports the rescattering interpretation of the strong OZI 
rule violation. As emphasized in [352], there are more processes in which OZI rule violating effects 
can be studied and possibly linked to a hidden ss component in the NN wave function. 

Decays of the J/tp can help to elucidate the problem further A sample is given in Table 46. The 
J/ip wave function has little hidden strangeness, and any tt^ production must be due to rescattering. 
In reversing the argument, the absence of in J/tp decays could possibly be interpreted as evidence 
that the large OZI rule violating in NN annihilation must come from the nucleon wave function. 

Table 46: Selected decays modes of J/V^ 

J/lp — S- rjLU TTLO TTcj) K*K 

Rate (10-4) 15.8 ±1.6 4.2 ± 0.6 <0.068 92 ± 6 

From the r]u}/ncj ratio, isospin violation can be estimated to occur at the 14% level. The K*K 
final state was not analysed to identify isospin breaking effects; however a K*K ( J = 1) contribution 
of more than lO"'^ can be expected. Rescattering may then lead to ncj) production at a detectable 
level. 

The problem was studied in [360] using a dispersion theoretical approach. The interpretation 
of large 07r frequency in pp annihilation as an effect of rescattering was confirmed. The low (pn 
production in J/tp decays is due to its larger phase space, and is compatible with the rescattering 
mechanism for OZI rule violating effects. 

8.8 Flavour flow and flavour symmetry 

Support for the use of quark-line diagrams can be found by comparing branching ratios for annihi- 
lation into specific two-body final states. In Table 47, we compare dynamically corrected annihila- 
tion frequencies DR, using Eq. (8.6). The numerators correspond to final states with one isovector 
and one isoscalar meson, and the denominators to annihilation into two isovector or two isoscalar 
mesons. 

The ratios require some corrections. In the process pp — !■ nrj, the initial pp system couples only 
to the nn component and not to the ss component. The OZI rule reduces the amplitude for this 
process by cos 9, and the branching ratio by cos^ 6. Correspondingly, pp ttt]' is reduced by 
sin^ 9, where 9 is defined in Eq.(3.7). We use cos^ 9 = 0.65 ± 0.07 and sin^ 9 = 0.35 ± 0.07. 

The table reveals a surprise: most results are about compatible with each other, except for those 
containing the DR for ivrj' (which are too large). The mean value of those ratios gives 0.29 ± 0.03. 
with a /Np = 2. There is no a priori reason why these ratios should be similar. The two final 
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Table 47: Ratios of dynamically corrected annihilation frequencies DR into an isovector and an 
isoscalar meson compared to production of two isoscalar or two isovector mesons. Bose symme- 
try is taken into account. 
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states TT^r] and 7r°7r° are produced from the ^Pq initial state, prj and lot] from "^Si, £2(1270)77" and 
a2(1320)7r° from ^Sq. Table 47 includes final states with two Ught mesons like tt and rj, or with two 
massive mesons like u! and p; the ratio is formed using pseudoscalar, vector and tensor mesons in 
arbitrary combinations. In some cases, broad mesons are only in the numerator, in others only in the 
denominator. Obviously, the flavour content and the flavour coupling are the decisive ingredients. 
Table 47 supports the conjecture made in Sec. 7.7 that the flavour flow has a decisive impact on 
annihilation dynamics. 

In Table 47, the two annihilation modes pp — > p^p^ and pp T/f2(1270) are excluded. The 
frequency for annihilation into pp p'^p'^ is compatible with zero, and there is no known reason for 
this effect. The reaction pp 77f2(1270) needs i = 2 and is very close to threshold; the measured 
frequency is compatible with zero but is also expected to be small. 

8.8.1 SU(3): quark-line rule, s-chaimel resonances and baryon exchange 

As we have seen pp annihilation can be discussed in rather different languages. Quark models 
describe annihilation in terms of planar and non-planar diagrams (often called annihilation and rear- 
rangement diagrams). On the other hand, pp annihilation may prefer to proceed via a few s-channel 
resonances (e.g., by mixing between the pp system and (qqqq) states close in mass). In this case, 
a description in terms of s-channel amplitudes may be more appropriate. Or, alternatively, baryons 
and mesons might be the relevant degrees of freedom, and pp annihilation could be most efficiently 
described by baryon exchange amplitudes. Figure 26 visuaUses the different approaches. The ques- 
tion to be addressed is which scheme is best suited to incorporate the most important aspects of 
annihilation, and in particular, whether dynamical selection rules find a natural interpretation in one 
of the three coupUng schemes. The three descriptions are related by unitary matrices which were 
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Figure 26: Annihilation of protons and antiprotons in different coupling schemes. From left to right: 
the quark line coupling depicting the flavour flux, the s-channel coupling with meson formation in the 
s-channel, and baryon exchange coupling where baryons are exchanged in the t-channel. 

developed in [361]. 

Figure 26 shows for each couphng scheme two basic diagrams. Their meaning is illustrated 
using the quark line coupling scheme. The two outgoing mesons can be symmetric or antisymmetric 
with respect to their exchange leading to a doubling of diagrams. In any diagram, one antiquark 
and one quark reach the final state as 'spectators'. The spectators can be either {u, u) or (d, d), 
again leading to a doubUng of diagrams*. Hence 8 diagrams are needed to describe the annihilation 
process. These are shown in figure 27. In case of two mesons in the final state belonging to the 
same multiple! (e.g. for annihilation into two pseudoscalar mesons), the generalised Pauli principle 
requires R3 to vanish. 

For a quantitative analysis we use the dynamically corrected branching ratios of Table 4 1 . These 
should be related to the squared SU(3) transition matrix elements. In the fits we use a strangeness 
suppression factor A = 0.7 to 0.8. The results do not depend critically on this assumption. A SU(3) 
error of 20% is introduced. The latter is added quadratically to the experimental errors of Table 41. 

AH three coupling schemes allow us to fit the data with identical y^. There are 10 measured 
branching ratios for annihilation into two pseudoscalar mesons and 1 1 for annihilation into a vector 
and a pseudoscalar meson, 7 amplitudes in the former case and 8 in the latter one. Since the pseu- 
doscalar mixing angle has been fixed from these data (Eq. 8.12), the number of degrees of freedom 
is 2 in both cases. These two sets of branching ratios are fitted with — 1.3 For annihilation 
into a pseudoscalar and a tensor meson, there are only 8 data points and 8 amplitudes; the data are 
reproduced with = 0. The three couphng schemes give identical descriptions of the data. Even 
within one coupling scheme, different solutions exist with identical or similar x^. Hence care has to 
be taken in interpreting the results. 

8.8.2 The quark coupling scheme 

Figure 27 shows the decomposition of the quark couphng scheme. There are 8 amplitudes which 
may be fit to the values given in Table 48, where a positive sign has been chosen for the amplitudes 
A^. The largest contribution is given by R^. Before discussing the meaning and significance 
of this result, a comment is made concerning the amphtudes A^ and A^. They have opposite 
signs and are about equal in amplitude. The approximate relations Aj*" + A^ hold for all five 
annihilation processes from initial S-states. For annihilation into two pseudoscalar meson, the data 
are also approximately compatible with this relation with an alternative solution where |A+| ^ 
I AJ I . Setting Af + A^ = substantially helps further discussion: annihilation into two isovector 
mesons (first line in Table 41) is given by the SU(3) matrix element |2Af + 2A^ + Rt?- 

'^More precisely, the spectator either belongs or not to an antisymmetric (u, d) pair. 
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Figure 27: Explicit representation of quarl<-line amplitudes of all independent couplings in the quark- 
line scheme. Pairs of (anti-)quarks antisymmetric with respect to their exchange are connected by 
small loops. 



Table 48: Best-fit values for the amplitudes in the quark coupling scheme. 
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R+ 
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-2.1 


9.7 


-0.3 




1.4 


-1.5 


1.2 


PS-V 


2.3 


-1.9 


10.9 


1.2 


-3.9 


1.1 


-1.8 


-1.6 


PS-T 


1.3 


-1.8 


8.9 


-2.8 


-1.0 


0.8 


-1.5 


-1.8 



can be chosen positive or negative; for + ^ two values are found; both are large but 
differ in the precise number. Only for A^ + A^ = 0, does R^,^ have a unique value. As long as 
the interest is in a quahtative understanding, it suffices to state that rearrangement diagrams give a 
substantial contribution to annihilation and should not be neglected. This statement is compatible 
with the findings presented above (see Sec. 8.5.3). 

Among the rearrangement diagrams, plays the dominant role. This result has to be inter- 
preted with some precaution. Annihilation into pw from the ^Si initial state is driven by ^R^", 
annihilation into tor] by ^Rf + Rj. Obviously, there is a small and a large negative value of R2' 
which can satisfy the relation. So we can conclude only that the data are compatible with Rj and 
R3 both being small but RJ could be large. The relation + Af ~ reproduces the dynamical 
selection rules governing strangeness production. The pn puzzle and the generahsation to nn and 
a2(1320) is related to the dominance of R^. 
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8.8.3 The s-chaimel coupling scheme 

The 8 SU(3) amplitudes can be decomposed into .s-ciiannel couplings denoted as si , sg^^ , , sSas ' 
S8„a > ^10) SfO) S27. The coupling scheme classifies the SU(3) structure of the intermediate state in 
the s-channel. The coupUngs of the pp system and of the two mesons to an octet intermediate state 
(qq) or (qqqq) can be symmetric or antisymmetric; this flexibility leads to four amplitudes. The 
intermediate state could also be flavour singlet, decuplet, {10} or a {27}-plet. 

As stated above, the fit with all amplitudes leads to identical results. But the s-channel couphng 
scheme leads to a different interpretation. The four annihilation diagrams go through a qq interme- 
diate state; obviously no s-channel decuplet, {10} or {27}-plet can contribute to these diagrams. 
More interesting is the observation that the s-channel decuplet, { 10} or {27}-plet amplitudes can all 
be set to zero without a significant deterioration of the fit. This fact is due to the dominance of Rf 
which can be decomposed into si , ss,, , S8^„ , Ss^^ , S8„„ only. If pp annihilation dynamics is inter- 
preted in the s-channel coupling scheme, there are no "exotic" components even though four-quark 
qqqq "crypto-exotic" components do play a very significant role. The (generalised) pn puzzle is 
thus interpreted as absence of exotic states in the s channel of pp annihilation. The strange pattern 
in strangeness production finds however no straightforward interpretation in the s-channel coupling 
scheme. 

8.8.4 The t-u-chaimel coupling scheme 

The 8 amplitudes of the i-u-channel coupling scheme are given by Bi, Bs^^, Bs^^, Bfi^^, B^^^, 
Bio, BfQ, B27 where Bs^^ e.g. stands for an octet baryon in the t- or w-channel with symmetric 
N-N-meson coupling to both mesons. These 8 amplitudes fit data again with a as given above; 
suspiciously, singlet baryon exchange plays the most important role. 

In the baryon-exchange picture, contributions due to {10} and {27}-plet exchanges could be 
expected to be small; further, the four N-N-PS and N-N-V couplings (with different symmetries ss, 
as, sa, ss) should be related by one F/D ratio for the N-N-PS coupling and one F/ D ratio for the N- 
N-V coupling. Predictions for these ratios are model dependent. SU(6) predicts F/D = 2/3 [362] 
while an analysis of hyperon decays yields F/D = 0.575 ± 0.016 [363]. The N-N-V coupling is 
less well established and considered as free parameter here. 

The data are completely incompatible with these constraints. Relaxing the F/ D ratio yields un- 
reasonably small F/D and still a bad fit. Also, there is no link of the amplitudes with the dynamical 
selection rules. In this simple form baryon exchange does not provide additional insight. Likely, 
the exchange of excited baryons would be needed to achieve a better understanding of annihilation 
dynamics using meson and baryons as fundamental actors. But this would be in conflict with the 
spirit of this study in which simple interpretations of the branching ratios are searched for 

Details of the fit method and results based on older data can be found in a thesis at Mainz [364]. 
The results are very similar to the ones obtained here. 

8.8.5 Discussion 

A similar range of energy release is involved in pp annihilation and m i/ip decays, and a compari- 
son between these two processes is instructive. Data on annihilation into a vector and a pseudoscalar 
meson are collected in Table 49. The branching ratios for J/i/) are taken from [237], for pp anni- 
hilation from Table 29 in Sec. 6.5. We use only pp annihilation from the ^Si initial state since its 
J^^ quantum numbers are identical to those of the J/V'. The branching ratio for pp p°7r° from 
file ^Si state is 21.2 x 10~^, whilst from J/V^ decay it is about 5 times smaller, 4.23 x 10~^. This is 
probably due to the larger phase-space available for J /ijj which gives access to more final states. To 
facihtate the comparison, the branching ratios for p°7r° are normalised to 100. 

Production of strange mesons is much larger in J/^/; decays than in j3p annihilation. The rate in 
J/^ is of the magnitude expected by simple SU(3) considerations. It is anomalously low in pp. This 
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Table 49: Annihilation fractions for pp annifiilatlon and brancfiing rations for J/ip decays into vector 
and pseudoscalar mesons. 

p°n^ ujr) ujr)' <pr] (j>r]' K^+K" K*'^W 

J/tp ICQ 37.3 ±3.8 3.9 ±0.6 15.4 ±1.7 7.8 ±0.9 59.1 ± 4.7 49.6 ±4.7 
pp 100 109.0 ±5.2 57.1 ±7.6 0.38 ±0.10 x 7.0 ±1.2 7.0 ±1.2 



reduction has often been interpreted as being due to the dramatic suppression of ss pair creation in 
low energy physics. In Sec. 8.6.3, the reduction was suggested to arise from large contributions from 
rearrangement processes in which no new quarks need to be created. These processes are, of course, 
absent inJ/tp decays where all quarks in the final state have to originate from the vacuum. This fact 
enhances cj) production in 3/ip decays when it recoils against the rjovq' . In pp annihilation, there is 
no or little ss in the wave function, and </> production is small. 

In J/'^ decay, the branching ratios decrease when the mass of the produced mesons increases, by 
the ordinary phase-space effect. A very interesting feature of pp annihilation is that this reduction is 
not observed. In annihilation dynamics, production of large masses is preferred over production of 
high momenta, as noted by Vandermeulen [316], and discussed in Sec. 8.3.1. In atomic physics, a 
similar effect is observed in the Auger effect, with low momenta highly preferred. This is understood 
by a better overlap of the wave function, and this overcomes phase-space considerations. Perhaps 
a derivation of the Vandermeulen effect is to be sought in the quark wave function of higher mass 
mesons, as compared to those of light mesons. 

Annihilation into two vector mesons was excluded in the SU(3) analysis. In Sec. 7.6 it was 
shown that the daughter mesons from the two vector mesons can interact before leaving the strong 
interaction region. This may be the reason that the dynamical selection rules are not observed in 
this case. In particular the strong ujuj channel, in absence of a strong counterpart, requires 
a large contribution of decuplet or 27-plet four-quark configurations in the intermediate state. It 
is interesting to note that in the 77 pp channel, isotensor interactions are required above the 
pp threshold [365] and this may be the reason why 'exotic' exchanges are realistic. The ratio of 
77 ^ P^P^ [366] and 77 — > p'^ p~ [367] measured at LEP by L3 support [368] the previous 
findings of the TASSO collaboration. 
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9 Conclusions and outlook 

In this review we have discussed how a nucleon and an antinucleon annihilate to create mesons in 
the final state. In NN annihilation, part of the incoming mass, 2 GeV for annihilation at rest, is 
transformed into pure energy that hadronises into mesons. The fraction presumably increases from 
a vanishing value for pure rearrangement to reach 1 or 2 GeV for processes involving annihilation 
of quark-antiquark pairs in the initial state. The experiments performed at LEAR have boosted our 
knowledge about mesons formed in the annihilation process; inclusive and exclusive final states are 
often known with good accuracy. The multiplicity distribution of charged and neutral mesons and 
their momentum distribution finds an interpretation in a thermodynamic picture. A fireball with 
120 MeV temperature and size 1/m^ annihilates into pions. The number of charged and neutral 
pions also follow a statistical distribution. This simple picture misses, however, the importance of 
meson resonances and their production in two-body annihilation modes. These reactions are very 
important for an understanding of the annihilation process. 

A systematic approach to annihilation dynamics requires annihilation frequencies to be measured 
for a large number of reactions. We now have 'complete' information on frequencies for annihilation 
into two pseudoscalar mesons, one pseudoscalar and one vector or tensor meson, and into two vector 
mesons; the information is complete in the sense that the frequencies of all kinematically allowed 
annihilation modes are known. Unfortunately, this is only true for the pp initial state; data on pn 
annihilation are still obscured by the spectator proton in pd annihilation, and data on low-energy np 
annihilation into exclusive final states are still scarce. 

A given two-body annihilation channel can be produced from different atomic states or partial 
waves of the pp system. The assignment of a fraction of an annihilation frequency to a specific 
partial wave requires further experimental input. This could come from experiments using polarized 
antiprotons and protons; such experiments were not carried out at LEAR. The initial atomic states 
can be restricted when events are tagged by observation of a coincident X-ray from the pp atomic 
cascade. In some cases a series of frequency measurement at different target pressures is sufficient 
to constrain the quantum numbers of the initial state from which annihilation occurred. 

Using these techniques, annihilation frequencies were recalculated to yield branching ratios for 
pp annihilation from a specified initial atomic state to two-meson final states. These are the numbers 
which should be compared to models of NN annihilation. 

Strong interaction physics, which is the domain of Quantum ChromoDynamics, has evolved 
considerably in recent years. However, precise predictions dealing with NN annihilation remain out 
of reach. Models need to be developed to identify the relevant degrees of freedom and the effective 
forces and symmetries. In an attempt to identify the leading mechanisms in annihilation dynamics, 
the experimental branching ratios were corrected dynamically using a model proposed by Vander- 
meulen and extended here to account for centrifugal barrier effects. The corrections parametrise the 
finite size of the pp source, and the preference for annihilation into mesons with a high mass. The 
parameters of the model were determined by the requirement that the pseudoscalar mixing angle be 
reproduced from branching ratios for annihilation into ?/ and rf mesons. 

The dynamically corrected annihilation branching ratios provide a surprise: there are annihilation 
modes which seem to be suppressed in comparison with other modes even though they are perfectly 
legitimate and compatible with all known selection rules. Here we quote Carl Dover [5]: 

The search for signatures of quark- gluon dynamics in NN annihilation is somewhat analogous 
to the search of the phase transition from a hadron gas to a quark- gluon plasma in relativistic ion 
transitions. The signal must be isolated from a background of statistical processes characteristic of 
a system with many degrees of freedom. For the NN system, an important role is played by quasi- 
two-body intermediate states, or "doorway states", which display directly the selection rules arising 
from baryon exchange or quark dynamics. ... These dynamical selection rules provide key signatures 
of the annihilation mechanism. 

It has been shown that the dynamical selection rules can be grouped into two classes of observa- 
tions. The first selection rule is observed in pp annihilation into two isovector mesons (tttt, ptt and 
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a2(1320)7r). These frequencies are large from initial states which allow annihilation into the two 
neutral (non-strange) mesons (i.e., into tt^tt^ etc), and they are small from initial states forbidding 
two neutral mesons by charge conjugation. The selection rule is known, at least for annihilation into 
a pseudoscalar and a vector meson, as the pn puzzle. It has found different dynamical interpreta- 
tions which however missed the generalisation to nn and a2(1320)7r. The analysis of the flavour 
flow in quark line diagrams linked the first class of selection rules to the suppression of decuplet or 
27-plet four-quark intermediate states. According to this interpretation, exotic (or non-octet) four- 
quark states could have a large mass. Annihilation modes requiring these intermediate states are 
then suppressed. Non-exotic four-quark configurations are not only allowed as intermediate states 
but provide even the leading contribution to annihilation dynamics. The large branching ratios for 
pw, iOT], a2(1320)7r and other meson pairs from isoscalar initial states are driven by this mecha- 
nism. This observation does not claim the existence of four-quark resonances in this mass range but 
states only that four-quark configurations can be formed at the NN mass if, and only if, they carry 
flavour-octet quantum numbers. 

The second selection rule is found in annihilation into two strange mesons. At least in annihi- 
lation from S-wave orbitals one isospin in the initial NN state makes a dominant contribution. In 
an analysis of the flavour flow, the selection rule can be traced to a symmetry property between 
quark line annihilation diagrams. Likely, the symmetry pattern reflects the symmetry properties of 
the quark pair annihilation/creation operator. It is often advocated that qq pairs couple to the gluon 
fields with vacuum quantum numbers; this is the basis of the ^Pq model. The annihilation and recre- 
ation of a qq pair with zero total angular momentum may also be governed by instanton-induced 
forces. However, the origin of the symmetry responsible for the second dynamical selection rule is 
so far unexplored. 

There have been attempts to apply the modem concepts of strong-interaction physics to NN 
annihilation, at an exploratory level; however the first results obtained are interesting [369, 370]. 

Lattice simulations, no longer restricted to the quenched approximation, have made dramatic 
progress. However, systems of two interacting hadrons are at the edge of current possibilities, as 
demonstrated by the somewhat contradictory results obtained by different groups on the pentaquark. 
For a review, see, e.g., [371]. 

The method of effective Lagrangians, which has a wide range of applications, is particularly 
suited for low-energy strong interactions, a domain where QCD can hardly be applied directly. In- 
stead, a Lagrangian having the appropriate symmetries, in particular chiral symmetry, can describe 
hadronic processes with only a few parameters. The field was at first restricted to mesons, and evalu- 
ating the TTTT scattering length was a typical challenge. The method is now applied to meson-baryon 
interaction and to nuclear forces. Recent reviews can be found in Refs. [372-374]. 

QCD itself is sometimes considered in a variant with a large number of colour degrees of 
freedom, instead of the actual = 3. It has been shown that the A^c ^ 3c limit provides a simpler 
picture, in which just a few diagrams contribute, the other being suppressed by powers of N^. For 
instance, many features of the charmed baryons are understood from considerations based on 1 /N^. 
It has also been emphasised [375, 376] that a well-understood large and a well-controlled chiral 
theory gives compatible results. Unfortunately, it is still difficult to apply large- A^c methods to NN 
annihilation. There are debates about the hierarchy of diagrams and the selection rules suggested by 
the Nc ^ oo view of annihilation. 

It is fair to say that the aim of most annihilation experiments was to study the spectrum of mesons 
and discover new meson resonances, rather than being solely aimed at the study of strong interaction 
dynamics in annihilation. Instead, annihilation dynamics was a side product of the experiments. The 
main objective was meson spectroscopy and the search for new forms of hadronic matter. These 
were quasinuclear and four-quark states when LEAR was started. Later, the fashion changed to 
glueballs or hybrids. 

The study of the meson spectrum has been rather successful, both in early bubble-chamber ex- 
periments and in recent LEAR experiments. A review on light-meson physics, and in particular the 
implications of LEAR results, is found, e.g., in [377]. For mesons that were already known, the 
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role of LEAR experiments was twofold: their existence was confirmed and their structure and decay 
properties were studied. But also a large number of mesons was discovered in NN annihilation. 
LEAR has provided evidence for several new meson resonances, and in some sectors such as the 
scalar sector J^'^ = 0^'^, there may be too many mesons to be accommodated by qq states, even 
when radial excitations are included. Some of the best experts on the physics of mesons have also 
proposed a multiquark interpretation of the excess mesons^. 

When this review was started, the authors envisaged some apologetic words to mention baryo- 
nium as the main motivation for building the LEAR facility. See Sec. 1 and Ref. [1]. Ironically, this 
review was finished at the end of 2004, and almost every day there is a paper pubUshed suggest- 
ing that one of the new mesons, such as X(3872), D,5(2632), etc., or new exotic baryons such as 
^+(1540) are of multiquark or hadron-hadron nature. See, e.g., [378-381]. 

Throughout this review, a comparison has sometimes been attempted between pp annihilation 
and J/ijj decay. The similarities and differences certainly deserve to be studied more closely. The 
study of NN annihilation does not aim at remaining an isolated field. Topics such as the topology 
of quark diagrams, the rate of hidden strangeness, the production of high-mass resonances occur in 
several processes; in particular the decay of particles containing heavy quarks. It is hoped that an 
unified picture of hadronisation will emerge when analysing the results collected at future antiproton 
facilities and heavy quark factories. 



'We benefited from many discussions with Lucien Montanet on this subject. 
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